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iThe  Individual  emissions  were  assumed  to  be  in  a  coherent,  chaotic,  or  n 

/ 

state.  The  statistical  nature  of  the  emission  times  results  in  fluctuations 
of  the  relative  contributions  of  different  emissions  at  a  given  time.  In  the 
limit  when  the  emissions  overlap  strongly,  the  field  exhibits  the  correlation 
properties  of  chaotic  light,  regardless  of  the  statistics  of  the  individual 
emissions.  This  is  an  important  result.  The  effect  of  Bernoulli  random 
deletion  and  additive  Independent  Poisson  noise  on  photon-counting  statistics 
has  been  examined.  Under  the  action  of  such  deletion  and/or  noise,  both 
super-Poissonlan  and  sub-Poissonian  distributions  move  toward  the  Poisson 
distribution  but  never  convert  from  one  form  to  the  other.  We  also  studied 
the  counting  statistics  for  stationary  and  nonstationary  cascaded  Poisson 
processes,  obtaining  a  simple  equation  for  the  variance-to-mean  ratio 
(noisiness)  in  the  limit  of  long  counting  times .(^-These  results  will  be 
Important  for  understanding  light  detectors,  such  as  avalanche  photodiodes. 

Signal-dependent  noise  is  manifested  in  many  physical  processes  and  often 
represents  a  fundamental  limitation  to  the  minimum  uncertainty  in  the 
magnitude  of  a  physical  quantity.  When  this  quantity  is  the  number  of  photons 
incident  on  a  photodetector,  the  uncertainty  in  the  number  of  photons  is 
proportional  to  the  total  energy  received.  This  signal-dependent  photon  noise 
places  a  limit  on  the  minimum  light  energy  that  can  be  received  with  a  given 
reliability.  When  the  quantity  of  interest  is  the  number  of  impurities 
implanted  in  a  semiconductor  material,  the  uncertainty  in  the  number  of 
impurities  la  proportional  to  their  average  number.  This  signal-dependent 
"Impurity-noise'*  places  a  limit,  in  certain  circumstances,  on  the  minimum 
device  dimensions  in  a  VLSI  system  that  will  function  with  a  given  reliability. 

Substantial  progress  has  been  made  in  understanding  the  role  of 
signal-dependent  noise  in  the  performance  of  source  modulation  formats  in 
fiber  optic  coamunlcations  and  on  the  performance  of  VLSI  systems  containing 
ion-implanted  bipolar  transistors. 

The  recombination  of  excess  minority  carriers  at  grain  boundaries,  or 
other  Interfaces  with  space-charge  regions,  has  been  treated  theoretically  for 
general  energy  distributions  of  Interface  states  (recombination  centers) *  The 
distinction  is  made  between  minority  carrier  recombination  velocity  at  the 
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Block  20  continued  -  Abstract 

(grain- boundary)  Interface  itself,  and  the  effective  recombination  velocity 
for  the  collection  of  these  carriers  by  the  adjacent  space-charge  region. 
Calculations  of  the  effective  recombination  velocity  have  been  made  as  a 
function  of  the  excess  minority-carrier  concentration  at  the  edge  of  the 
space-charge  region  since  this  is  the  quantity  of  most  convenience  for  device 
modelling. 

Results  of  theoretical  and  experimental  investigations  of  the  effects  of 
high  temperature  (27°C  to  300°C)  on  N  and  P  channel  MOS  transistors  have  been 
obtained.  At  the  device  level,  the  temperature  dependences  of  the  electrical 
parameters  are  characterized  individually;  they  Include  the  threshold  voltage, 
the  surface  mobility,  and  the  junction  leakage  currents.  Current-voltage 
characteristics  were  obtained  for  the  subthreshold,  nonsaturation,  and 
saturation  regions  of  operation.  Transistor  performance  limitations  were 
determined  from  the  observed  degradation  mechanisms.  The  conventional  large 
and  small  signal  models  of  the  MOSFET  were  modified  so  as  to  embody  the 
effects  of  high  temperature.  The  static  and  switching  high  temperature 
performances  of  a  CMOS  inverter,  a  basic  logic  cell,  were  also  obtained.  Some 
implications  of  the  experimental  result  on  high  temperature  analog  circuit 
design  are  discussed. 

A  general  theory  of  phonon-assisted  transition  of  majority  carriers 
across  a  grain  boundary  (GB)  has  been  formulated  with  the  GB  angle  as  a 
variable.  The  bulk  relaxation  time  in  the  crystalline  region  and  the 
transition  time  at  the  GB  have  been  found  to  be  the  Important  parameters  that 
determine  the  transmission  coefficient  of  carriers  across  the  GB.  Using  the 
proposed  model,  an  n-type  polycrystalline  silicon  is  employed  as  an  example  in 
a  numerical  calculation. 

Using  the  classical  diffusion  equation  for  an  amorphous  (a)-Sl-Al-Sl 
sandwich  structure,  the  diffusivity  of  a-Si  in  A1  between  470-S70°C  has  been 
measured  and  the  mass-transfer  coefficient  under  the  equilibrium  solid-phase 
epitaxial  growth  condition  has  been  derived.  Our  data  can  be  used  to  explain 
recent  results  of  junction  formation  by  solid-phase  epitaxy.  The  activation 
energy  of  this  process  is  found  to  be  0.80  eV. 


Block  20  continued  -  Abstract 


Thin  films  have  been  grown  on  <111>  germanium  by  thermal  nltrldatlon. 

The  surface  morphology,  composition,  thickness,  and  refractive  Index  of  the 
films  have  been  investigated.  Compositional  analysis  was  performed  by 
electron  beam  x-ray  mlcroanalysls  which  Indicates  that  the  films  contain 
oxygen  as  well  as  nitrogen  and  have  a  nitrogen  to  oxygen  ratio  on  the  order  of 
2.5:1.  Film  thickness  and  refractive  Index  were  determined  by  ellipsometry. 
The  thicknesses  ranged  from  60  to  150  X,  and  the  refractive  indices  ranged 
from  1.5  to  2.3  depending  on  growth  parameters.  A  universal  chart  for 
evaluating  the  thickness  and  refractive  index  of  dielectric  films  on  germanium 
from  the  elllpsometrlc  parameters  ^  and  A  is  also  presented. 

An  experimental  technique  has  been  developed  to  study  the  electrical 
properties  of  semiconductor  grain  boundaries  (GB's)  by  a  focused  laser  beam. 
The  laser  beam  is  trained  on  a  GB  while  the  photoconductivity  of  the  sample  is 
measured.  This  technique  allows  a  localized  region  of  an  individual  GB  in 
semiconductors  with  multiple  grains  to  be  examined.  The  measurement  of  the 
steady  state  and  transient  signals  as  a  function  of  temperature  determines  the 
grain  boundary  barrier  height,  trap  energy,  and  capture  cross  section. 

A  study  of  the  UV  photochemistry  of  organometalllc  molecules  in  the  vapor 
and  adsorbed  phases  has  been  made.  The  gas-phase  spectra  is  Interpreted  on 
the  basis  of  molecular-orbital  theory  and  experimental  data.  The  absorption 
spectra  of  the  surface  adlayers  are  determined  by  computer  subtraction  of  the 
gas-phase  spectra  and  an  interpretation  is  provided. 

Submicrometer  gratings  have  been  etched  in  GaAs  and  CdS  crystals  which 
have  been  immersed  in  an  oxidizing  etch  and  illuminated  with  interfering  laser 
beams.  A  resolution  of  130  nm  was  obtained.  The  use  of  ln-situ  optical 
measurements  of  grating  parameters  allows  ready  optimization  of  the  grating 
fabrication  process. 

Both  pulsed  exclmer  lasers  and  focused  cw  ion  lasers  have  been  used  to 
form  SiO,  patterns  from  spun-on  silicates  on  silicon  substrates.  The  patterns 
formed  using  the  Ar  laser  have  -1-2  urn  resolution.  This  resolution  allows 
direct  patterning  of  SIO2  on  various  substrates  for  both  microelectronics  and 
integrated  optics. 
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Block  20  continued  -  Abstract 

The  deposition  of  thin  films  with  UV  laser  light  and  the  physical 
processes  which  contribute  to  the  deposit  alcrostructure  are  being 
investigated.  TEM  studies  of  the  early  stages  of  fila  growth  show  the 
iaportance  of  the  light  Interaction  with  the  deposit  nuclei. 

Observation  of  the  laser  enhancement  of  the  etching  rate  in  a  planar 
plasma  etcher  has  been  made.  The  enhancement  has  been  aeasured  as  a  function 
of  laser  intensity. 

Photon-assisted  dry  etching  of  GaAs  has  been  achieved.  The  process  is 
based  on  photodissociation  of  halogen-containing  organic  molecules,  using  UV 
laser  radiation. 

The  characteristics  of  the  capillaritron,  a  high  Intensity  source  of  ions 
and  radicals,  are  being  Investigated.  Mass  spectrometry  and  optical  detection 
suggest  that  many  species  are  present  in  the  beam.  Time  of  flight  studies 
Indicate  that  while  the  ions  exhibit  a  wide  spread  in  velocity,  the  neutral 
species  have  a  much  narrower  distribution. 

Work  has  been  completed  on  the  study  of  the  photophysics  of  isolated 

in  rare  gas  matrices.  Following  excitation  with  a  dye  laser,  emission  was 

3+3  3 

seen  from  the  three  low  lying  electronic  states:  B  irO  ,  A  ir,  ,  and  A*  it*  . 

°  u*  lu*  2u 

The  last  state  emits  only  when  the  host  matrix  is  Xe.  The  results  are 
discussed  in  terms  of  the  Importance  of  the  solvent  cage  effect  and  potential 
Implication  for  the  development  of  Iodine  atom  chemical  lasers. 

Low  temperature  matrices  of  1^  in  CCl^  have  been  studied.  Iodine  was 
photodlssoclated  by  excitation  with  a  dye  laser  at  500-540  nm.;  the  fluores¬ 
cence  obtained  subsequently  was  characterized  temporally  as  well  as  spec¬ 
trally.  Spectral  resolution  of  the  fast  component  of  emission  was  achieved, 
but  its  origin  could  not  be  determined.  Comparison  with  the  absorption 
spectrum  of  iodine  in  liquid  CCl^  is  being  sought  to  elucidate  the  nature  of 
the  emitting  species. 

The  rates  of  lntermolecular  vibrational  energy  transfer  among  the  Vj 
bending  fundamental,  the  2v2  overtones,  and  the  symmetric  stretch  funda¬ 
mental  of  C(>2  have  been  measured.  The  system  was  treated  as  a  4-level  system 
and  the  3  measured  eigenvalues  were: 
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The  ladder  climbing  rate  for  the  mode  corresponded  to  the  first  measured 
eigenvalue  of  390  +  80  msec  %orr  ,  in  good  agreement  with  previously 

2 

measured  values.  The  fundamental  appears  to  couple  directly  to  the  (02  0) 
level  and  only  indirectly  to  the  Fermi  coupled  lower  laser  levels,  the  (10%) 
and  (02%)  states.  Work  is  in  progress  to  establish  the  remaining  two 
eigenvalues  to  better  accuracy  and  to  establish  the  kinetic  rates  constants 
derived  from  these  eigenvalues . 

T-V  excitation  transfer  by  hot  atoms  has  been  studied  in  the  HBr-N^O  and 
HC1-CH^  systems  where  channels  for  both  Inelastic  and  reactive  scattering  are 
available.  In  the  HBr-N^O  system  at  least,  both  kinds  of  scattering  are 
observed.  Also,  mode-specific  excitation  of  N^0  by  hot  hydrogen  atoms  appears 
to  be  a  distinct  possibility. 

Emission  from  the  vibrational  mode  of  the  ground  electronic  state  of 
S02  has  been  observed  after  excitation  to  low-lying  singlet  electronic  states. 
Preliminary  results  Indicate  that  the  mode  is  initially  excited,  followed 
by  equilibration  with  the  mode. 

Collision  Induced  multiple-photon  absorption  by  SF  of  cw  C0_  laser 

O  4 

radiation  has  been  carried  out  inside  a  transparent,  capillary  nozzle  beam 

source.  Resonant  laser  excitation  of  the  gas  just  upstream  of  the  nozzle  exit 

enhances  the  internal  energy  by  0.38  eV  molec  1  over  its  value  at  T  -  290K. 

nz 

In  comparison,  simple  thermal  heating  of  the  nozzle  to  544K,  which  provides  a 

velocity  distribution  in  the  beam  similar  to  that  found  with  laser  excitation, 

produces  a  beam  enhanced  in  Internal  energy  by  only  0.03  eV  molec-^.  A 

relative  (laser  on/laser  off)  cross  section  has  been  determined  for  crossed- 

beam  reactions  of  SFt  with  K- .  This  chemi-lonization  cross  section,  o  ,, 

6  2  *  rel* 

decreases  with  vibrational  excitation. 
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The  dominant  pathway  In  the  photodissociation  of  an  anthracene 
endoperoxlde  yielding  excited  singlet  oxygen  plus  an  aromatic  fragment  has 
been  found  to  Involve  an  intermediate  chemical  structure.  The  decay  time  (45 
ps)  of  the  Intermediate,  believed  to  be  either  a  Dewar  like  isomer  or  an 
excited  biradical,  was  observed  by  a  picosecond  laser  induced  fluorescence 
method . 

The  dynamics  of  generation  of  diphenylcarbene  from  a  dlphenyldlazomethane 
precursor  and  energy  relaxation  in  the  diphenylcarbene  fragment  was  carried 
out  using  picosecond  fluorescence  methods.  Information  on  the  dynamics  of 
singlet-triplet  interconversion,  of  key  importance  in  the  chemistry  of 
carbenes,  and  the  energy  separation  between  these  states  was  obtained  from 
these  studies. 

The  stabilization  of  a  twisted  internal  charge  transfer  state  in  the 
molecule  p-dimethylamino  benzonitrile,  which  is  responsible  for  a  dual 
fluorescence,  was  discovered  to  be  due  to  a  local  Interaction  with  polar 
solvent  molecules.  The  dynamics  of  the  twisting  process  were  not  found  to 
correlate  with  the  dielectric  relaxation  times  of  the  solvent  as  had  been 
proposed  by  various  research  groups. 

The  Billiard  Ball  Echo  Model  was  extended  to  include  the  effects  of  using 
sub-Doppler-llnewidth  lasers  in  photon  echo  experiments.  The  resulting  Ellip¬ 
tical  Billiard  Ball  Echo  Model  develops  the  theory  in  a  more  rigorous  manner 
and  explains  such  echo  phenomena  in  gases  as  echo  duration  and  degradations  of 
echo  intensity  with  excitation  beam  angling.  Photon  echo  experiments  are 
described  supporting  the  Elliptical  Billiard  Ball  Model.  These  experiments 
have  detected  echoes  over  a  dynamic  range  exceeding  10^  produced  by  atoms 
excited  as  much  as  15  lifetimes  earlier. 

Two-photon  tri-level  echoes  have  been  generated  which  enable  the  measure¬ 
ment  of  colllsional  relaxation  of  D-P  superposition  atom  states.  Tri-level 
echo  studies  of  S-S  and  S-D  superposition  state  relaxation  in  Li  vapor  has 
begun.  A  search  for  the  first  detection  of  tri-level  echoes  in  solids  is 
underway . 
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I.  QUANTUM  DETECTION  AND  SENSING  OF  RADIATION 

A.  NOISE  IN  THE  GENERATION,  PARTITION,  AND  DETECTION  OP  LIGHT* 

(B.E.A.  Saleh,  D.  Stoler,  P.R.  Prucnal,  K.  Matsuo,  M.C.  Telch) 
(JSEP  work  unit  1,  1982-1985) 

(Principal  Investigator:  M.  C.  Telch  (212)  280-3117) 


1.  Noise  in  the  generation  of  light 

Since  1956,  when  Hanbury-Brown  and  Twiss  observed  correlation  in  the 

fluctuations  of  two  photoelectric  currents  induced  by  thermal  light^^  and  by 
(2) 

starlight,  the  coherence  properties  of  optical  fields  have  been  studied 
intensively,  from  both  a  theoretical  and  an  experimental  point  of  view.^3^ 

The  usual  kinds  of  light  that  have  been  Investigated  are  chaotic  (thermal) 
light,  coherent  (laser)  light,  and  mixtures  of  both.^^  More  recently, 

the  fluctuation  properties  of  antibunched  light  have  received  considerable 
attention.'7^9* 

Fluctuations  in  the  overall  number  of  active  radiators  in  a  source  of 
light  can  be  an  important  determinant  of  its  coherence  properties,  as  pointed 
out  by  Forrester, and  discussed  by  Loudon. This  effect  is  of 
central  Importance  for  scattered  light,  where  the  number  of  active  radiators 
is  a  stochastic  quantity. Such  fluctuations  also  play  a  role  in  the 
generation  of  antibunched  resonance  fluorescence,  as  discussed  by  Carmichael 
et  al./l2),<13)  Jakeman  et  al.,(14)  and  Mandel  et  al.(9),(l5) 

We  have  examined  the  coherence  and  fluctuation  properties  of  optical 
fields  when  the  times  of  emission  of  the  Individual  radiators  are  describable 
by  a  homogeneous  Poisson  process .  In  particular,  we  obtained  the  first- 

and  second-order  field  correlation  functions  for  such  light  in  the  framework 
of  semlclassical  theory  and  quantum  electrodynamics.  The  mean  and  variance  of 
the  photon  count  were  also  obtained.  In  certain  limits,  expressions  for  the 
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photon-counting  distributions  were  derived.  We  considered  in  detail  a  number 
of  special  cases,  Including  individual  atomic  emissions  modeled  by  coherent, 
chaotic,  and  number-state  descriptions. 

Our  model  is  expected  to  play  an  important  role  in  a  number  of  diverse 
applications.  Consider,  for  example,  the  light  emitted  when  a  stream  of 
energetic  electrons  impinges  on  a  betaluminescent  material.  If  the  times  at 
which  the  electrons  strike  the  material  occur  in  accordance  with  the  Poisson 
process,  this  randomness  will  be  imparted  to  the  emitted  optical  field.  This 
effect  will  be  particularly  evident  when  the  electron  current  is  low  (its 
fluctuations  are  then  the  largest).  These  fluctuations  are  in  addition  to 
those  intrinsic  to  the  optical  field. v  ' 

In  previous  work,^^  we  presented  a  semiclass  leal  analysis  of  phenomena 
similar  to  those  described  above,  under  the  assumption  that  the  irradlance 
could  be  modeled  as  shot  noise.  We  have  now  provided  a  more  complete  analysis 
of  this  class  of  processes  when  the  optical  field,  rather  than  the  Irradlance, 
is  modeled  as  a  shot-noise  process.  Thus,  intrinsic  field  as  well  as 
intensity  fluctuations  are  Incorporated  into  our  model. 

Our  results  lndlcate^^^  that  Che  underlying  Poisson  randomness  of  the 
emission  times  imparts  additional  fluctuations  to  the  radiated  field.  When 
the  optical  field  at  a  given  time  is  a  result  of  contributions  from  a  large 
average  number  of  radiators,  the  light  becomes  chaotic,  whether  the  individual 
emissions  are  themselves  coherent  or  chaotic.  This  is  a  consequence  of  the 
central  limit  theorem.  When  the  average  number  of  radiators  contributing  to 
the  optical  field  at  a  certain  time  is  not  large,  the  deviations  from  chaotic 
behavior  have  been  calculated.  These  deviations  are  characterized  by  an 
increase  in  the  normalized  second-order  correlation  function,  which 
corresponds  to  an  Increase  in  the  variance  of  the  number  of  photons  counted  in 
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» 

a  fixed  time  Interval,  and  to  additional  photon  bunching. 

t 

The  excess  variance  Is  proportional  to  the  mean  number  of  photons, 
indicating  that  these  excess  fluctuations  are  particle-like  in  nature. 
Furthermore,  we  have  shown  that  the  corresponding  excess  fluctuations  of  the 
number  of  photons  in  a  given  time  interval  are  enhanced  by  an  increase  of  the 
time  Interval.  This  is  unlike  the  excess  noise  due  to  wave  fluctuations, 
which  are  known  to  be  averaged  by  an  increase  in  the  counting  time. 

It  has  also  been  shown  that,  when  the  lifetime  of  the  individual 
emissions  is  so  short,  or  when  their  rate  is  so  low,  that  overlap  is  unlikely, 
no  interemission  interference  takes  place,  and  the  light  Intensity  is 
described  by  a  shot-noise  stochastic  process.  The  photons  are  then  described 
by  the  shot-noise-driven  doubly  stochastic  Poisson  point  process,  which  also 
exhibits  excess  bunching  of  a  particle-like  nature.  If  counted  over  a 
counting  time  longer  than  the  lifetime  of  a  single  emission,  the  excess 
fluctuations  due  to  random  emission  times  exhibit  themselves  fully.  In  this 
case,  when  each  single  emission  is  coherent,  and  when  it  is  chaotic,  the 
photons  are  described  by  the  Neyman  type-A  and  the  generalized  Polya-Aeppli 
counting  distributions,  respectively. 

We  have  also  formulated  a  general  quantum-electrodynamic  model  for  an 
optical  field  generated  by  a  sequence  of  emissions  at  random  times. This 
enabled  us  to  examine  the  effect  of  random  emission  times  in  cases  when  the 
emissions  cannot  be  described  classically.  For  example,  when  the  individual 
emissions  are  in  the  highly  antibunched  Fock  state,  we  find  that  the 
randomness  of  the  emission  times  results  in  particle-like  bunching  which 
^  quickly  overpowers  the  inherent  quantum-mechanical  antibunching,  as  the  rate 

of  emissions  Increases  and  the  emissions  overlap.  In  the  limit  of  a  large 
number  of  overlapping  emissions,  the  second-order  correlation  function 
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eventually  exhibits  the  usual  (bunched)  chaotic  behavior.  The  Important  role 

played  by  the  randomness  In  the  number  of  radiators  for  the  observation  of 

(14) 

antibunched  light  In  resonance  fluorescence  is  well  recognized. 

Vhile  the  analysis  we  carried  out  has  been  limited  to  the  more  common 
statistical  models  which  individual  emissions  may  satisfy  (coherent,  chaotic, 
coherent-chaotic  mixture),  other  statistical  models  may  be  analyzed  through 
the  use  of  the  same  methods. 

Spatial  effects  have  also  not  yet  been  considered.  These  effects  may  be 
included  by  assigning  positions  to  the  radiators,  and  by  making  the  radiated 
field  a  function  of  position  as  well  as  time.  By  assuming  that  the  positions 
of  the  radiators  in  the  source  volume,  and  their  emission  times,  are  random 
points  in  four-dimensional  time-position  space,  in  accordance  with  a  4D 
Poisson  point  process,  we  may  proceed  to  determine  the  temporal  and  spatial 
coherence  properties  of  the  radiated  field  through  generalizations  of  the 
methods  we  employed. 

2.  Noise  in  the  partition  of  light 

Point  processes  often  undergo  random  deletion  or  partition.  An  obvious 

example  of  Importance  in  optics  is  the  usual  case  of  photodetection,  in  which 

the  quantum  efficiency  is  invariably  less  than  unity.  Another  case  in  point 

is  optical  absorption.  It  has  long  been  known  that  the  Poisson  process,  which 

is  probably  the  most  ubiquitous  of  all  point  processes,  remains  Poisson  under 

(18) 

the  action  of  such  Bernoulli  selection.  ' 

He  have  explicitly  considered  the  effects  of  random  deletion  on  the 
photon-counting  statistics  for  a  number  of  cases  that  appear  repeatedly  in 
optics. The  doubly  stochastic  Poisson  point  process  (DSPP),^®^  which 
always  gives  rise  to  a  super-Polssonian  counting  distribution  (count  variance 
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greater  than  count  mean)'  has  been  treated  In  detail.  The  results  for 

(22)  (17) 

thermal  light  and  multlplled-Polsson  light  '  are  special  examples. 

Since  antibunched  llght^^’^^*^^’^^  appears  to  be  playing  an 

ever-increaslngly  Important  role  In  optics,  we  also  treated  several  cases  In 

which  the  count  variance  Is  less  than  the  count  mean.  Finally  we  studied  the 

combined  effects  of  random  deletion  and  additive  independent  Polsson-noise 

counts . 

Consider  a  random  number  n  of  events,  and  let  each  event  be  multiplied  by 
a  discrete  multiplication  or  reduction  factor  *  0,1,2,...  for  k  * 
l,2,...,n.  Then  the  total  number  of  multiplied  events  Is 


m  -  Z  x,  . 

k-1  k 


If  the  multiplication  factors  (x^)  are  statistically  Independent,  then  the 

random  variable  m  has  a  moment-generating  function  (mgf)  Q^s)  ■  <exp(-sm)>, 

(18) 

which  Is  related  to  the  mgf' a  of  n  and  x  by'  ' 


V8*  ■  Q„(_ln  V8)) 


The  means  and  variances  are  related  by 


<m>  -  <n><x> 


Var(m)  -  <x>  Var(n)  +  <n>Var(x), 
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When  the 


(25)  (26) 

Equation  (4)  is  known  as  the  Burgess  variance  theorem.  1 ’ v  1 

multiplication  factors  {x^}  are  Bernoulli  distributed,  i.e.,  x^  -  {1,0}  with 
probabilities  (n,  1  -  n). 


Qx(s)  -  1  -  n+  ne'8, 

<x>  ”  n ,  Var(x)  *  n(l  -  n).  (5) 

Substituting  Eqs.  (5)  into  Eqs.  (2)-(4),  we  obtain 

Q^s)  -  Qn(-ln(l  -  n  +  ne“8)), 

<m>  »  n<n>, 

Var(m)  ■  n^Var(n)  +  n(l  -  n)<n>.  (6) 

Equations  (6)  relate  the  properties  of  the  number  of  counts  after  Bernoulli 
deletion  to  those  before  Bernoulli  deletion,  in  terms  of  the  deletion 
parameter  n. 

By  repeated  use  of  Eqs.  (6),  it  can  be  easily  shown  that  successive 

random  deletions,  with  deletion  parameters  h  jhj,...,  are  ®<iuiv*l®nt  to  a 

single  process  of  random  deletion  with  parameter h  » 

Ve  consider  the  effect  of  random  deletion  on  a  representative  sub- 

(19) 

Polssonian  distribution.  Consider  the  counting  distribution 

11,  n-N 

P(n)  »  {  (7) 

(0,  n  4  N 

associated  with  the  number  state  (Fock  state)  |N>.  This  corresponds  to 

Qn(s)  -  exp(-sN), 
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<n>  “  N,  Var(n)  ■  0. 


(8) 


» 


t 


Using  Eqe.  (6),  we  obtain 

(^(s)  -  (1  -  ti  +ne‘8)N, 

<n>  •  nN, 

Var(a)  -  n(i  -  n)K  -  <n>(l  -  n),  (9) 


so  that  the  randomly  deleted  number-state  counting  distribution  is 
binomial. It  is  clear  from  Eqs.  (8)  and  (9)  that  both  the  number  state  and 
the  binomial  distributions  are  sub-Polssonian. 

Note  that,  if  a  binomial  distribution  of  parameters  (n,N)  is  further 
Bernoulli  selected  with  a  deletion  paramter  h* ,  the  resultant  distribution 

remains  binomial,  with  paramters  (nn  \  N). 

(19) 

We  have  also  shown'  that  neither  additive  independent  Poisson  noise 
events,  nor  Bernoulli  random  deletion,  alters  the  sub-  or  super-Poissonlan 
nature  of  the  light.  Indeed,  both  effects  are  quite  similar;  Increasing  the 
amount  of  either  drives  the  counting  distribution  toward  Poisson.  A  moment's 
thought  provides  the  reason:  The  deletion  serves  to  reduce  correlated  (or 
anticorrelated)  event  occurrences,  thereby  bringing  the  distributions  close  to 
the  zero-memory  Poisson  distribution. 

The  result  is  of  particular  importance  in  the  current  effort  to  produce 
sub-Polssonian  light.  It  shows  that  although  loss  (e.g.,  low  quantum 
efficiency)  and  additive  Poisson  noise  serve  to  reduce  the  observability  of 
the  sub-Polssonian  character  of  the  light,  they  do  not  destroy  it. 
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3.  Noise  In  the  detection  of  light 

Optical  detectors  often  operate  by  random  multiplication  to  increase  the 

responsivity .  Uell  known  examples  are  the  photomultiplier  tube  and  the 

avalanche  photodiode.  In  the  process  of  multiplying  the  signal,  the  detectors 

r  27) 

introduce  noise.  The  purpose  of  the  study  described  here  is  to 
investigate  this  noise. 

Multiplication,  reduction,  and  branching  processes  have  been  examined  in 
a  broad  variety  of  contexts.  '  '  '  Applications  range  from  astrophysics  to 

biological  information  transmission.  In  the  majority  of  mathematical 
treatments,  the  multiplication  or  branching  is  treated  as  an  Instantaneous 
effect.  However,  in  many  physical  systems,  a  random  time  delay  (or  spatial 
dispersion)  is  inherent  in  the  multiplication  process.  He  have  carried  out  an 
analysis  of  a  cascade  of  Poisson  multiplications  that  includes  such  time 
effects. Quit  results  reduce  to  previously  known  descriptions,  in  the 
limit  of  instantaneous  multiplication. 

He  developed  the  statistics  of  a  point  process  generated  by  a  cascade  of 
independent  Poisson  processes,  and  found  the  moment  generating  function,  as 
well  as  the  counting  and  time  statistics  when  dynamics  are  included.  Both  the 
stationary  and  nonstationary  cases  were  considered.  A  simple  expression 
for  the  variance-to-mean  ratio  at  the  mth  stage  has  been  obtained.  We  have 
carried  out  a  parametric  study  of  the  counting  distributions,  by  employing  the 
DEC  PDP  11/60  and  IBM  4341  computers. 

In  some  of  the  aforementioned  applications  of  cascaded  Poisson  processes, 
a  statistically  independent  additive  Poisson  point  process  may  also  be 
present,  representing  for  example,  broadband  background  light  and/or 
thermionic  emission  in  a  photomultiplier  tube.  The  counting  statistics  for 
the  superposition  process  can  be  simply  determined  by  the  rise  of  numerical 
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discrete  convolution.  We  are  in  the  process  of  considering  related  branching 
processes,  in  which  all  initiating  events  are  included  in  the  final  point 
process.  In  the  future,  this  work  will  be  applied  to  the  noise  properties  of 
real  detection  devices. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-ENG78-26498. 
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B.  SIGNAL-DEPENDENT  NOISE* 

(P.  R.  Prucnal,  M.  C.  Teich) 

(JSEP  work  units  1  and  A,  1982-1985) 

(Principal  Investigators:  M.  C.  Teich  (212)  280-3117, 

E.  S.  Yang  (212)  280-3120, 

W.  Hwang  (212)  280-3115) 

Signal-dependent  noise  is  manifested  in  many  physical  processes  and 
often  represents  a  fundamental  limitation  to  the  minimum  uncertainty  in  the 
magnitude  of  a  physical  quantity.  When  this  quantity  is  the  number  of  photons 
incident  on  a  photodetector,  the  uncertainty  in  the  number  of  photons  is 
proportional  to  the  total  energy  received.  This  signal-dependent  photon  noise 
places  a  limit  on  the  minimum  light  energy  that  can  be  received  with  a  given 
reliability.  When  the  quantity  of  Interest  is  the  number  of  impurities 
Implanted  in  a  semiconductor  material,  the  uncertainty  in  the  number  of 
impurities  is  proportional  to  their  average  number.  This  signal-dependent 
"impurity-noise"  places  a  limit,  in  certain  circumstances,  on  the  minimum 
device  dimensions  in  a  VLSI  system  that  will  function  with  a  given  reliability. 

Substantial  progress  has  been  made  in  understanding  the  role  of  signal- 
dependent  noise  in  each  of  these  situations:  the  performance  of  source  modula¬ 
tion  formats  in  fiber  optic  communications;  the  performance  of  VLSI  systems 
containing  ion-implanted  bipolar  transistors.  These  results  are  summarized 
below  and  reported  in  detail  in  the  literature. 


1.  Effect  of  signal-dependent  noise  on  VLSI  yield 

An  expression  was  derived  for  the  doubly-stochastic  distribution  of  the 
number  of  impurities  in  the  base  region  of  a  bipolar  transistor;  the 
distribution  results  from  signal-dependent  uncertainty  in  ion  implantation 
parameters.  Expressions  were  derived  for  device  yield,  and  VLSI  (very  large 
scale  Integration)  chip  yield  with  an  N-bit  parity  check.  These  derivations 
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can  be  extended  to  other  devices  in  a  straightforward  manner. 

The  following  notation  was  used: 

L  emitter  dimension,  defined  as  the  fundamental  dimension 

in  the  technology 

C^(x)  emitter  donor  impurity  concentration 

Al*  rl»  8i  amplitude,  range,  and  straggle  of  Cdl(x) 

Ca2(x)  base  acceptor  impurity  concentration 

A2*  r2’  S2  amplitude,  range,  and  straggle  of  Cg2(x) 

amplitude  of  collector  donor  impurity  concentration 
Gfa,  Gg  base,  emitter  Gummel  number 

k  scaling  factor 

n  number  of  base  impurities 

Pn  average  number  of  base  impurities 

0  current  gain 

Dnb*  Dpe  diffusion  coefficient;  electrons  in  base,  holes  in  emitter 

A  tolerable  fractional  range  of  variation  of  0  about  G{0} 

Nchip  number  of  devices  per  chip 

R  number  of  s-dependent  devices  in  a  given  region  of  chip 

N  number  of  bits  per  word 

Yj  D  device  yield 

*chip  D  chip  yield,  no  parity  check 

Y_  „  yield  of  N-bit  word  with  parity  check 

Ychip,D;N  chlp  yleld»  N-blt  Parlty  check 
blnm  (),  poim  ()  pmf  of  binomial  and  Poisson  distributions 
The  following  assumptions  were  made: 

1.  The  device  structure  is  typical  of  recent  1  urn  bipolar  technology. 
See  Fig.  1,  for  details. 
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2.  The  collector  region  has  uniform  (substrate)  doping  concentration.  In 
a  real  device  this  doping  is  not  uniform  but  Increases  proceeding  deeper  Into 
the  substrate  (away  from  the  base)  due  to  the  presence  of  a  burled  N+  layer. 
The  collector  doping,  Is,  however,  constant  over  the  region  of  Interest  here. 

3.  The  base  region  Is  formed  by  implanting  acceptor  (boron)  Impurities. 
The  emitter  region  Is  formed  by  Implanting  donor  (phosphorous  or  arsenic) 
impurities. 

4.  The  Impurity  concentrations  are  Gaussian  in  shape.  The  effects  of 
channeling  on  the  distribution  (which  modify  the  Gaussian  shape)  are  minimized 
by  implanting  off-axis  of  the  crystal. 

5.  The  doping  is  uniform  in  the  plane  parallel  to  the  silicon  surface. 

6.  For  the  narrow  base  regions  considered,  the  base  transport  factor  is 
unity. 

7.  6  »  1. 

8.  The  n  is  a  (signal-dependent)  Poisson  r.v.,  provided  the  impurities 
are  s-independent,  that  is,  the  presence  of  an  impurity  does  not  influence  the 
presence  of  an  adjacent  impurity. 

9.  Any  of  the  parameters  r^,  r2,  s^,  s2  may  be  stochastic.  The  Gaussian 
distribution  is  used  as  an  example,  although  another  distribution,  such  as  the 
uniform,  could  be  used  as  easily.  For  a  1  ia  structure  E{r^}  -  0,  Efr^)  ■ 

0.05  urn  EtSj}  -  E{s2}  -  0.02  ym,  A^  -  1020cm"3,  A2  -  1018cm’3,  A3  -  1016cnT3. 

Calculations  have  been  performed  using  specific  parameters,  which  have 
led  to  the  following  observations:  (see  Figs.  1  to  9  for  details) 

1.  The  doubly  stochastic  effect  is  most  sensitive  to  uncertainty  in  the 
straggle  (standard  deviation)  of  the  emitter  impurity  distribution. 

2.  Uncertainty  of  the  order  of  5Z  in  an  implantation  parameter  causes 
substantial  broadening  of  the  distribution  of  impurities,  for  the  case 
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Fig. 


COLLECTOR  BASE  EMITTER 


Fig.  1  Ion- implanted  bipolar  transistor  structure. 


2  Dashed  curve  represents  singly-stochastic  distribution  of 
base  impurities  for  emitter  dimension  L  *>  1.  Solid  curve 
represents  doubly-stochastic  distribution  of  base  impurities 
corresponding  to  Gaussian  fluctuations  in  device  parameter. 
Cases  r2,  S ^ ,  and  S2  are  shown.  Three  standard  deviations 
of  device  parameter  equals  1%  of  mean.  Here  E{rj}  *  0  pm, 
E{r2)  -  0.05  um,  EtS^  -  E{S2}  -  0.02  pm,  hi  -  108/(pm)3, 

A2  -  106/(ym)3,  A3  -  104/(um)3. 
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Fig.  3  Same  as  Fig.  2,  except  for  fluctuations  in  Sj  only. 

Three  standard  deviations  of  equals  1/2%,  1%,  2%, 
5%  of  mean. 


tMITTtH  DIMENSION  l  <Mm) 


Fig.  A  Device  yield  versus  emitter  dimension,  where  tolerable  range 
of  variation  of  g  is  E{g}  +  2.4%.  Here  E{ri}  **  0  ym,  E{r2)  * 

C.U5  L  ym.  E{Si)  »  E{S2 J  *  0.2L  ym,  Aj  ■  10®/ [L^(ym) 3] ,  A2  * 
10®/(L2(ym)3] ,  A3  «  10^/[L2(ym)3] .  Solid,  dotted,  and  dashed 
curves  correspond  to  fluctuations  in  Sj,  r2  and  S2  respectively. 
Three  standard  deviations  of  device  parameter  equals  1/2%  of  mean. 


Lo«  nCmi, 


Fig.  8  Minimum  device  yield  necessary  for  chip  yield  greater  than 
0.6  versus  log  Nchip.  for  L  -  1.  Cases  of  no  parity  check 
and  N  «  2,4,8,16. 


EMITTER  DIMENSION  Ll^nO 

Fig.  9  Chip  yield  versus  emitter  dimension  for  16-bit  parity  check. 

Here  E{r,}  -  0,  E{r2}  -  0.5  L  um,  E{S, }  -  E{S2}  =  0.02L  urn, 

A^  ■  10®/ [L2(pm) 3] ,  a2  *  i0°/[L2(ym)^J ,  and  A3  *  10^/[L^(pm)3] . 
Cases  Sj,  r2,  and  S2  shown.  Three  standard  deviations  of  device 
parameter  equals  1 / 2%  of  mean.  Tolerable  range  of  variation  of 
8  is  E{8}  +  2.4%.  Number  of  devices  equals  100  for  solid  curves 
and  100/L2~for  dashed  curves. 


considered. 


3.  Device  yield  decreases  rapidly  for  dimensions  less  chan  a  well-defined 
threshold  (~  0.75  pm  for  the  case  considered). 

4.  Chip  yield,  without  a  parity  check,  exhibits  a  threshold  effect  at 
device  yield  *  ^“^^chip*  ^chip  *  nun,ber  of  devices  per  chip.)  The  device 
yield  must  exceed  this  threshold  to  produce  large  chip  yields. 

5.  The  use  of  a  parity  check  reduces  the  device  yield  threshold  to 
1-10/N^^.  Use  of  fewer  bits  per  parity  check  reduces  the  threshold  further. 

6.  For  the  example  considered,  the  minimum  device  dimensions  for  large 
chip  yields  is  of  the  order  of  1  to  1.5  um,  using  a  16-bit  parity  check. 

The  minimum  device  size  for  reliable  system  performance  for  other  cases 
will  depend  upon  specific  device  parameters.  However,  it  is  apparent  that 
chip  yield  is  extremely  sensitive  to  individual  device  yield,  and  that  it 
decreases  with  increasing  uncertainty  in  device  parameters,  smaller  device 
dimensions,  greater  number  of  devices  per  chip,  larger  number  of  bits  per 
parity  check,  and  with  smaller  allowed  ranges  of  tolerable  device  behavior. 

2.  Signal-dependent  noise  in  frequency-shift  keyed  fiber  communication 

An  optical  fiber  communication  system  was  implemented  using  frequency- 
shift  keying  (FSK)  and  amplitude-shift  keying  (ASK).  The  effect  of  photon 
noise  on  performance  was  studied  theoretically  and  experimentally.  The 
bit-error  rates  of  the  FSK  and  ASK  systems  were  calculated  as  follows.  The 
mean  voltage  at  the  avalanche  photodiode  receiver  output  is 

<Vij>  ”  K1<G>Er/hf  (1> 

where  v^  is  the  voltage  in  channel  i  given  that  bit  j  was  transmitted,  K^  is 
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» 

an  Irrelevant  constant,  hf  is  the  photon  energy,  and  E^  is  the  received  pulse 
t  energy,  taking  into  account  the  extinction  ratios  of  the  lasers,  the  cross¬ 

talk,  and  all  optical  losses.  The  variance  of  the  voltage  is 

*  var  v1J  -  K1<GXviJ>F  +  K2Z2  (2) 

where 

F  -  k<G>  +  (l-k)(2-l/<G»  (3) 

is  the  APD  excess  noise  factor,  and 

1/2 

Z  -  (NEP)BA/  <G>T/hf  (4) 

is  the  preamplifier  noise  parameter  at  room  temperature,  with  T  the  bit 
interval  and  B  the  amplifier  bandwidth. 

Using  a  Gaussian  approximation  for  the  probability  density  function  of 
the  voltage,  tho  FSK  bit-error  rate  is 

*  BERFSK  -  (l/4)erfc[2~1/2«v11>-<v£)1»/(var  vn+var  vQ1)1/2] 

+  (l/4)erfc{2  1/2(<v0Q>-<v10>)/( var  vQ()+var  v1Q)1/2].  (5) 

*  The  ASK  bit-error  rate  is 

BERASK  -  (l/4)erfc{2”1/2(<v11>-6)/(var  vu)1/2] 
f  +  (l/4)erfc[2_1/2(9-<v10»/(var  v1Q)1/2]  (6) 

where  9  is  the  threshold. 

» 

The  bit-error  rate  for  FSK  and  ASK  corresponding  to  theory  (curves)  and 
data  (dots)  is  shown  in  Fig.  10.  Both  the  theory  and  the  data  in  Fig.  10  show 
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approximately  a  3  dB  improvement  of  FSK  over  ASK.  The  advantage  of  FSK 
increases  slightly  as  the  optical  power  increases.  This  effect  can  be 
explained  by  the  performance  of  FSK  relative  to  ASK,  as  the  contribution  of 
signal-dependent  noise  increases. 

Indeed,  the  theoretical  curves  in  Fig.  10  reflect  the  contributions  of 
both  shot  noise  and  thermal  noise.  For  the  parameters  of  this  system,  the 
shot  and  thermal  noise  powers  in  Eq.  (2)  are  equal  when  the  received  optical 
power  is  1.2  nW.  Thus,  the  shot  noise  is  comparable  to  the  thermal  noise  at 
the  right  boundary  of  Fig.  10  and  dominates  the  thermal  noise  at  the  left 
boundary. 


OPTICAL  LOSS  ;dB) 

52  5  54  55  56  $?  58 


Figure  10 
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That  the  advantage  of  FSK  Increases  as  the  shot  noise  increases,  is  seen 
as  follows.  Since  the  FSK  receiver  subtracts  the  output  of  two  detectors, 
both  the  signal  energy  and  thermal  noise  power  are  doubled,  compared  to  ASK. 

The  advantage  of  FSK  in  the  thermal  noise  dominated  regime  is  then  1.5  dB.  On 
the  other  hand,  the  shot  noise  in  the  FSK  receiver  is  the  same  as  for  a 
transmitted  "1"  in  the  ASK  receiver.  The  advantage  of  FSK  in  the  shot  noise 
dominated  regime  is  then  6  dB,  for  a  transmitted  ”1”.  For  a  transmitted  ”0", 
ASK  has  no  shot  noise,  and  FSK  may  or  may  not  have  an  advantage,  depending  on 
the  relative  magnitudes  of  the  shot  and  thermal  noise  powers.  It  is  pointed 
out  that  the  SNR-type  arguments  above  are  not  necessarily  sufficient  to 
describe  system  performance,  and  are  used  here  only  to  give  an  intuitive 
understanding  of  the  results. 

It  is  apparent  from  the  above  considerations  that  the  advantage  of  FSK 
may  be  greater  in  the  shot-noise  dominated  regime.  In  fact,  other  types  of 
noise,  with  stronger  signal  dependence  than  shot  noise,  could  theoretically 
result  in  a  greater  advantage  for  FSK.  This  is  illustrated  by  considering  a 
general  class  of  signal-dependent  noise  in  which  the  variance  is  related  to 
the  mean  by  the  power  law 

var  v  -  c<v>2p  (7) 

where  c  is  a  nonnegative  constant  and  0<p<l.  Equating  the  SNRs  for  a 
transmitted  "1"  for  FSK  and  ASK  and  simplifying  yielcs 

2<VFSK>  P  "  <VASK>  P' 

The  decrease  in  minimum  optical  power  required  by  FSK  for  the  same  performance 


The  decrease  in  minimum  optical  power  required  by  FSK  for  the  same  performance 
as  ASK  is 


D  -  3/(l-p)  dB. 


(9] 


When  p  ■  1/2,  the  variance  is  proportional  to  the  mean,  as  for  shot  noise, 
and,  D  *  6  dB,  in  agreement  with  the  above  discussion.  As  p  Increases  in  the 
region  0<p<l,  D  also  increases.  Therefore,  greater  signal  dependence 
increases  the  advantage  of  FSK.  This,  and  other  advantages  of  FSK,  such  as 
Immunity  to  threshold  variations,  are  discussed  in  reference  (2). 


*Thi8  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-ECS80-06881 . 

(1)  P.  R.  Prucnal,  "Wideband  FSK  fiber-optic  communication  experiment,"  J. 
Opt.  Comm.,  issue  // 2 ,  in  press.  (1983). 
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II.  PHYSICAL  AND  PHOTOCHEMICAL  PROPERTIES  OF  ELECTRONICS  MATERIALS 

A.  CARRIER  RECOMBINATION  AT  GRAIN  BOUNDARIES  AND  THE  EFFECTIVE  RECOMBINATION 
VELOCITY* 

(U.  Hwang,  E.  Poon,  H.  C.  Card+) 

(JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  V.  Hwang  (212)  280-3115 

E.  Yang  (212)  280-3120) 

In  problems  involving  the  recombination  of  excess  minority  carriers  at 

semiconductor  interfaces  such  as  grain  boundaries  and  oxidized  surfaces,  there 

often  arises  a  confusion  between  the  actual  minority  carrier  recombination 

velocity  at  the  interface  Itself  and  the  effective  collection  velocity  for 

minority  carriers  at  the  edge  of  the  adjacent  space  charge  region.  In 

modelling  device  behavior  it  is  usually  convenient  to  consider  the  latter 

mechanism,  since  under  these  conditions,  the  carrier  removal  can  be  considered 

in  terms  of  an  'effective'  recombination  velocity  at  the  boundary  of  a 

field-free  region.  One  specific  example  involves  the  disappearance  of  the 

photogenerated  minority  carriers  at  the  grain  boundary  in  photovoltaic  solar 

cells.  A  phenomenon  of  some  practical  importance  in  the  design  and  simulation 

of  solar  cell  devices  is  the  appropriate  description  of  the  dependence  of  the 

effective  recombination  velocity  upon  the  photogeneration  rate,  or  optical 

illumination  Intensity.  It  is  customary  in  such  cases  to  associate  a  constant 

recombination  velocity  with  the  boundary  of  the  field-free  region.  Recent 

(1  2) 

experimental  studies  ’  indicate,  however,  that  this  velocity  is  not  a 

constant,  but  rather  shows  a  substantial  dependence  upon  the  photogeneration 

rate,  or  in  other  words,  on  the  recombination  current  density. 

A  specific  example  will  clarify  the  distinction  between  the  true 

interface  recombination  velocity  and  the  effective  recombination  velocity, 

(3) 

which  in  many  cases  is  a  more  convenient  variable. 
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Ac  Che  grain  boundary  (GB)  interface,  Che  recombinaCion  currenC  densiCy 

(4) 

is  as  described  in  a  previous  paper.  The  inCerface  recombinaCion  velociCy, 
S,  aC  Che  GB  is  defined  in  Cerms  of  Che  recombinaCion  currenC  density,  JR,  and 
Che  excess  minoricy  carrier  concenCradon  aC  Che  GB  (x*0)  for  n-Cype  silicon  as 

S(0)  -  JR/qtp(0)-p0]  (1) 

where  p^  is  Che  equilibrium  concenCradon  of  Che  minoricy  carriers  aC  Che  GB. 
For  device  calculaCions ,  Che  "effecdve"  recombinaCion  velociCy  aC  the  edge  of 
the  depletion  region,  S(w),  is  given  by  ^ 

S(w)  -  JR/2q[p(w)  -  pnQ]  (2) 

where  pnQ  refers  to  the  equilibrium  concentration  of  holes  in  the  neutral 
region  and  p(w)  is  given  by 


p(w) 


p(0)e 


-qVd/kT 


+  bx 


/qV./kT 

;  d 

o 


,-y 


dy 


(3) 


under  assumption  JR(0)  *  2JR(w).  (No  recombinaCion  or  generation  within  the 


space  charge  region).  is  Che  GB  diffusion  potential. 


2D  q 
P 


and  x 


/  E  D 

^  °pTr  /q  Ndup 


At  the  grain  boundary  plane  are  localized  interface  sCates,  which  in 
general  are  statistically  distributed  in  energy  over  the  energy  gap.  These 
states  have  a  distribution  determined  by  the  nature  of  the  disorder,  by 
dangling  bonds,  and  by  the  local  electronic  potential  fluctuations  producing 
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stress  fields  In  the  regions  of  structural  irregularities.  They  aay  also 

arise  from  chemical  impurities  trapped  between  adjacent  grains.  Together 

/ 

dangling  bonds  and  local  structural  strain  can  produce  an  inhomogeneous 
broadening  of  the  discrete  levels  in  the  gap  if  their  effects  are  sufficiently 
severe. 

As  a  feasible  approximation  to  the  energy  spectra  for  these  localised 
states,  Gaussian  energy  distributions  are  adopted.  These  states  are 
classified  into  two  types,  "acceptor-like”  and  "donor-like"  Interface  states 
htA(E)  and  and  are  expressed  by 


htA<E> 


(3) 


and 


(4) 


where  HtA  and  are  the  total  density  of  the  localized  states  per  unit  area, 
8^  and  8q  are  distribution  parameters  and  E£A  and  E£p  are  the  mean  value  of 
the  acceptor-like  and  donor-like  distributions,  respectively.  The  exact 
energy  and  type  of  traps  which  are  present  in  polysilicon  are  not  completely 
understood.  It  is  generally  agreed  that  the  trapping  levels  are  near  midgap 
on  the  basis  of  the  observed  activation  energies  and  bias  dependence  from 
transport  measurements.^  An  Important  feature  of  the  "deep-level"  Gaussian 
distribution  model  suggested  here  is  the  energy  position  of  the  mean  value  of 
the  "acceptor-like"  and  "donor-like"  distributions.  We  take  these  to  be, 
respectively,  slightly  above  and  slightly  below  midgap  as  in  the  case  of 
amorphous  Si.  '  The  density  and  distribution  (broadening)  parameters  of  the 


] 
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GB  trapping  states  are  expected  to  depend  on  the  conditions  of  sample 
preparation,  orientational  mismatch,  passivation  and  annealing  processes. 
Typical  results  are  presented  for  the  following  cases  of  the  GB  interface 
state  distribution: 

Case(a)  :  E£A  -  0.5  eV,  sA  ■  10  kT 
Case(b)  :  E  -  0.56  eV,  sA  -  10  kT 
Case(c)  :  EfcA  “0.5  eV,  sA  ■  kT 
Case(d)  :  E  -  0.56  eV,  sA  ■  kT. 


In  all  cases  considered  here,  the  donor-like  states  are  assumed  to  be 

filled  with  electrons  and  therefore  neutral.  Other  values  used  in  the 

11  -2 

calculations  are:  H£A  *  10  cm  ,  doping  concentration  within  the  grain  Nd 

■  2x10*^  cm  ^  e  »  11.8  e  ,  V..  •  10^  cm/sec,  a  **  10  cm^  and  o„  »  10  *** 

8  u  tn  c  N 

2 

cm  . 

The  latter  two  parameters  correspond  to  capture  cross-sections  associated 


with  coulomblc  and  neutral  capture  respectively,  of  carriers  by  the  Interface 


states 


(4), (6) 


Figures  1  and  2  show  the  dependence  of  the  Interface  recombination 


velocity  S(0)  and  the  effective  recombination  velocity  S(w)  on  the  value  of 

(4) 

AEj,  respectively.  These  results,  as  reported  earlier,  show  that  S(0) 
exhibits  an  Increase  with  increasing  AE^  for  small  values  of  AE^,  then  reaches 
saturation  at  high  AE^. 

The  velocity  S(w)  should  not  be  taken  to  be  constant,  but  rather 


Increases  with  increasing  AE^,  reaches  a  peak  at  a  critial  value  of  AE^,  and 
then  decreases  with  increasing  AE^  (with  high  illumination  levels)  in  all  four 

cases. 
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In  Figs.  3  and  4  we  show  the  interface  recombination  velocity  and  the 
effective  velocity  as  a  function  of  the  recombination  current  J^.  Fig.  5 
shows  the  dependence  of  the  effective  velocity  S(w)  on  p(w),  the  excess 
minority  carrier  concentration  at  the  edge  of  the  neutral  region.  Figs.  4  and 
5  are  of  most  utility  to  device  modelers  because  they  simply  express  a 
nonlinear  boundary  condition  for  effective  recombination  rate  at  the  edge  of  a 
field-free  region.  The  solution  to  the  continuity  equation  in  the 
quasi-neutral  region  is  obtained  with  a  nonlinear  boundary  condition,  but  one 
which  does  not  contain  additional  unknowns,  such  as  AE^  at  the  grain  boundary. 

These  results  are  qualitatively  consistent  with  those  of  Ref.  (9),  which 
considered  the  effective  recombination  velocity  for  grain  boundaries  with 
monoenergetic  (midgap)  recombination  centers. 

It  has  been  demonstrated  in  the  above  example  that  in  grain  boundary 
recombination  processes,  it  is  possible  for  the  interface  recombination 
velocities  to  increase  over  a  range  of  recombination  rates  (illumination 
intensities),  while  at  the  same  time  the  effective  recombination  velocity  of 
the  minority  carriers  decreases  with  the  recombination  rate.  This  has  arisen 
because  the  diffusion  potential  created  by  the  interface  charge  decreases  with 
recombination  rate  in  a  way  which  overcompensates  for  the  increasing 
recombination  velocity. 

It  is  therefore  essential  that  in  modelling  interface  recombination  in 
devices,  whether  this  be  at  surfaces  or  at  grain  boundaries,  one  distinguishes 
between  the  true  recombination  velocity  which  depends  only  upon  the 
recombination  process,  and  the  effective  velocity  at  the  boundary  of  the 
field-free  region,  which  depends  also  upon  the  minority  carrier  transport 
through  the  adjacent  space-charge  region. 
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Recombination  Current  Density,  JR  (amp/cm 


Figure  5.  Effective  recombination  velocity  at  the  edge  of  the  neutral  region  as 

a  function  of  excess  minority  carrier  concentrations  at  this  location  (cases  a-d) 
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B.  ELECTRICAL  CHARACTERISTICS  OF  LARGE  SCALE  INTEGRATION  (LSI)  MOSFETS 

AT  VERY  HIGH  TEMPERATURES 

(F.  Shoucair,  W.  Hwang) 

( JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  W.  Hwang  (212)  280-3115 

E.  Yang  (212)  280-3120) 

The  main  objective  of  this  work  has  been  the  study  of  the  effects  of  high 
temperature  on  the  electrical  behavior  of  silicon  MOS  transistors  and 
circuits. ^)*(^)»(3)  formuiation  of  the  theory  was  based  on  the  physical 

structure  of  a  commercially  available  CMOS  process  with  junction  Isolated  N 
and  P  channel  MOSFETs,  but  the  results  obtained  are  valid  for  dielectrically 
Isolated  transistors  (DIMOSFETs)  and  silicon  on  Insulator  (SOI)  devices  as 
well,  especially  with  regard  to  the  modelling  of  leakage  currents. 

The  theoretical  approach  consisted  in  examining  the  MOSFETs'  temperature 
dependent  electrical  parameters  individually.  They  were  found  to  Include  the 
threshold  voltage,  the  channel  mobility,  and  the  channel  length  modulation, 
the  latter  parameter  being  relevant  only  to  devices  operated  in  the  saturation 
region.  Necessary  and  sufficient  conditions  for  the  existence  of  Zero- 
Temperature-Coefficient  (ZTC)  gate  bias  voltages  have  been  developed.  The 
procedure  followed  for  this  purpose  relied  on  the  mutual  cancellation  of  the 
threshold  voltage  and  the  channel  mobility  temperature  derivatives  over  a 
specified  temperature  range.  It  has  been  found  that  a  given  transistor  has 
one  ZTC  gate  bias  voltage  in  each  one  of  the  nonsaturation  and  the  saturation 
regions  of  operation,  and  that  the  two  gate  voltages  V„c  (ZTC)  only  depend  on 

vd 

the  parameters  p^  and  of  the  threshold  voltage  temperature  dependent 
function  of  the  device. The  small  signal  parameters  of  the  MOSFET  biased 
at  its  ZTC  points  were  obtained  analytically,  and  it  was  shown  that  small 
signal  capacitances  are  very  weakly  dependent  on  temperature. 
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Experimental  data  collected  on  long  N  and  P  channel  MOSFETa  in  the  range 

t  25°C  to  300°C  has  revealed  that  (i)  the  threshold  voltage  magnitude  |v_u(T)[ 

decreases  linearly  with  increasing  temperature,  (11)  the  channel  mobility  p(T) 

displays  a  T  dependence,  and  (ill)  the  drain  to  substrate  leakage  currents 

exhibit  proportionality  to  n^(T),  the  intrinsic  carrier  concentration,  and  to 
2 

n^(T),  in  the  lower  and  upper  parts  of  the  temperature  range  under 
consideration,  respectively. 

Current  voltage  characteristics  have  been  obtained  In  the  subthreshold, 

nonsaturation,  and  saturation  regions  of  operation.  Subthreshold  currents 

have  been  found  to  increase  exponentially  with  temperature,  thereby  severely 

limiting  the  range  over  which  MOSFETs  can  be  used  in  this  bias  regime.  The 

nonsaturation  and  saturation  region  current  characteristics  were  shown  to 

present  ZTC  bias  points  to  gate  voltages  V  (ZTC)  in  close  agreement  with  the 

G5 

( 7 ) 

theoretical  predictions. 

The  experimental  data  has  been  used  to  extend  the  validity  of  the 

conventional  (room  temperature),  large  and  small  signal  models  of  the  MOSFET 
I  o 

'  up  to  300  C.  It  has  been  shown  that  the  effects  of  high  temperature  on  the 

terminal  characteristics  of  MOSFETs  could  be  adequately  modelled  by  the 
Inclusion  of  a  reverse  biased  diode  between  the  drain  and  substrate  terminals 

*  of  these  devices.  This  diode  conducts  the  leakage  currents  of  the 
drain-substrate  Junction,  which  increase  exponentially  with  temperature.  The 
high  temperature  small  signal  model  was  found  to  be  identical  to  the  room 

* 

temperature  one,  provided  the  correct  temperature  dependences  of  the  various 

(3) 

electrical  components  are  used. 

The  effects  of  increasing  temperature  on  the  dc  characteristics  of  a  CMOS 
Inverter  were  shown  to  be  (1)  decreasing  noise  margins  and  transition  region 
slope  between  2S°C  and  270°C,  and  (11)  the  triggering  of  a  pnpn  latchup 
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phenomenon  near  270°C.  At  270°C,  the  switching  speed  of  such  Inverters  was 
found  to  incur  a  67Z  reduction,  on  average,  with  respect  to  that  at  25°C. 

Resistive  load  inverters  biased  near  (N  channel  MOSFETs)  or  at  (P  channel 
MOSFETs)  their  saturation  region  ZTC  points  have  been  shown  to  present  a 
decrease  in  small  signal  gain  with  increasing  temperature,  in  close  agreement 
with  the  predictions  of  the  models  developed.  Some  implications  of  these 
models,  on  the  design  and  analysis  of  analog  MOS  circuits  have  been  discussed, 
whereby  it  has  been  shown  that  the  main  effect  of  high  temperature  is  the 
introduction  of  a  constraint  on  the  choice  of  MOSFETs'  gate  widths.  This 
constraint  arises  because  of  leakage  current  considerations,  and  is  necessary 
in  order  to  avoid  the  drift  of  bias  quantities. 

The  results  presented  in  this  work  have  demonstrated  that  silicon  MOSFETs 
can  be  used  in  very  high  temperature  environments  for  digital  and  analog 
applications.  In  view  of  the  fact  that  the  CMOS  process  investigated 
experimentally  was  designed  for  room  temperature  operation,  this  process  shows 
great  promise  as  a  high  temperature  silicon  technology. 
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Experiment,"  (submitted  for  publication). 

(3)  F.  Shoucair,  W.  Hwang,  and  P.  Jain,  "Electrical  Characteristics  of  Large 
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Modeling  and  Circuit  Behavior,”  (submitted  for  publication). 
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c.  PHONON-ASSISTED  CARRIER  TRANSPORT  ACROSS  A  GRAIN  BOUNDARY* 

(C.  M.  Wu,+  E.  S.  Yang) 

(JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  E.  S.  Yang  (212)  280-3120 

W.  Hwang  (212)  280-3115) 

The  electrical  properties  of  polycrystalline  silicon  are  extensively 
studied  because  of  its  importance  in  the  development  of  Integrated  circuits. 
The  polysilicon  prepared  by  chemical  vapor  deposition  (CVD)  usually  has  a 
grain  size  of  1  pm  or  less.  The  control  of  conductivity  of  the  CVD  poly¬ 
silicon  is  important  when  it  is  utilized  as  an  interconnection  or  a  distri¬ 
buted  resistor.  Seto^^  proposed  a  model  Introducing  the  trapping  states  to 
describe  the  GB  electrical  properties.  Following  the  same  approach,  other 
researchers^^  ^  have  studied  samples  of  both  large  grain-size  wafers  and 
small  grain-size  polyslllcon  films.  Lu  et  al.  used  an  attenuation  factor  in 
the  form  of  a  scattering  potential  to  describe  the  phonon  scattering  at  the 
GB.  The  reported  data  show  that  the  majority-carrier  conduction  is  signifi¬ 
cantly  less  than  that  predicted  by  the  thermionic  emission.  This  conclusion 
is  reasonable  since  the  transport  of  carriers  from  one  grain  to  the  other 
require  not  only  enough  energy  to  overcome  the  potential  barrier,  but  also  an 
appropriate  momentum  in  the  second  grain.  In  other  words,  the  quantum 
mechanical  effect  does  not  allow  the  incident  carrier  to  exist  in  the  second 
grain  if  no  suitable  state  is  available.  This  selection  rule  causes  the 
reduction  of  carrier  transmission.  Since  the  current  conduction  due  to  the 
emission  of  carriers  trapped  by  the  GB  states  is  limited  (5.5  percent  of  total 
carriers  leaving  the  GB),  we  will  only  discuss  the  carrier  transmission 
directly  from  one  grain  to  the  other.  In  addition  to  phonon  scattering, 
Coulomblc  Interaction  with  interface  charge  can  be  significant,  and  it  has 
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been  reported  separately 


(6) 


Thermionic  emission  has  been  discussed  In  detail  for  carrier  transport 
across  a  metal-semiconductor  junction. ^  In  some  cases,  the  effects  of 
tunneling  and  quantum  mechanical  reflection  may  also  Influence  the  current 
conduction.  Since  the  metal  is  always  degenerated  and  the  electron  density  is 
very  high,  a  carrier  can  be  emitted  freely  from  the  semiconductor  without  any 
limitation  on  its  momentum  as  long  as  the  kinetic  energy  is  large  enough  to 
overcome  t>e  potential  barrier.  But  this  is  not  true  for  a  carrier  travelling 
from  one  semiconductor  to  another,  1  where  the  probability  of  the 

scattered  carrier  to  occupy  an  appropriate  state  in  the  second  semiconductor 
or  grain  must  be  taken  into  account.  We  consider  the  case  of  electron 
transport  from  one  conduction-band  valley  In  the  first  grain  (valley  a)  to  one 
of  the  valleys  in  the  second  grain  (valley  b).  This  transition  from  the  wave 
vector  k^  to  the  other  wave  vector  kfe  usually  occurs  with  either  absorption  or 
emission  of  a  phonon.  Thus  the  thermionic  emission  current  from  the  first 
grain  (z  <  0)  to  the  second  grain  (z  >  0)  is  given  by 

J1~2  ■  *1  T<iv*)  <*> 

x  z 

where  the  summation  is  over  the  ellipsoids. 

Since  T  is  normally  less  than  unity,  (1)  explains  why  an  attenuation 
factor  is  needed  to  fit  the  experimental  data  in  the  previous  work.^^  Let 

us  compare  (1)  to  the  classical  form  of  the  thermionic  emission  current  given 

by 


jjjdSjiirt,) 


(2) 
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It  is  now  evident  that  the  average  effect  of  the  transmission  coefficient  T 
gives  rise  to  the  attenuation  factor  which  is  defined  as 


» 


n 


(3) 


It  is  possible  to  calculate  the  current  in  (1)  numerically  if  T(k^,  K)  is 
given  as  a  function  of  the  kinetic  energy  and  the  assisting  phonon  energy 
4u»l. 

If  the  first  grain  is  located  on  the  Cartesian  coordinate  system,  and  the 
orientations  of  (100),  (010),  and  (001)  are  along  the  x^-,  and  x^~  axes, 

respectively,  three  parameters  are  needed  to  specify  the  orientation  of  the 
second  grain. 

For  a  cubic  crystalline  lattice  such  as  silicon,  the  k-axes  have  the  same 
directions  as  the  real-space  axes.  Although  the  origins  of  the  real  space  in 
the  two  grains  may  be  displaced,  the  major  axes  of  the  two  grains  in  the 
k-space  have  the  same  origin.  Thus  the  carrier  at  the  Interface  should  have 
zero  momentum  in  the  second  k-space  if  it  has  a  zero  momentum  in  the  first 
k-space.  Based  on  this  argument,  the  k-space  of  the  second  grain  is  a  system 
obtained  by  the  three  rotations  of  the  k-space  of  the  first  grain.  Fig.  1(a) 
shows  the  three  rotations.  The  system  is  first  rotated  along  the  x^-axis, 
then  the  x^-axis.  Finally,  the  system  is  rotated  along  the  x^-axis.  The 
rotation  angles  are  9,  4>,  and  respectively. 

As  discussed  previously,  the  current  is  caused  by  the  transition  of 
carriers  from  one  of  the  energy  valleys  in  the  first  k-space  to  the  energy 
valley  in  the  second  k-space  (Fig.  1(b)).  The  wave  vector  of  the  phonon 
ranges  from  0  to  (4/3)  (ir/a).  It  could  be  either  an  acoustic  phonon  or  an 
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optical  phonon.  The  phonon  frequency  is  in  the  order  of  10  Ha  and  the 
phonon  energy  can  be  as  high  as  0.06  eV  in  silicon,  which  is  larger  than  its 
thermal  energy  at  room  temperature.  Therefore,  at  large  angles,  the  number  of 
phonons  available  is  so  limited  that  the  carrier  transition  is  constrained. 

By  examining  the  rotated  axis  in  the  octant  of  the  original  vector  space, 
(Fig.  1(b)),  we  can  always  find  an  original  axis  which  has  an  angle  less  than 
54.7°.  In  other  words,  the  carrier  in  the  valley  of  the  second  grain  can 
always  make  a  transition  with  the  assistance  of  less  than  0.919  kQ  that  cor¬ 
responds  to  the  phonon  energy  of  30  meV.  The  interference  Integral  is  approx¬ 
imately  0.75.  At  room  temperature,  the  average  transition  probability  is 
estimated  to  be  30  percent.  Once  we  have  a  valley  with  an  angle  less  than 
54.70°,  the  transition  probability  is  greater  than  30  percent,  and  the  n-value 
is  higher  than  0.3  at  room  temperature.  When  we  study  silicon,  a  multivalley 
semiconductor,  the  overall  transition  should  be  taken  into  account.  Thus  the 
numerical  method  has  been  employed  to  calculate  the  attenuation  factor  to  show 
the  effect  of  the  GB  angle. 

The  GB  angle  of  two  neighboring  grains  Is  described  by  three  rotation 
angles  6,  4,  and  <1*.  To  conserve  computer  time,  we  limit  the  rotation  angles 
with  the  conditions  0  <  9  £  45°,  0  <  $  <  45°,  and  $  ■  0  in  the  calculation 
of  the  attenuation  factor  n.  The  results  are  shown  in  Fig.  2(a)  and  (b)  for  n 
at  300  and  77  K.  The  minimum  n-  value  at  I  ■  300  K  is  0.8,  but  the  minimum 
n-value  at  T  ■  77  K  is  0.05.  The  transition  is  influenced  strongly  by  the 
misorlentation  of  the  two  grains  at  low  temperatures,  but  not  at  room 
temperature.  This  is  because  the  phonon-assisted  transition  depends  strongly 
on  the  number  of  phonons  even  though  the  phonon  energy  does  not  vary  with  the 
temperature.  This  is  consistent  with  Einstein's  theory,  which  specifies  the 
phonon  energy  as  the  activation  energy  of  phonon-assisted  transition  when 
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300*  K 


»  kT.  Ue  know  that  we  can  at  leaat  find  a  transition  froa  a  certain 
valley  requiring  a  phonon  energy  leas  than  300  meV.  Thus  the  activation 
energy  is  less  than  30  meV.  Of  course,  this  condition  does  not  apply  when  the 
relaxation  time  is  dominated  by  other  mechanisms  such  as  the  impurity 
scattering.  The  activation  energy  due  to  the  number  of  phonons  remains, 
however,  although  the  transition  is  scaled  down  by  the  low  relaxation  time. 
Generally  speaking,  the  scattering  potential  can  be  explained  by  an  average 
phonon  energy  required  for  the  carrier  transition  between  energy  valleys  in 
different  grains. 

+Xerox  Corporation,  El  Segundo,  CA. 
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D.  DIFFUSIVITY  AND  GROWTH  RATE  OF  SILICON  IN  SOLID-PHASE  EPITAXY  WITH  AN 

ALUMINUM  MEDIUM 

(Q.  Hua,  E.  S.  Yang) 

( JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  E.  S.  Yang  (212)  280-3120 

W.  Hwang  (212)  280-3115) 

Solid-phase  epitaxy  (SPE)  is  a  relatively  new  crystal-growth  process 
whose  importance  should  not  be  underestimated.  The  advantage  of  this 
technique  is  derived  from  the  low  temperature  at  which  single  crystal  layers 
are  grown.  For  example,  when  using  A1  as  the  transport  medium  in  Si  SPE,  the 
growth  temperature  is  between  A00-550°C  so  that  an  abrupt  doping  profile  at 
the  epitaxial  layer-substrate  interface  can  be  easily  obtained.  Moreover,  the 
width  of  the  epitaxial  layer  depends  only  on  the  deposited  amorphous  (a)  Si 
thickness,  which  can  be  controlled  during  evaporating;  thus,  suboicron  thin 
layers  can  be  readily  grown.  Tsaur  et  al.^1^  have  employed  this  technique  for 
the  formation  of  p-n  junction  solar  cells  which  minimize  the  adverse  effects 
of  high  temperature,  and  they  have  produced  devices  with  high  efficiency. 

Since  the  emergence  of  SPE,  a  large  number  of  papers  have  been  published 
on  this  topic.  Most  of  them  describe  the  growth  phenomena  by  using  various 
techniques  such  as  transmission  electron  microscopy,  Rutherford  back- 
scattering  and  channeling,  and  Auger  spectroscopy  to  analyze  the  structural 
character  of  the  growth  layer.  However,  the  number  of  experimental  studies  on 
the  growth  mechanisms  and  related  parameters  is  quite  limited.  In  this  work, 
we  attempt  to  determine  the  dlffuslvlty  and  mass  transfer  coefficient  and  to 
estimate  the  growth  rate  of  SPE. 

The  classical  diffusion  equation  is  used  to  determine  the  diffusivity  of 
A1  in  silicon  in  the  sample  structure  shown  in  the  insert  of  Fig.  1.  When  the 
sample  is  heated  in  a  furnace  the  a-Si  will  gradually  dissolve  into  the 
underlying  A1  layer  until  the  a-Sl  disappears  entirely  from  the  top  surface. 
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Fig.  1  Diffusivity  and  growth  rate  as  a  function  of  temperature 


GROWTH  RATE  R  (ft/sec) 


It  Is  implicitly  assumed  that  the  transport  of  Si  across  the  a-Si-Al  interface 
is  fast.  As  it  will  be  explained  later  this  assumption  is  consistent  with  the 
activation  energy  obtained  from  our  experimental  data.  If  the  growth  process 
is  interface  reaction  limited,  the  activation  energy  would  have  to  be 
approximately  equal  to  the  bond-breaking  energy.  In  other  words,  in  a 
structure  with  a  thin  a-Si  layer  and  relatively  thick  A1  layer,  the  transport 
of  Si  in  A1  may  be  considered  as  diffusion  limited.  Since  we  are  studying  the 
initial  stage  of  silicon  diffusion  up  to  the  point  that  the  surface  silicon 
starts  to  disappear,  we  assume  that  the  interface  silicon  concentration  is 
constant.  The  solution  of  the  diffusion  equation  with  constant  diffusivity  at 
a  given  temperature  is  therefore: 

C(x, t)  -  C^erfcf x/2(Dt)1/2]  (1) 


where  is  the  Si  molar  solubility  in  A1  which  specifies  the  boundary 

condition  at  the  a-Si  and  A1  Interface.  D  is  the  diffusivity  of  Si  In  Al,  x 

is  the  distance  from  the  a-Si-Al  interface,  and  t  is  the  diffusion  time. 

Integrating  Eq.  (1)  from  x  »  0  to  ®  ,  we  obtain  the  relationship  between  the 

total  amount  of  deposited  a-Sl  and  the  diffusivity.  A  similar  study  has  been 

(2 ) 

done  by  Immorllca  et  al.  to  measure  the  diffusivity  and  solubility  of  Ge 
and  GaAs  in  Pb  in  liquid-phase  epitaxy,  and  the  resulting  equation  for  D  is 
modified  as  follows: 


D 


it  /  dSi  pSi  MA1  \2 

'  XSi  pAlMSi/ 


(2) 


where  d^  is  the  thickness  of  a-Si,  Xgi  i3  the  equilibrium  solubility 
expressed  in  atomic  fraction,  which  can  be  obtained  from  the  silicon-aluminum 
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phase  diagram  at  a  given  temperature.  Pg^  and  are  the  densities  of  Si  and 

3  3 

A1  (Pg^  ■  2.33  g/cm  and  p^  •  2.70  g/cm  ).  and  Mg^  are  the  molecular 

weights  of  A1  and  Si.  The  parameter  t  is  the  diffusion  time  which  is 
specified  by  the  disappearance  of  the  a-Si  on  the  top  surface  of  the  sample. 
Because  the  color  of  the  thin  a-Si  film  is  easily  distinguishable  from  A1 
under  an  optical  microscope,  t  can  be  determined  by  observing  the  change  of 
color  of  the  sample. 

The  samples  were  prepared  by  evaporation  in  a  vacuum  of  2  x  10  ^  torr  on 

a  (100)  n-Si  substrate.  The  thickness  of  the  A1  layer  was  1  urn,  and  it  was 

covered  by  ~100  %  of  e-beam  evaporated  a-Si.  The  samples  were  heated  in  a 

furnace  in  a  high-purity  hydrogen  environment.  The  diffuslvlty  D  was 

calculated  from  Eq.  (2)  using  the  measured  t.  The  experimental  results  are 

shown  in  Fig.  1.  Using  the  equation  D  ■  Dq  exp(-EA/kT),  we  find  » 

-7  2 

1-.25  x  10  cm  /sec  and  Z^  “  0.8  eV.  It  is  noted  that  the  value  of  Si 

diffuslvlty  in  A1  in  this  structure  is  several  orders  of  magnitude  less  than 

(3) 

that  in  the  structures  used  for  Al-Si  metallization  '  although  it  is  a  factor 
of  two  greater  than  Tsaur's  results  (Fig.  1).  However,  the  activation  energy 
is  almost  identical  to  that  of  the  Al-Sl  system.  We  believe  that  the 
difference  in  the  values  of  diffuslvlty  is  caused  by  different  experimental 
conditions.  For  example,  a  thin  oxide  layer  at  the  a-Si-Al  Interface  may 
reduce  the  diffusion  process  significantly. 

By  using  the  results  of  the  above  measurement,  another  Important  growth 
parameter,  the  mass  transfer  coefficient  h,  can  be  deduced.  The  process  of 
SPE  may  be  divided  into  three  steps.  In  the  first  step,  the  a-Sl  is  dissolved 
Into  A1  and  then  diffused  through  the  A1  layer  to  ’’each  the  Al-Si  substrate 
interface.  Subsequently,  the  migrated  Si  atom  nucleates  epitaxially  at  the  Si 
substrate  surface  through  a  surface  reaction.  It  has  been  established  that 
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Che  activation  energy  for  surface-reaction-controlled  vapor-phase  epitaxy  is 
(U) 

approx.  1.9  eV.  Therefore,  the  surface  reaction  does  not  seen  to  be  the 

dominant  mechanism  in  SPE  because  of  the  lower  activation  energy  according  to 

(3) 

the  above  measurement  and  previous  experimental  values.  Therefore,  the 
diffusion  flux  in  the  A1  layer  (J^)  can  be  expressed  as  follows: 


J 


1 


-D 


crc2 

dAl 


(3) 


where  is  the  Si  concentration  at  the  Al-Si  substrate  Interface,  and  dA1  is 
the  thickness  of  the  A1  layer.  The  mass  transfer  process  in  the  epitaxial 
growth  can  be  expressed  as: 


J2  *  h(CL  -  C2) 


(4) 


where  J2  is  the  Si  mass  flux.  When  the  equilibrium  growth  condition  is 
reached,  we  have  ■  J2,  so  that 


h 


(5) 


2 

The  total  growth  time  may  be  divided  into  two  parts:  (1)  t^  •  S^/D  *8  t*,e 
time  required  to  reach  the  equilibrium  growth  condition;  (2)  t2  »  dg^/h  is  t*,e 
growth  time  after  the  equilibrium  condition.  In  order  to  compare  our  results 
with  Tsaur's  data  on  growth  rate,  we  have  added  the  vertical  scale  on  the 
right  of  Fig.  1  with  the  specific  conversion  parameters  of  2000  X  A1  and 
2000  X  aSi.^  It  is  evident  that  the  diffusion-limited  mechanism  proposed 
here  is  consistent  with  available  data.  In  practice,  the  growth  time  usually 
should  be  longer  to  produce  higher  quality  growth  layers. 
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In  conclusion,  the  diffusivity  of  a-Si  in  A1  in  the  structure  for  SPE  has 

*  «  Q 

been  measured  between 470-570  C,  and  the  mass  transfer  coefficient  and  growth 
rate  have  been  obtained.  These  data  can  be  used  to  explain  the  results 

reported  by  Tsaur  et  al.^  The  activation  energy  of  this  process  is  found  to 

f  /  3) 

be  0.80  eV,  which  is  consistent  with  the  literature. v  * 

*Thi8  research  was  also  supported  by  the  Solar  Energy  Research  Institute 
contract  XW-1-1272-1 . 
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E.  THERMAL  NITRIDATION  OF  GERMANIUM* 

(Q.  Hua,  J.  Rosenberg,  E.  S.  Yang) 

(JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  E.  S.  Yang  (212)  280-3120 

W.  Hwang  (212)  280-3115) 

The  high  mobility  of  carriers  in  germanium  could  make  it  an  Important 

material  for  realizing  high-speed  digital  or  microwave  metal-insulator- 

semiconductor  field-effect  transistors  if  a  suitable  gate  Insulator  could  be 

fabricated.  Although  germanium  nitride  powder  was  prepared  by  direct  reaction 

of  Ge  powder  with  NH^  as  early  as  1930,^  very  little  Information  exists  in 

(2) 

the  literature  on  thin  germanium  nitride  films.  Igarashl  et  al.  report  the 
growth  of  partially  crystallized  films  by  direct  reaction  of  germanium  wafers 
with  ammonia  (NH^  gas)  above  700°C.  Several  authors^^  report  on  chemical 
vapor  deposition  of  Ge^N^  on  a  variety  of  substrates  (Si,  Ge,  GaAs,  InP)  and 
on  the  excellent  insulating  properties  of  these  films. 

Since  it  has  been  found  that  thermal  nitrldatlon  provides  a  simple  and 
convenient  method  for  fabricating  dense,  uniform,  and  pinhole-free  silicon 
oxynitride  films,^’^  we  have  investigated  the  growth  of  films  by  a  similar 
technique  on  germanium  substrates. 

The  system  used  for  film  growth  is  illustrated  in  Fig.  1.  Nearly 
intrinsic  <111>  germanium  wafers  were  used  as  substrates.  Wafer  preparation 
consisted  of  degreasing  In  organic  solvents  followed  by  a  chemical  etch.  Two 
different  etching  techniques  were  studied  In  order  to  assess  their  impact  on 
the  resultant  film  morphology.  One  etch  consisted  of  a  two-minute  dip  in  30Z 
BjO^,  followed  by  a  deionized  water  rinse  and  immersion  in  buffered  KF  until 
the  wafer  became  hydrophobic.  The  other  etch  was  a  NaOCl  (4Z):  HFtH^O  -  1:1:5 
solution.  The  second  etch  is  much  faster  than  the  first,  and  tends  to  pre¬ 
ferentially  attack  (and  remove)  areas  of  crystal  damage  or  imperfection. 
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Fig.  1  Schematic  diagram  of  the  experimental  system,  (a)  Gas  cylinders, 

(b)  quartz  tube,  (c)  rolling  movable  furnace,  (d)  rotary  pump,  and  (e)  gas  trap. 


Immediately  after  etching,  the  wafers  were  loaded  Into  the  quartz-tube, 
■ovable  furnace  reactor,  which  was  then  evacuated  to  less  than  0.1  Torr.  The 
reactor  was  then  purged  with  high  purity  nitrogen.  After  a  thirty  minute 
purge,  a  mixture  of  NH^  and  Nj  was  Introduced  into  the  reactor  and  the 
preheated  furnace  was  moved  into  place.  After  a  prescribed  growth  time,  the 
furnace  was  removed,  and  the  wafers  were  allowed  to  cool  for  ten  minutes 
before  the  ammonia  flow  was  stopped.  The  reactor  was  purged  with  nitrogen 
before  the  wafers  were  removed. 

Ellipsometry  provides  a  means  for  determining  tiro  of  the  most  fundamental 
parameters  of  a  thin  dielectric  film:  the  thickness  and  the  refractive  index. 

For  our  study,  we  used  HeNe  laser  illumination  ( A  *»  6328  X),  at  an  angle 
of  Incidence  of  70°.  At  the  wavelength  6328  X,  the  complex  index  of  refrac¬ 
tion  of  germanium  is  given  by  5.2  -  0.8i.^^  The  film  was  assumed  to  be 
nonabsorbing  (l.e.,  nf  was  assumed  to  be  real). 

Following  the  procedure  used  by  Archer, v  we  present  in  Fig.  2  a 

universal  chart  for  determining  the  thickness  and  refractive  index  of  any 

nonabsorbing  film  on  germanium.  This  chart  Is  applicable  for  ellipsometry 

systems  utilizing  a  wavelength  of  6328  X  and  an  angle  of  incidence  of  70°. 

The  different  curves  each  correspond  to  a  different  value  of  n^.  Distance  along 

each  curve  gives  the  value  of  6,  which  is  Indicated  by  hash  marks  at  intervals 

of  w/9  (20°)  in  Fig.  2(a).  Film  thickness  is  Indicated  directly  in  Fig. 

2(b).  The  chart  is  used  by  determining  the  point  on  a  curve  which  most 

closely  corresponds  to  the  measured  values  of  ij/  and  A.  This  determines  n^  and 
5>(8), (9) 

Experience  has  shown  that  some  oxygen  Incorporation  in  silicon  nitride  is 
unavoidable  when  the  film  is  grown  thermally.  Consequently,  we  investigated 
oxygen  incorporation  in  our  germanium  nitride  films  utilizing  electron  beam 
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x-ray  .oanalysis.  This  method  also  provides  an  estimate  of  the  film 
thickness  which  is  independent  of  the  optical  measurements.  In  addition, 
carbon  contamination  (from  sources  such  as  organic  solvent  residue  or  pump 
oil)  can  be  Investigated. 

Several  samples  were  analysed  by  this  method.  In  all  cases,  carbon 
contamination  was  at  or  below  the  detection  limit  (a  fraction  of  monolayer). 
The  nitrogen  to  oxygen  ratios  were  in  the  range  of  2:1  to  3:1. 

The  film  thickness  estimates  were  made  by  computing  the  thickness  of  a 
composite  film  incorporating  the  nitrogen  as  Ge^N^  and  the  oxygen  as  GeO^  and 
using  the  bulk  densities  for  these  materials.  Despite  the  crudeness  of  these 
estimates,  they  correspond  to  the  elllpsometry  measurements  within  10X  for 
films  with  smooth  surface  morphology.  It  would  be  difficult  to  make  an 
estimate  of  the  refractive  indices  of  the  films  from  the  compositional 
analysis  since  the  function  dependence  of  the  index  on  composition  is  unknown. 

A  wide  range  of  surface  morphologies  results  from  different  growth 
conditions.  Any  leak  or  contamination  in  the  reactor  gives  rise  to  a  surface 
with  tetrahedral  crystallites  and  whiskers.  Even  under  clean  and  leak-free 
conditions,  the  film  will  tend  to  crystallize  if  either  the  growth 
temperature  la  too  high  or  the  growth  time  is  too  long,  as  observed  for 
silicon  nltrldation.^^  In  our  case,  the  germanium  nitride  crystallites  are 
triangular  and  oriented  along  crystallographic  directions.  Under  appropriate 
growth  conditions  (2-4  h  at  600-650°C  on  a  substrate  which  was  etched  using 
the  treatnent  described  earlier,  and  grown  in  a  1:1  mixture  of  NH^  and 

N^),  very  smooth  and  uniform  films  can  be  grown.  A  preliminary  electron 
diffraction  study  indicates  that  the  smooth  looking  films  are  indeed 
amorphous . 

The  results  of  the  elllpsometry  measurements  are  shorn  in  Figs.  3-5. 
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Fig.  3  Refractive  Index  of  germanium  nitride  films  vs  growth 
time  for  films  grown  In  pure  NH^  at  650°C. 
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Fig.  5  Film  thickness  and  refractive  index  vs  gas  composition  for 
films  grown  at  650°C  for  2h. 
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Fig.  3  shows  the  variation  of  the  refractive  index  of  the  germanium  nitride 

t  films  with  growth  time.  The  average  film  Index  Increases  from  1.5  at  one  hour 

to  1.8  at  four  hours.  Thereafter,  there  is  no  further  increase.  It  should  be 
pointed  out  that  ellipsometric  determination  of  the  refractive  index  of  such 

*  thin  films  is  usually  not  accurate,  so  the  numbers  presented  here  should  be 
regarded  primarily  in  terms  of  the  trend  which  they  indicate.  The  increase  of 
the  refractive  index  with  growth  time  suggests  that  an  initial  oxide  layer 
(which  has  a  lower  index  than  the  nitride)  is  formed  on  the  germanium  surface, 
which  is  subsequently  nitrided  during  the  growth.  Some  corroboration  of  this 
hypothesis  is  provided  by  the  electron  beam  microprobe  analysis  which 
indicates  that  films  which  have  a  lower  refractive  index  have  a  lower  nitrogen 
to  oxygen  ratio.  The  refractive  indices  of  bulk  GeO^  and  Ge^N^  are  1.6  and 
2.1,  respectively. 

*  Fig.  4  shows  film  thickness  as  a  function  of  growth  time  for  films  grown 
in  pure  NH^  on  germanium  at  600  and  650°C  and  on  silicon  at  1200°C.  The  data 
for  the  silicon  growth  curve  is  due  to  Ito  et  al.^  and  is  Included  for 

*  comparison  purposes.  In  all  cases,  the  films  show  rapid  Initial  growth, 
followed  by  a  period  of  very  slow  growth.  The  growth  rate  increases  again 
for  low  growth  times. 

*  Although  the  evidence  is  not  conclusive,  the  fact  that  the  film 
composition  changes  as  a  function  of  growth  time  (i.e.,  the  N:0  ratio 
Increases)  suggests  that  the  growth  mechanism  involves  the  diffusion  of 
nitrogen  and  oxygen  ln*a  and  through  the  film,  which  would  yield  time  varying 
concentrations  of  nitrogen  and  oxygen  in  the  film  and  at  the  film-substrate 
Interface.  If  the  growth  mechanism  involved  the  diffusion  of  germanium  to  the 
surface  of  the  film,  the  reacting  germanium  would  be  exposed  to  an  essentially 
time  invariant  environment,  and  one  would  expect  the  film  composition  to  be 
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essentially  independent  of  growth  time.  We  propose  that  the  initial  growth  is 
surface-reaction  rate  limited.  In  this  regime,  there  may  be  significant  oxide 
formation  from  residual  oxygen  (as^both  0 2  and  water)  in  the  reactor.  As  the 
film  becomes  thicker  and  the  rate  of  nitrogen  and  oxygen  diffusion  through  the 
film  becomes  the  limiting  factor  in  the  growth  rate,  nitride  formation  will 
most  likely  dominate  due  to  the  relative  concentrations  of  nitrogen  and 
oxygen.  Some  replacement  of  oxygen  by  nitrogen  in  the  film  probably  occurs  as 
well.  Germanium  nitride  is  substantially  denser  than  germanium  oxide,  and  in 
the  form  of  a  continuous  amorphous  film,  it  probably  presents  a  strong  barrier 
against  the  diffusion  of  both  oxygen  and  nitrogen  (as  does  silicon  nitride). 
This  would  account  for  the  regime  of  very  slow  growth.  The  onset  of  more 
rapid  growth  at  long  growth  times  corresponds  roughly  to  the  onset  of 
crystallization  of  the  films.  On  the  basis  of  this,  we  believe  that  struc¬ 
tural  changes  which  accompany  crystallization  of  the  films  allow  Increased 
diffusion  of  nitrogen  through  the  film  (perhaps  along  grain  boundaries),  which 
results  in  the  Increase  of  the  growth  rates. 

Film  growth  in  mixtures  of  and  NH^  was  also  studied.  Figure  5  shows 
film  thickness  and  refractive  index  as  a  function  of  mixture  composition  for 
two  hours  growth  at  650°C.  Several  authors  have  reported  growing  nitride 
films  on  silicon  wafers  in  pure  Ng,  but  we  observed  no  measurable  growth  on 
germanium  in  pure  Nj  at  temperatures  up  to  650°C.  This  is  consistent  with 
Johnson*  s^^  results  on  reaction  of  germanium  powder  with  •  The  absence  of 
nitride  formation  on  germanium  in  pure  probably  is  due  to  insufficient 
generation  of  dissociated  nitrogen  from  at  a  growth  temperature  of  650°C 
(as  contrasted  with  1200  to  1300°C  used  for  formation  of  silicon  nitride  from 
pure  N2). 

Several  nitride  films  were  grown  in  pure  NR.  for  greater  than  4  h,  which 


had  thicknesses  greater  than  100R  and  refractive  indices  greater  than  2. 
Unfortunately,  results  from  these  thicker  filas  wer*  not  reproducible  and  are 
not  Included  in  Figs.  3-5. 

Several  saaples  vere  prepared  for  electrical  aeasureaents  by  evaporating 
an  ohaic  contact  on  the  wafer  backside  (Sn)  and  1-aa  aluminum  dots  on  the 
fila.  The  electric  field  which  the  films  can  sustain  varies  widely  with 

surface  morphology.  Films  with  smooth,  noncrystallized  surfaces  show  the 

—  13  2  6 

onset  of  substantial  leakage  (10  A/ pm )  at  fields  of  the  order  of  10  V/cm, 

and  destructive  breakdown  at  fields  in  excess  of  3  x  10**  V/cm. 

Germanium  nitride  films  with  thicknesses  less  than  150  X  can  be  grown  by 

thermal  reaction  with  ammonia.  The  best  surface  morphology  is  obtained  using 

NH^N^  mixtures  with  ratios  in  the  range  of  1:1  to  2:1,  with  growth 

temperatures  between  600  and  650°C.  Reactor  cleanliness,  gas  purity,  and 

wafer  surface  preparation  all  have  significant  impact  on  the  resulting  film 

morphology. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  no.  ECS-81-16163. 
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F.  GRAIN  BOUNDARY  INTERFACE  STATES  MEASUREMENT  BY  A  LASER  BEAM* 

(E.  Poon,  E.  S.  Yang,  H.  L.  Evans,  W.  Hwang,  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  4,  1982-1985) 

(Principal  Investigators:  E.  S.  Yang  (212)  280-3120 

W.  Hwang  (212)  280-3115) 

Traditionally,  polycrystalline  materials  have  been  used  in  Integrated 

circuit  interconnections,^*^  varistors,  and  positive  temperature  coefficient 
(2) 

devices.  In  the  last  few  years,  new  interest  has  been  centered  on 

(3) 

recrystallized  semiconductors  for  microelectronics.  The  earlier 
applications  made  use  of  the  unusual  electrical  properties  derived  from  grain 
boundaries  (GB's)  whereas  the  new  research  treats  the  GB  as  a  structural 
defect  giving  rise  to  poor  electrical  performance.  In  either  case,  an 
understanding  of  the  GB  is  of  fundamental  importance  for  device  applications. 
The  GB  barrier  height  is  the  most  important  parameter  in  determining  the 
conductivity.  Until  now,  because  of  the  changing  structural  properties  along 

(  A) 

the  GB  between  two  adjacent  grains,  7  accurate  measurement  of  the  barrier 
height  has  not  been  possible.  For  example,  large  area  capacitance  and 
conductance  techniques  can  only  yield  electrical  parameters  averaged  over  an 
uncertain  cross-sectional  area.  In  this  report  we  present  a  method  to 
overcome  these  difficulties. 

The  schematic  diagram  of  the  experimental  setup  is  shown  in  Fig.  1.  The 
chopped  He-Ne  laser  beam  is  focused  on  a  grain  boundary  of  a  p-type  poly¬ 
crystalline  silicon  sample  which  rests  on  the  X-Y  state.  The  sample,  approxi¬ 
mately  2  x  10  x  0.2  mm,  contains  two  electrically  active  parallel  grain 
boundaries  that  run  across  the  whole  width.  The  external  dc  bias  V_  is  kept 

O 

constant  throughout  the  experiment.  The  chopped  light  beam  generates  a 
periodic  photoconductance  signal  which  is  measured  across  the  load  resistor 
R^  by  a  boxcar  averager.  A  graph  of  photoconductance  versus  temperature  is 
then  plotted  by  an  X-Y  recorder  in  real  time. 


63 


OBJECTIVE 

CHOPPER  LENS 


Fig.  1  Schematic  diagram  of  the  laser-induced  photoconductivity 
experiment. 


The  photoconductance  signal  AV^  across  the  load  resistor  is  given  by 


VkT 

AVL  “  AA*qT  6 


(1) 


if  the  intensity  of  the  incident  beam  is  high  enough  to  reach  the  flatband 
condition.  Since  the  pre-exponential  factor  is  relatively  insensitive  to 
temperature,  the  slope  of  the  plot  of  In(AV^)  vs  1/T  gives  the  barrier  height 
in  the  dark  4*^* 

The  laser  Intensity  under  typical  conditions  is  estimated  to  be 

12.7  x  10^  mW/cm^  (4  x  10*®  photons/cm^  s).  is  150  kft  and  is  set  at  1.5 

V.  The  sample  has  a  dark  resistance  of  2  kfi.  From  the  C-V  measurement,  the 

15  -3 

doping  concentration  is  found  to  be  2  x  10  cm  .  It  is  apparent  that  the 

flatband  condition  is  established  at  this  intensity  since  the  signal 

saturates.  The  signal  starts  to  decrease  when  the  laser  Intensity  is  at  a 

value  equal  to  1/5  of  the  typical  intensity.  Figure  2  shows  the  plot  of 

ln(AV  )  vs  1000/T  at  the  typical  Intensity.  The  barrier  height  is  found  to  be 
L 

0.37  eV. 

A  typical  decay  time  is  shown  in  Fig.  3  in  which  two  time  constants  are 

seen.  Our  experimental  setup  limits  our  ability  to  measure  the  fast  time 

constant  so  that  our  analysis  is  based  on  the  slower  decay  time  (5.44  ms)  in 

Fig.  3.  The  photoconducting  decay  time  (t.)  is  related  to  the  trap  emission 

d 

rate  (eQ)  by 


where  v  ^  is  the  thermal  velocity,  oq  is  the  thermal  capture  cross  section,  N 
is  the  effective  density  of  states  of  the  conduction  band,  is  the 
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conduction-band-edge,  is  the  trap  energy,  and  T  is  the  temperature.  Using 
Eq.  (2)  and  the  data  in  Fig.  4  we  obtain  an  activation  energy  of  0.5  eV. 

Since  v  -  107  cm/s  and  N  *  lO^-9  cm  ^  at  300°K,  we  find  a  ■  2  x  10  ^  cm^, 
which  is  a  reasonable  value  for  a  deep  trap  in  silicon. 

From  these  experimental  data,  we  have  obtained  a  dominant  trap  of  the 
donor  type  at  0.5  V  from  the  conduction-band  edge.  In  a  thermal  activation 
experiment  performed  recently  Redfleld  found  that  the  GB  potential  barriers 
increase  markedly  after  annealing  at  900  -  1000°C.^^  He  concluded  that 
oxygen  Impurity  at  the  GB's  is  the  only  acceptable  explanation  because  oxygen 
requires  lov-temperature  exposure  before  segregation  at  defects  occurs. 

The  known  energy  level  of  oxygen  in  silicon  is  a  donor  level  at  0.51  eV  below 
the  conduction-band  edge,^7^  in  remarkable  agreement  with  our  measured  value. 
The  oxygen  segregation  at  the  GB  would  also  explain  why  the  measured  GB 
interface  states  are  not  continuous  as  in  the  case  of  SiO^-Sl  Interface.  As 
for  the  barrier  height,  our  result  is  in  the  range  determined  from  earlier 
experiments. 

An  Important  feature  of  the  present  technique  is  its  ability  to  measure  a 
localized  region  of  defects.  In  contrast  with  the  capacitance  and  conductance 
methods  which  are  capable  of  measuring  materials  with  a  single  GB  average  over 
a  large  area,  the  focused  laser  beam  can  probe  individual  GB  in  a  wafer  with 
many  grains  and  its  spatial  resolution  is  limited  only  by  the  light  beam.  In 
comparison  with  the  electron-beam-induced  current  (EBIC)  technique,  the 
present  approach  does  not  require  a  p-n  junction  so  that  materials  can  be 
measured  without  going  through  a  high-temperature  process.  Thus,  Influence  of 
device  processing  can  be  determined  before  and  after  a  particular  step. 
Further,  the  light  absorption  mechanism  is  well  characterized  whereas  the 
electron  beam  impact  on  semiconductor  remains  to  be  clarified. 
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In  conclusion,  the  method  Is  simple  and  has  good  spatial  resolution  which 


is  limited  by  the  region  activated  by  the  focused  laser  beam.  By  observing 
the  steady  state  transient  response,  we  show  for  the  first  time  that  the 
photoconductance  can  be  used  to  calculate  the  grain  boundary  trapping  states 
parameters  in  an  unambiguous  way. 
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G.  ETCHING  OF  SINGLE  CRYSTAL  GaAs  SUBSTRATES  BY  UV  LASER-INITIATED  GAS-PHASE 
PHOTODISSOCIATIONS  OF  CF^Br  AND  CH^Br 

(P.  Brewer,  S.  Halle,  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.  M.  Osgood,  Jr.  (212)  280-4462) 

The  rational  behind  the  use  of  uv-laser,  gas-phase  photodissociation  In 
the  heterogeneous  etching  of  semiconductor  substrates  lies  in  the  fact  that 
the  initiation  of  the  gas-solid  reactions  occurs  only  after  the  absorption  of 
a  uv  photon  by  an  inert  parent  molecule.  Various  inert-halogen-contalnlng 
parent  molecules  have  been  employed  to  release  reactive  halogen  atom  fragments 
by  uv-photons  for  etching  of  amorphous  GaAs  and  InP.^^  In  these  experiments 
due  to  the  laser-substrate  configuration,  the  distinction  between  pure  photo¬ 
chemical  reactions  and  laser-surface  interactions  could  not  be  totally 
separated.  Another  approach  to  ga8-phase  etching  has  been  in  the  use  of 

vlbratlonally  excited  molecules  containing  fluorine  (F)  atoms  in  etching  Si 

(2\  (3)  2  o  2  o 

substrates.  '  In  the  present  study  bromine  atoms  (  p  and  P 3/2^  are 

used  in  the  etching  of  single  crystal  GaAs  by  uv  laser  initiated  photo- 

dissociation  of  CF^Br  and  CH^Br. 

The  nature  of  the  photodlssoclatlon  of  CF^Br  and  CH^Br  in  the  uv  region 
(4) 

has  been  summarized.  The  two  primary  processes  which  occur  in  this  region 
are 


CX^Br  +  hv 

-► 

CX3  +  Br(2P1/2) 

X  -  F  or  H 

CXjBr  +  hv 

-► 

CX3  +  Br(2P3/2) 

Leone  et  al.^  have  found  that  the  quantum  yield  of  the  electronically 
2  o 

excited  Br(  P^2)  «toms  produced  in  the  photodlssoclatlon  at  193  nm  of  CF^Br 
to  be  0.50.  In  the  case  of  CR^Br  the  indrease  of  the  population  as  a 

function  of  exposure  time  at  193  nm  was  measured  with  an  initial  pressure  of 


I 
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10  torr  of  CF^Br.  Shown  in  Fig.  1  is  a  plot  of  the  CjF^  signal  with  exposure 
tiae  indicating  a  significant  portion  of  the  CF^  radicals  are  eliminated  by 
blaolecular  collisions. 


cf3  +  cf3 


C2P6 


Early  work  on  gas-phase  etching  of  Si  substrates  by  ir-aultlphoton 

dissociation  of  CF3Br  indicated  that  the  heterogeneous  gas-solid  reaction 

(2) 

causing  etching  was  the  vibrationally  excited  CF3  radicals.  The  evidence 
from  the  present  experiments  indicate  that  the  etching  reactions  involve  the 
Br  atoas  with  the  GaAs  substrate.  This  assumption  is  based  on  the  similar 
reactivities  of  the  photodissociate  fragments  froa  both  CFjBr  and  CH3Br  and 
the  likely  dissimilar  reactivities  of  CH3  and  CF3  radicals. 

The  products  of  the  gas-surface  reaction  of  bromine  atoms  and  the  GaAs 
substrate  fora  liquids  on  the  substrate  at  rooa  temperature.  The  effect  of 
the  liquid  is  to  coat  the  surface  with  a  protective  film  as  seen  in  Fig.  2 
where  the  etch  depth  as  a  function  of  laser  power  is  plotted.  In  order  to 
eliminate  the  liquid  formation  or  the  GaAs  substrates  the  samples  were  raised 
to  T  ^  100°C. 

Two  types  of  experiments  were  performed  with  respect  to  the  laser- 
substrate  configuration,  positioning  in  the  parallel  (indirect/  and  perpen¬ 
dicular  (direct)  orientations.  The  parallel  configuration  used  a  focused 
exciaer  laser  beam  positioned  between  0.1  tai  and  5  am  above  the  GaAs  surface. 
The  etching  rate  was  measured  by  monitoring  the  Ga  atomic  population  in  the 
gas  phase  using  a  Ga  fluorescence  line  at  417  na.  The  detection  scheme 
employed  was  based  on  the  assumption  that  the  focused  exciaer  laser  beam  after 
generating  the  reactive  broaine  atoas  with  one  pulse  would  probe  (in  the  next 
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pulse)  all  Gs  and  As  containing  con pounds  leaving  the  surface  entering  the  gas 
phase  by  nultiphoton-dlsaoclatlng  the  product  compound  and  exciting  the  Ga  or 
As  atoms  produced.  Thus  by  using  a  static  cell,  a  “steady  state"  population 
of  Ga  or  As  atoms  could  be  monitored.  In  Fig.  3  the  population  of  Ga  atoms  in 
the  gas  phase  is  measured  as  a  function  of  laser  exposure  time.  Fig.  4  shows 
the  increase  of  Ga  “steady  state"  population  with  temperature  giving  an 
activation  energy  of  14.9  Real /mole.  The  etched  atom  population  was  insensi¬ 
tive  to  the  position  of  the  laser  beam  over  the  GaAs  surface.  This  set  of 
data  confirms  the  "cloud  or  diffusion"  model  for  the  etching  reaction.  The 
etch  profile  of  a  GaAs  (100)  surface  is  characterized  by  a  rough,  sawtooth 
pattern  indicative  of  anisotropic  effects.  Typical  conditions  for  the 

indirect  (parallel)  configuration  were:  surface  temperature  "•60oC;  laser  20Hz, 
2 

350mJ/cm  (focused),  15  min  exposure;  and  7  torr  of  CF^Br. 

The  direct  (perpendicular)  configuration  produced  a  much  more  uniform 

etch  with  smaller  microscopic  structures  (grain  size) .  The  experimental 

details  for  a  typical  direct  etch  are  as  follows:  surface  temperature  ~60°C; 

2 

laser  50Hz,  35mJ/cm  (mild  focus),  5-10  min  exposure;  and  7  torr  of  CF^Br. 

The  etching  of  the  GaAs  was  confined  mainly  to  the  illuminated  area  of  the 
sample.  Fig.  5  shows  a  plot  of  the  etch  depth  as  a  function  of  exposure  time 
(laser  flux).  Strict  control  of  the  laser  power  was  required  in  order  to 
avoid  any  laser-surface  etching  by  oblation.  For  long  etching  times  >10  mins 
the  non-llluminated  parts  of  the  substrate  show  etch  pits. 

This  report  shows  the  promise  of  gas  phase  etching  of  single  crystal  GaAs 
by  uv-laser  initiated  photochemistry.  The  method  of  preference  is  the  direct 
(perpendicular)  configuration  due  to  the  small  surface  microstructure  and  its 
applicability  to  large  area  etching  of  GaAs. 
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Go  Atom  Populotion  vs.  Loser  Exposure  Time 
GaAs  T  =  60°C/  Laser  Indirect  20  Hz,  350mJ/cm2) 


Fig.  3  Ga  atom  steady  state  population  in  gas  phase. 
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Figure 


Exclmer  Loser  (I93nm)  Photochemical 
Large  Scale  Direct  Etching  of  GoAs 

Gas:  CF5Br,  7  Torr 


78 


EXPOSURE  TIME  (min.) 

Figure  5 
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H.  DIRECT  OBSERVATION  OF  A  SURFACE-ENHANCED  CHEMICAL  REACTION* 

(C.J.  Chen,  R.  M.  Osgood,  Jr.) 

(JSEF  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.  M.  Osgood,  Jr.  (212)  280-4462) 

Ue  have  recently  published  a  theoretical  study  of  the  optical  field  near 
the  surface  of  microscopic,  ~200  X,  conducting  spheres.  The  calculations 
showed  that  the  overall  electric  field  is  strongly  enhanced  if  the  optical 
frequency  is  close  to  a  resonant  frequency  of  the  sphere.  In  addition,  the 
electric  field  is  strongest  at  the  poles  of  the  sphere,  where  the  orientation 
of  the  poles  is  determined  by  the  polarization  vector  of  the  light;  see  Fig. 

I.  Both  the  overall  and  directional  enhancements  are  dependent  on  the  shape, 
size,  and  composition  of  the  particle.  The  growth  of  the  particle  by 
photodeposition  is  dependent  on  the  electric  field  at  its  surface.  Thus, 
there  is  a  coupling  between  the  particles'  shape  and  size  and  the  photo- 
dlssoclatlng  electric  field.  The  theoretical  study  indicates  that  a  spherical 
particle  should  grow  most  rapidly  at  its  poles,  transforming  itself  into  an 
ellipsoid  in  the  process.  The  final  shape  of  the  particle  for  a  relatively 
short  growth  period  will  depend  on  how  close  the  particle's  resonant  frequency 
matches  that  of  the  optical  field  during  growth.  The  final  size  of  the 
particle  is  governed  not  only  by  the  match  in  frequencies,  but  also  by  the 
optical  losses  for  both  very  large,  500-1000X,  and  very  small,  0-50  X, 
particles. 

In  order  to  observe  these  effects,  we  used  a  collinear  observing  and 
illumination  microscope,  which  allowed  us  to  accurately  position  the  sample 
and  observe  the  deposited  film  thickness.  With  this  system  we  observed 
deposition  of  Cd  from  Cd^H^g  molecules  adsorbed  on  Cd  spheroids,  10  -  500X 
in  diameter.  These  small  metal  spheres  had  a  resonant  frequency  close  to  the 
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optical  frequency  of  the  UV  laaer  at  257  na. 

The  metal  spheres  were  made  by  vacuum  depositing  the  metal  In  an  evapora¬ 
tor  filled  with  a  relatively  high  pressure  of  argon  to  encourage  gas-phase 
nucleatlon.  The  spheroids  were  round  for  small  diameters,  but  lnclplently 
crystalline  for  large  diameter;  see  Fig.  2.  Prior  to  Irradiation  with  the  UV 
laser  light,  the  spheres,  which  were  mounted  on  a  thin  carbon  film,  were 
positioned  on  a  fine  copper  mesh  for  examination  with  a  TEM.  This  sample  was 
then  placed  In  a  cell  chamber  and  exposed  to  sufficent  CdCCHj)^  gas  to  form  an 
adsorbed  film  of  these  molecules  on  the  spheres. 

The  results  of  the  irradiation  are  shown  In  Fig.  3.  First,  It  Is 
apparent  that,  as  predicted,  the  small  Cd  spheroids  grow  into  ellipsoids. 

This  growth  pattern  Is  In  accord  with  the  high  electric  field  at  the  poles  of 
the  particles.  Second,  although  not  clearly  shown  In  the  figure,  the  smaller 
particles  do  not  exhibit  rapid  growth  until  after  they  grow  by  a  much  slower 
process  to  a  critical  size  of  -50  X.  This  latency  period  Is  In  accord  with 
the  optical  losses  of  the  small  particles.  Similarly  we  find  that  very  large 
particles,  which  are  lossy  because  of  reradiated  light,  do  not  grow  rapidly. 
Pig.  4  shows  a  recent  calculation  of  this  effect  based  on  the  theoretical 
model  described  above. 

In  this  experiment,  we  have  studied  film  growth  from  separately  fabri¬ 
cated  tiny  spheres.  However,  It  Is  clear  that  the  results  obtained  are  also 
applicable  to  the  evaluation  of  small  nuclei  Into  macroscopic  grains  during 
photodeposition. 

*This  research  was  also  supported  by  the  Army  Research  Office  under  Contract 
DAAG29-82-R-0009 . 


I.  LASER-CONTROLLED  PLASMA  ETCHING* 

(G.  Reksten,  W.  Holber,  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.  M.  Osgood,  Jr.  (212)  280-4462) 

Introduction  of  light  into  a  semiconductor  processing  chamber  allows  one 
to  control  the  spatial  extent  and  chemical  reactions  in  the  chamber.  One  can 
envision  experiments  with  laser  controlled  CVD,  MBE,  etc.  We  have  recently 
begun  experiments  to  determine  the  extent  to  which  visible  laser  light  can  be 
used  to  control  etching  of  semiconductor  wafers  in  a  plasma  etching  machine. 

The  enhancement  can  be  based  on  surface  or  gas-phase  heating,  a  laser-beam 
generated  photovoltage,  local  ionization,  or  surface  desorption. 

A  commercial  plasma  etcher  (Technics  Model  PE  II)  was  used  in  the  experi¬ 
ments  described  below.  The  etch  rate  of  n-type  phosphorous  doped  Si  In  a  CF^ 
plasma  was  measured  as  a  function  of  CF^  pressure  and  r.f.  power.  A  2W  argon 
ion  laser  (C.W.)  was  then  used  in  conjunction  with  the  plasma  etcher  and  the 
effect  of  the  laser  beam  on  the  etching  process  was  investigated. 

a)  Experimental  arrangement 

A  schematic  of  the  plasma  etcher  and  the  optical  path  of  the  Ar+  laser 
beam  arrangement  is  shown  in  Fig.  1.  The  r.f.  power  applied  to  the  electrodes 
could  be  varied  from  0-300  W  and  the  r.f.  frequency  of  this  particular  model 
was  30  kHz.  The  wavelength  of  the  cw  argon  ion  laser  was  514  nm  (green)  and 
the  maximum  power  available  was  2  W. 

In  order  to  find  the  optimum  CF^  pressure,  several  Si  samples  were  etched 
at  CF^  pressures  ranging  from  50-500  mTorr.  The  r.f.  power  was  kept  constant 
at  180  W,  an  arbitrarily  chosen  value.  The  resulting  curve  of  etch  rate 
versus  pressure  is  shown  in  Fig.  2.  The  maximum  etch  rate  was  obtained  at  a 
CF^  partial  pressure  of  100  mTorr. 
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b)  Etch  rate  vs.  r.f.  power 

The  above  experiaent  was  repeated,  this  tlae  keeping  the  CF^  pressure 
constant  at  100  aTorr  and  ranging  the  r.f.  power  froa  50-300  W.  The  resulting 
etch  rate  versus  r.f.  power  curve  is  shown  in  Fig.  3.  At  powers  >  100W  the 
etch  rate  increases  linearly  with  r.f.  power. 

c)  Laser-plasma  etching 

The  effect  of  the  argon  ion  laser  beam  on  the  Si  etch  rate  was 
investigated  using  the  experimental  arrangement  shown  in  Fig.  1.  The  laser 
beam  was  focused  onto  the  Si  sample  with  the  aid  of  a  beam  collimator  and  a 
lens  placed  inside  the  etching  chamber,  but  outside  the  plasma.  Ideally  the 
angle  of  incidence  should  be  closer  to  90°,  this  was  however  not  feasible  due 
to  the  geometry  of  the  Technics  plasma  etcher.  For  every  measurement  made,  a 
second  Si  sample  was  placed  near  the  laser-irradiated  sample  for  comparison 
purposes.  The  Si  wafers  were  masked  with  an  aluminum  dot  pattern  (1.5  pm 
thick.)  This  pattern  stayed  Intact  after  the  etching  process. 

Twenty  minute  runs  were  then  made  with  the  r.f.  power  constant  at  180  W 
and  the  CF^  pressure  constant  at  100  mTorr  -  while  the  laser  power  was  varied 
from  0.25-2.0  watts.  The  aluminum  mask  was  then  removed  and  the  etch  depth 
measured  on  both  the  irradiated  and  unirradiated  sample. 

The  data  resulting  from  this  experiment  is  plotted  in  Fig.  4.  The  plot 
shows  that  at  high  laser  powers  the  etch  rate  increases  dramatically  due  to 
the  laser.  It  should  be  noted  that  the  laser  enhances  the  etch  rate  over  the 
entire  Si  sample  (7  mm  X  7  mm)  even  though  the  laser  was  focused  to  a  spot  of 
1  am  diameter.  Further  experiments  are  needed  to  explain  this  and  to  obtain 
more  data  for  a  more  complete  etch  rate  versus  laser  power  curve. 
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A  second  "bulk  effect"  was  also  observed*  At  high  laser  powers  (above 
1W)  the  surface  roughness  over  the  entire  etched  SI  piece  decreased  markedly 


Surface  roughness  of  SI  etched  with  CF^  has  been  reported  by  other  workers 
when  reactive  Ion  etching  Is  employed  and  from  our  results  we  see  that  laser 
irradiation  may  have  a  positive  effect  on  this  problem. 

At  high  laser  powers,  an  etch  fit  of  the  focused  laser  could  be  seen. 

The  microscope  revealed  that  this  mask  had  a  fringe  structure  which  serves  to 
Illustrate  the  potential  resolution  of  this  etching  method. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research/DARPA  under  Grant  F-49620-82-K-0008 . 
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J.  OPTICAL  TECHNIQUES  FOR  INSULATOR  FORMATION* 

(R.  Krchnavek,  H.  H.  Gllgen,  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.  M.  Osgood,  Jr.  (212)  280-4462) 

Low- temperature  oxide  formation  is  a  desirable  capability  in  planar 
processing  of  microelectronics.  The  normal  high- temperature  growth  used  in 
standard  silicon  processing  results  in  the  bowing  of  silicon  wafers  and 
prevents  the  use  of  photolithographic  lift-off  techniques  for  oxide  pat¬ 
terning.  For  direct  writing,  low-temperature  formation  is  also  an  important 
capability  because  it  prevents  thermal  alteration  or  damage  In  a  region 
adjacent  to  that  being  processed. 

He  have  recently  demonstrated  that  it  is  possible  to  use  a  pulsed  laser 
to  form  an  SiO^  structure  from  a  thin  layer  of  silicate  which  is  spun  on  the 
wafer.  This  approach  was  begun  in  collaboration  with  Allied  Corporation  in 
Morristown,  NJ,  which  provided  the  silicate.  The  technique  works  by  spinning 
on  the  silicate  which  is  solvated  in  an  organic  liquid;  a  patterned  oxide  film 
is  then  produced  by  exposing  the  film  to  approximately  0.3  J  of  UV  light  at 
193  nm. 

The  writing  apparatus  used  was  that  described  in  the  previous  section. 

The  quality  of  the  oxide  line  pattern  depends  sensitively  on  the  laser  power 
and  techniques  used  for  removing  the  untreated  spin-on  material  from  the  wafer. 

Infrared  spectroscopy  has  been  used  to  study  the  oxide  layers  produced  on 
single-crystalline  silicon  using  a  spin-on  silicate  subsequently  irradiated 
with  an  ArF  pulsed  exclmer  laser  (193  nm).  Figure  1  shows  the  IR  spectra  of 
thermally  prepared  silicates. 

Curve  A  is  obtained  with  the  fresh  silicate  as  spun  onto  an  Si  wafer;  the 
reference  signal  is  from  a  bare  Si  wafer.  Curve  B  shows  the  coating  after  a 
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150°C-15  min  prebake;  it  ahovs  a  reduction  in  the  CH~  and  CH_  Uganda  which  la 
t  1  i 

due  primarily  to  volatlzlng  of  the  aolvent.  After  baking  in  an  oven.  Curve  C, 

the  CH^  and  CH^  peaks  are  completely  gone.  Ue  also  see  a  decrease  in  the  SiOH 

(allanol)  peak  although  it  is  somewhat  hidden  because  of  the  S10  bonds  in 

i 

curves  A+B. 

Fig.  2  shows  the  results  after  laser  irradiation  (curve  A)  and  laser 
irradiation  with  subsequent  baking  (curve  B).  Close  examination  of  these 
curves  shows  that  the  excimer  irradiation  acts  as  a  low  temperature  cure. 
Repeated  exposure  times  did  not  change  the  spectrum.  Further  work  will 
Include  a  measurement  of  higher  UV  exposures. 

Breakdown  voltage*  tests  on  irradiated  samples  show  values  similar  to 
thermally  prepared  organosilicate  samples  (1-2  10^  V/cm)  and  slightly  less 
than  the  published  value  for  thermally  grown  Si02  (107  V/cm). 

*Thi8  research  was  supported  by  the  Air  Force  Office  of  Scientific 
Research/DARPA  under  Grant  F-49620-82-K-0008. 

» 

» 

> 


t 


93 


SILABOL  SOLVENT 


K.  MASKLESS  FABRICATION  OF  SUBMICROMETER  GRATINGS  * 

(D.  Podlesnlk,  H.  H.  Gilgen,  «nd  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  5,  L982-I98S) 

(Principal  Investigator:  R.  M.  Osgood,  Jr.  (212)  280-4462) 

Recent  Interest  in  guided-wave  optics  and  distributed  feedback  lasers  has 
prompted  the  need  for  producing  gratings  with  periods  on  the  order  of  100  nm 
on  semiconductor  materials.  Other  applications  include  surface  acoustic  vave 
devices  and  negative-resistance  electron  devices.  Currently,  there  are  two 
conventional  methods  for  grating  fabrication:  holographic  lithography  and 
electron-beam  lithography.  Uith  the  former  method,  gratings  with  spacing  as 
small  as  110  nm  have  been  produced;  however,  this  technique  is  cumbersome,  and 
grating  periods  are  generally  restricted  to  be  greater  than  180  nm.  Uith  the 
electron-beam  technique,  gratings  with  periods  less  than  100  nm  can  be 
produced.  But,  this  method  is  limited  by  the  difficulties  of  producing 
uniform  large-area  structures  and  the  long  exposure  times,  which,  for  a  300  nm 
grating,  can  be  as  long  as  one  hour  per  square  millimeter.  Both  methods 
Involve  multiple  steps.  In  addition,  for  both  methods,  photoresist  residues 
cause  contamination  problems,  limiting  their  use  in  electro-optical 
applications. 

In  the  present  paper  we  describe  a  technique  for  producing  submicrometer 
gratings  by  a  holographic  method  with  direct  chemical  etching  in  an  oxidizing 
etch.  The  resolution  of  this  technique  is  comparable  to  those  of  other 
methods.  Furthermore,  this  method  is  a  single,  one-step  process  using  very 
low  laser  Intensity.  Both  uniform  large-area  gratings  and  special  grating 
profiles,  e.g.  blazed,  can  be  easily  fabricated. 

The  anodic  dissolution  of  a  semiconductor  in  an  oxidizing  etch  la 
generally  Influenced  by  the  number  of  carriers  present  at  the  surface.  Since 
light  incident  on  the  semiconductor  surface  produces  electron-hole  pairs. 
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a  pattern  of  light  produces  a  corresponding  pattern  of  available  carriers. 

The  photogenerated  carriers,  typically  holes  for  n-type  GaAs,  control  the 
etching  rate  of  the  sealconductor,  producing  an  engraved  pattern  in  Its 
surface.  The  resolution  of  the  etching  process  nay  be  shorn  to  be  a  direct 
consequence  of  rapid  nomal  electric  field,  due  to  the  sealconductor  band 
bending  at  the  Interface.  Our  results  show  that  170  na  gratings  are  easily 
obtained  and  we  anticipate  that  with  an  appropriate  optical  arrangeaent,  a 
resolution  below  100  na  is  possible. 

a)  Experiment 

The  actual  experlaental  arrangeaent  used  for  the  holographic  exposure  Is 
shown  in  Fig.  1.  A  laser  beam  froa  an  argon-ion  laser  emitting  at  514.5, 
488.0,  or  457.9  na  passes  through  a  spatial  filter.  The  filtered  bean  is 
incident  on  a  right  angle  corner  on  which  a  mirror  and  a  sealconductor  sample 
are  aounted.  The  corner  is  Inside  an  optical  cell  containing  an  etching 
solution.  The  Interference  pattern  is  produced  by  the  superposition  of  the 
direct  and  the  reflected  beams  inside  the  solution.  The  grating  period,  s,  is 
given  by 

X 

*"  2n  sin9^ 

where  A  is  the  free-space  wavelength  of  incident  light,  n  is  the  index  of 
refraction  of  the  solution  and  6^  is  the  incident  angle  on  the  sample 
surface.  The  incident:  light  is  s-polarized .  After  passing  through  the 
spatial  filter,  the  laser  beaa  goes  through  the  window  of  the  optical  cell, 
and  it  is,  therefore,  desirable  to  aatch  as  much  as  possible  the  indices  of 
refraction  of  the  solution  and  the  window  in  order  to  reduce  spurious 
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TOP  VIEW: 


ETCHING  SOLUTION 


Fig.  1  Experimental  arrangement  for  grating  fabrication 
by  direct  chemical  etching. 


Interference  patterns. 

The  experiments  were  done  on  various  GaAs  crystals.  A  comparison  was 

18  ~7 

made  of  etching  gratings  on  n-type  (10  cm  )  and  semi-insulating  GaAs  of 
(100)  and  (111)  orientations.  On  all  samples,  good  quality  and  small-period 
gratings  were  produced.  Figure  2  shows  an  SEM  photograph  of  200  nm  grating  on 
n-type  GaAs.  Using  the  457 .9  nm  line  we  have  recently  obtained  gratings  with 
a  170  nm  period.  In  a  separate  experiment,  we  obtained  the  blazed  grating 
shorn  in  Fig.  3.  The  approach  was  to  tilt  the  substrate  so  that  the  bisector 
of  the  angle  of  interference  was  at  an  angle  with  respect  to  the  surface 
normal.  This  tilting  of  the  interfering  beam  in  conjunction  with  the 
continuous  dark  etching  effect,  which  progresses  normal  to  the  surface, 
produces  the  clear  blazed  profile  shown  in  the  figure.  The  profile  Is 
triangular  or  cusped  depending  on  the  relative  magnitude  of  the  dark  and 
light-enhanced  etching  rate. 

Although  several  different  etching  solutions  (e.g.  NH^OH,  H202,  and  I^O, 
or  KC1,  HC1,  and  H^O)  were  used,  most  of  the  work  was  done  with  the  following 
solution:  H^SO^,  and  H^O.  Since  we  were  primarily  concerned  with  grating 

periods  less  than  400  nm,  it  was  necessary  to  choose  an  etch  composition  with 
a  sufficiently  slow  dark  etching  rate,  so  that  the  products  of  the  etching 
reaction  would  not  Interfere  with  the  development  of  the  grating.  Based  on  an 
Initial  measurement  of  the  dark  and  light-enhanced  etching  rates  for  different 
compositions,  we  selected  the  following  composition  for  our  experiments: 

H2S04:  H202:  H20  -  1:1:50. 

Finally,  in  collaboration  with  a  project  at  Lincoln  Laboratory,  headed  by 
Dr.  Antonio  Sanchez,  we  have  explored  applying  the  same  etching  techniques  to 
CdS.  One  of  the  more  important  findings  has  been  that  photochemical  grating 
formation  is  most  effective  when  etching  a  crystalline  plane  parallel  to  a 
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Fig. 


The  SEM  photograph  of  the  profile  of 
200  run  grating  on  n-type  GaAs. 
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Fig.  3  SEM  photograph  of  500  nm  blared  grating 
on  n-type  GaAs. 
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prismatic  face  of  the  CdS.  In  the  Lincoln  Laboratory  project,  the  grating  is 
necessary  to  produce  distributed  feedback  to  a  blue-green,  electron-beam- 
pumped  laser.  As  a  result,  the  requirement  for  orientation  of  a  specific 
crystal  to  form  the  grating  vill  also  determine  the  crystal  orientation  of  the 
laser . 

b)  Metallic  Gratings 

In  addition,  we  have  used  the  interferometric  technique  described  above 
in  conjunction  with  the  organometalllc  gas-handling  system,  built  as  part  of 
a  separate  JSEP  project,  to  demonstrate  the  direct  deposition  of  metallic 
gratings.  Even  in  the  preliminary  experiments  the  resolution  of  these 
gratings  is  comparable  to  best  values  obtained  by  state-of-the-art  photo¬ 
lithographic  techniques.  Experiments  to  further  explore  the  fabrication  of 
these  gratings  are  currently  being  planned. 

*This  research  was  also  supported  by  the  Air  Force  Office  of  Scientific 
Research/DARPA  under  Grant  F-49620-82-K-0008. 


L.  SPECTROSCOPY  OF  ORGANOMETALL1C  MOLECULES  ON  SURFACES* 

(C.  J.  Chen,  and  R.  M.  Osgood,  Jr.) 

(JSEP  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.M.  Osgood,  Jr.  (212)  280-4462) 

Since  the  development  of  laser  photochemical  direct  writing,  there  has 
been  an  Increasing  interest  in  understanding  the  UV  photodissociation 
mechanism  of  organometallic  molecules.  Uhile  gas-phase  ultraviolet  absorption 
spectra  can  provide  significant  Information  on  the  dissociation  process,  the 
existing  spectroscopic  data  are  not  complete  and  an  analysis  of  the  excited 
states  has  not  been  made.  Although  one  observation  of  the  chemisorbed 
monolayer  of  DMCd  (dimethyl  cadmium)  has  been  reported  the  spectrum  of  the 
physlsorbed  molecular  layer  has  not  been  measured.  The  difficulty  of 
measuring  the  spectrum  of  the  physlsorbed  layer  stems  from  the  fact  that  at 
high  temperatures  the  physlsorbed  layers  always  coexist  with  a  relatively 
dense  ambient  gas.  As  a  consequence,  the  spectrum  of  the  adlayer  is  obscured 
by  that  of  the  gas. 

We  have  recently  begun  a  study  of  the  UV  spectroscopy  of  adsorbed 
metal-alkyl  molecules.  The  spectra,  which  include  both  the  chemisorbed  and 
physlsorbed  layers,  can  be  interpreted  by  comparison  with  our  recent 
measurements  of  gas-phase  metal  alkyl  molecules  as  well  as  by  an  analysis 
based  on  molecular-orbital  considerations. 

Essentially  the  apparatus  consisted  of  a  carefully  designed  UV 
monochromator  and  a  precision  multisurface  cell.  A  nitrogen-purged,  0.3-m 
scanning  monochromator  with  a  1200-groove-per-millimeter  grating  and  a 
deuterium  lamp  provides  a  spectrally-tunable  UV  probe  with  0.6&  resolution. 

The  power  supplies  for  the  anode  and  the  filament  of  the  deuterium  lamp  are 
both  highly  stabilized  to  provide  a  constant  light  Intensity.  A  PDP-11 
microcomputer  controls  the  scanning  of  the  monochromator,  collects  and  stores 
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the  readings  of  the  photomultiplier,  and  then  performs  scheduled  computation. 

The  absorption  cell  is  made  of  316  stainless  steel  with  both  ends 
precisely  ground.  Ten  to  thirty  suprasil  (fused  silica)  discs  can  be  fit  in 
the  cell,  separated  with  C-shaped  spacers  of  25-m  thickness.  The  cell  is 

t 

enclosed  in  a  copper  thermal  enclosure  with  thermal  insulation.  The 
temperature  of  the  cell  is  regulated  to  within  +  0.2°C  by  a  dry-ice  bath  and 
an  electronic  thermal  controller. 

The  gas-handling  system,  with  a  diffusion  pump,  is  made  of  stainless 
steel.  It  was  baked  to  200°  C  in  vacuum  before  use.  The  cell  is  directly 
connected  with  the  gas-handling  system.  The  pressure  is  read  with  a  digital 
capacitance  manometer  to  0.01  torr. 

a)  Gas  Spectra 

The  gas  spectra  of  DMZn,  DMCd  and  DMHg  are  shown  in  Fig.  1.  In  general, 
the  overall  structure  of  each  spectrum  is  similar;  each  consists  of  a  weak 

absorption  knee  and  a  strong  peak  with  almost  equally  spaced  bands.  However, 

f 

for  DMCd,  another  weak  absorption  peak  in  the  short  wavelength  region  is 
apparent.  We  will  not  discuss  the  detailed  assignment  of  the  electronic 
states  here;  instead,  only  our  general  conclusions  will  be  summarized. 

♦ 

Absorption  in  weak  feature  at  the  long  wavelength  end  leads  to  the  immediate 
formation  of  a  metal  atom  in  the  ground  state  and  two  methyl  radicals.  The 
strong  absorption  feature,  whose  cross  section  is  typical  for  an  allowed 
dipole  transition,  produces  a  methyl  radical  and  an  excited  metal  monomethyl 
which  dissociates  after  radiative  emission  to  the  ground  state.  A  third 
feature  is  just  observable  on  the  blue  side  of  the  strong  absorption  band.  In 
this  UV  band,  the  photon  energy  exceeds  the  sum  of  the  dissociation  energy  of 
the  molecule  into  an  excited  metal  atom  and  two  methyl  radicals. 


» 


103 


40,000  45,000  50,000 


Figure  1 


55,000  cm-l 


104 


I 


r 


♦ 


i 


f 


» 


b)  Spectra  of  the  Chemisorbed  Monolayer 

In  the  case  of  DMCd,  a  chemisorbed  monolayer  forms  first  on  the 
substrate.  Since  this  layer  cannot  be  removed  by  pumping,  the  measurement  is 
simple,  requiring  no  subtraction  of  the  gas  spectra  (Fig.  2).  Note  that  there 
is  no  chemisorbed  layer  for  the  related-compound  DKHg.  This  la  consistent 
with  the  fact  that,  unlike  DMCd,  DMHg  is  air  and  moisture  stable.  Ue  believe 
that  the  DMCd  reacts  with  the  OH  radical  on  the  silica  surface  and  produces 
cadmium  monomethyl  with  its  cadmium  atom  firmly  bonded  on  the  surface. 

If  the  monolayer  is  Irradiated  with  an  193  nm  laser  beam,  the  spectnmi 
changes  radically,  thus  indicating  that  the  monolayer  is  activated  or 
dissociated  by  ultraviolet  light.  Exposure  to  air  also  changes  the  spectrum 
strongly.  Note  that  at  257  nm,  the  absorption  strength  of  the  chemisorbed 
layer  is  about  7.8  times  greater  than  for  an  equivalent  quantity  of  molecules 
in  the  gas  phase. 

c)  Spectra  of  the  Fhyslsorbed  Layer 

The  spectra  of  the  physlsorbed  DMCd  layers  at  0°  C  is  shown  in  Fig.  3. 
Twenty  surfaces  were  used  to  obtain  this  spectra  in  a  cell,  with  a  free  path 
of  0.40  nm.  Under  these  conditions  the  absorbtion  due  to  the  gas  is  only 
about  30Z  of  the  total  absorption  signal.  Thus,  even  the  unprocessed  cell 
spectra  are  radically  different  from  the  gas  spectra.  The  gas  portion  of  the 
spectra  is  then  subtracted  accurately  by  an  on-line  computer.  Several 
features  of  the  absorption  spectra  of  the  physlsorbed  layers  are  notable. 
First,  the  strong  absorbtion  peak  is  slightly  redshifted  and  its  vibrational 
structure  is  almost  completely  washed  out.  Note,  however,  that  there  still 
appears  to  be  some  residual  vibronlc  structure.  Second,  there  is  a  strong 
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pressure  and  temperature  dependence  to  the  adsorption  depth  of  the  adlayer. 
Third,  there  Is  a  long  wavelength  tall  extending  directly  to  visible 
wavelengths.  The  wavelength  dependence  of  this  tail  is  very  flat  and  the 
pressure  dependence  of  it  is  radically  different  from  that  of  the  main  peak. 

At  low  pressures  the  dependence  of  this  tail  is  nonlinear;  in  contrast,  the 
linear  dependence  of  the  absorption  in  the  main  peak  is  consistent  with 
measurements  of  the  DMCd  Isotherm. 

The  interpretation  of  the  physisorbed  layer  spectra  is  as  follows.  Most 
of  the  molecules  in  the  physisorbed  layer  are  in  a  liquid-like  environment,  in 
which  the  dense  molecular  environment  can  broaden  spectral  features  due  to 
shift  in  the  intramolecular  potential  of  the  upper  and  lower  energy  states. 

The  long-wavelength  tall  is  due  to  the  Rayleigh  scattering  by  the  nonunl- 
formity  thickness  of  the  adsorbed  layer.  This  nonuniformity  increases  near 
the  condensation  point.  Note  that  the  high  surface  density  and  the  slight 
spectral  redshlft  of  the  adsorbed- layers  may  play  an  important  role  in 
photochemlcally  producing  the  first  metal-a'-c.j  clusters  on  the  surface. 

*Thl8  research  was  supported  by  the  Army  Research  Office  under  Grant 
DAAG29-82-K-0009 . 
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III.  CHARACTERISTICS  AND  PROPERTIES  OF  ENERGETIC  ION  AND  RADICAL  BEAMS 

A.  THE  CAPILLARITRON :  A  LASERLESS  MEANS  OF  PRODUCING  COLD  RADICALS* 

(B.  Brady,  T.  Allik,  G.  W.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  V.  Flynn  (212)  280-4162) 

Ions  and  radicals  are  important  In  varied  chemical  applicationa.  In  the 
field  of  microelectronics,  ions  and  radicals  of  several  kinds  are  Important  in 
the  etching  and  implantation  of  semiconductor  surfaces.  Ions  and  radicals  are 
also  important  in  many  gas  phase  reactions  and  their  spectra  are  often  poorly 
understood.  Since  ions  and  radicals  are  generally  created  with  large  amounts 
of  excess  energy,  there  are  usually  a  large  number  of  different  states 
populated.  This  large  variety  of  states  gives  rise  to  congested  spectra  which 
are  difficult  to  interpret  and  also  provides  for  a  myriad  of  competing 
channels  in  the  reactions  of  these  substances.  If  one  were  able  to  create 
ions  and  radicals  in  Just  a  few  of  their  low  energy  states,  it  would  simplify 
the  spectroscopy  and  allow  for  the  preparation  of  specific  states  with 
subsequent  study  of  state-to-state  reaction  rates.  It  was  with  this  goal  in 
mind  that  we  first  began  to  look  at  the  capillarltron. 

The  capillarltron  consists  of  a  capillary  nozzle  held  at  high  positive 

voltage  and  an  annular  ground  plate  a  few  millimeters  in  front  of  the  nozzle. 

In  operation,  a  gas  is  fed  into  the  nozzle  at  a  moderate  to  high  pressure  and 

allowed  to  expand  into  a  chamber  at  a  significantly  lower  pressure.  When 

voltage  is  applied  to  the  nozzle,  a  plasma  forms  between  the  nozzle  tip  and 

the  ground  plate  and  the  gas  is  ionized  and  sometimes  fragmented.  In  our 

application,  the  typical  backing  pressure  is  one  to  five  atmospheres  with  a 

-3  -4 

source  chamber  pressure  of  10  to  10  torr.  The  voltage  applied  is  usually 
between  5 00  and  1000  volts  and  the  ion  current  at  the  nozzle  is  a  few 
mllllamps. 


The  experimental  setup  has  two  additional  differentially  pumped  chambers 
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downstream  of  the  source  chamber.  After  a  pathlength  of  about  a  centimeter  in 
the  first  chamber,  the  beam  from  the  nozzle  is  skimmed  and  collimated.  At  the 
end  of  the  third  chamber  is  a  quadrupole  mass  spectrometer.  Pressures  in  the 
second  and  third  chambers  are  generally  kept  at  about  10  2  torr. 

When  no  voltage  is  applied  to  the  capillaritron,  the  nozzle  and  chamber 
arrangement  function  as  a  supersonic  molecular  beam  apparatus.  As  a  gas 
expands  into  the  chamber,  it  is  cooled  because  the  molecules  in  the  beam 
collide  and  excess  energy  is  removed  by  the  molecules  which  travel  out  of  the 
beam.  This  process  occurs  until  all  the  molecules  on  the  beam  axis  are  moving 
at  roughly  the  same  speed  at  which  point  the  beam  becomes  collisionless.  The 
velocity  distribution  along  the  axis  at  this  point  can  be  described  by  a 
Boltzmann  equation: 

f<v>  a  (~  )  exP  t  <v-vs> 2/°ts  1 
'“s  ' 

where  vg  is  the  stream  velocity  and  is  significantly  higher  than  those  in  bulk 
samples  at  room  temperature.  ag  is  the  width  parameter  and  is  used  to 
determine  the  effective  translational  temperature  by  the  equation: 

as  -  (2kTs/m)1/2 

This  effective  translational  temperature,  Tg,  is  very  low  when  the  beam 
becomes  collisionless.  Other  degrees  of  freedom  of  the  expanding  gas  are  also 
cooled  in  this  process,  but  only  to  the  extent  that  they  couple  to  the 
translational  energy  distribution.  The  ions  and  radicals  from  the  capillari¬ 
tron  are  created  early  in  this  expansion,  and  the  primary  goal  of  our  charac¬ 
terization  of  the  capillaritron  is  to  determine  if  these  species  are  in¬ 
fluenced  enough  by  the  supersonic  expansion  under  normal  operating  conditions 

no 
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Co  cool  significantly.  Our  efforts  to  determine  the  energy  distribution  of 
the  capillaritron  output  can  be  divided  into  two  parts:  spectroscopic  measure¬ 
ments  and  time-of-f light  (TOF)  measurements. 

The  spectroscopic  measurements  made  on  the  capillaritron  to  date  Include 
observation  of  visible  emission  from  atoms,  molecules  and  ions  excited  in  the 
plasma,  and  measurement  of  laser  Induced  fluorescence  from  these  species 
further  downstream.  The  emission  at  the  nozzle  tip  was  initially  observed  with 
an  optical  multichannel  analyzer  (OMA).  This  device  consists  of  an  array  of 
five  hundred  detectors  directly  coupled  to  the  output  end  of  a  monochromater . 

A  wide  range  of  wavelengths  can  be  surveyed  with  the  OMA  at  the  same  time  but 
since  there  is  no  exit  slit  on  the  monochromator,  the  ultimate  resolution  is 
fixed  by  the  width  of  the  detector.  Vhen  a  rare  gas  such  as  Ar  or  He  is  used 
in  the  capillaritron,  we  see  only  atomic  emission  lines. ^  The  emission  from 
ions  is  expected  to  be  much  weaker;  none  was  observed.  When  carbon  tetra- 
fluorlde  (CF^)  is  used  as  the  fuel  gas,  we  obtain  the  spectrum  shown  as  the 
lower  trace  in  Fig.  1.  Sharp  lines  can  be  seen  from  both  carbon  ion  (C  ) 
and  neutral  fluorine  and  diffuse  bands  are  observed  centered  at  roughly  570  nm, 
500  nm,  300  nm,  and  230  nm. 

The  first  task  in  analyzing  these  spectra  is  to  determine  the  emitting 
species,  after  which  information  may  be  obtained  about  their  energy  distribu¬ 
tions.  The  bands  at  300  nm  and  570  nm  have  been  seen  in  CF^  discharges  by 
other  workers.  The  emission  at  570  nm  has  also  been  detected  recently 

in  the  photodlssoclatlon  of  CF^Br  and  CF^Cl  and  assigned  to  the  CF^ 
radical. Unfortunately,  since  this  band  is  a  continuous  emission, 
thought  to  end  on  a  dissociative  state,  it  indicates  nothing  about  the  energy 
distribution;  it  does,  however,  confirm  the  presence  of  CF^  in  the  discharge. 
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The  ultraviolet  emission  hands  could  be  from  either  OT  CP, . 

investigators  have  ruled  out  CF^  as  an  emitting  species  in  an  RF  discharge 

(3) 

on  the  basis  of  lifetime  Information;  however,  we  were  not  able  to 
measure  lifetimes.  The  CF^  radical  also  shows  sharp  band  structure  when  it 
emits, the  absence  of  this  structure  either  indicates  that  the 
emitting  species  has  more  atoms  or  that  it  is  very  hot.  When  CF^  was  seeded 
in  helium  and  the  emission  observed  with  a  photomultiplier  tube  and  a 
monochromator  with  adjustable  slits,  some  band  structure  was  observed  as  shown 
in  the  upper  trace  in  Fig.  1.  The  shape  of  the  band  and  spacing  of  the  peaks 
is  similar  to  that  found  for  CF^  in  Ref.  5  and  6,  but  the  band  is  centered  at 
slightly  longer  wavelengths.  The  position  of  this  band  could  be  explained  if 
the  emission  were  from  higher  vibrational  states  than  in  Ref.  5,  but  the 
assignment  remains  uncertain. 

The  emission  band  centered  at  400  nm  has  not  been  observed  in  other 
discharges.  The  small  peaks  on  this  band  do  seem  to  correspond  well  to  a 

series  of  lines  seen  in  the  source  region  by  Van  Sprang,  et  al.  The  average 

-1  -1 
spacing  in  our  case  is  680  cm  and  in  their  spectrum  it  was  718  cm  .  This 

was  said  to  correspond  well  to  a  vibrational  progression  in  CF^+  at  730  cm  * 

and  assigned  accordingly.  This  assignment  is  in  doubt,  however,  because  other 

+  (9) 

investigators  have  found  that  CF^  decays  non-radiatively. 

One  reaso.i  why  the  emission  spectra  do  not  show  much  cooling  is  because 
they  sample  molecules  close  to  the  nozzle  and  in  electronically  excited 
states.  One  technique  which  allows  us  to  probe  molecules  in  their  ground 
state  at  variable  locations  with  very  good  spatial  resolution  and  sensitivity 
is  laser-induced  fluorescence  (L1F).  Another  advantage  of  LIF  is  that  it 
allows  us  to  probe  one  type  of  molecule  at  a  time  instead  of  collecting 
fluorescence  from  every  species  that  might  be  present  (using  the  OMA).  We 


used  Che  KrF  line  of  an  excimer  laser  to  excite  CFg  radicals  in  their  ground 
vibrational  and  electronic  states  and  resolved  the  resulting  fluorescence. 

The  spectrum  obtained  is  shown  in  Fig.  2  and  is  characteristic  of 

Since  only  one  state  was  probed,  we  do  not  know  what  the  energy 
distribution  is  but  this  does  confirm  the  presence  of  CF^  in  the  beam  and 
provides  a  basis  for  future  experiments  which  will  determine  the  internal 
energy  distribution. 

The  second  type  of  measurement  done  on  the  beam  was  a  TOF.  A  chopper 
wheel  with  a  small  slit  was  placed  in  chamber  2.  When  the  wheel  was  spun  a 
segment  of  the  beam  was  chopped  out  and  its  flight  from  the  chopper  whee’  to 
the  mass  spectrometer  was  timed.  In  this  manner,  the  velocity  distribution  of 
the  various  species  in  the  beam  was  determined.  When  measuring  the  TOF  of  the 
ions,  the  ionizer  on  the  mass  spectrometer  could  be  turned  off  and  the  mass 
spectrometer  used  to  separate  the  various  fragment  ions,  CF^+,  CF^*,  CF+,  and 
C  to  determine  their  individual  time  of  flights.  For  the  neutral  time  of 
flights,  the  ionizer  had  to  be  turned  on,  and  ions  from  the  capillarltron 
deflected.  Unfortunately,  CF^  fragments  substantially  in  the  ionizer  and 
fragments  created  there  could  not  be  distinguished  from  those  created  in  the 
capillarltron  discharge. 

When  At  is  put  in  the  capillarltron  the  energy  of  the  Ar+  ions  produced 

ranges  from  0  to  1000  eV.  CF^  produces  CF^+,  CF2*,  CF*  and  C+,  as  mentioned 

before,  and  each  of  these  exhibit  a  narrower  distribution  than  the  Ar.  Part 

of  the  reason  for  the  narrower  distributions  from  CF,  is  that  it  breaks  down 

4 

at  half  the  voltage  required  for  Ar.  The  CF*  ions,  however,  seem  to  exhibit  a 
very  narrow  distribution.  We  cannot  presently  determine  exactly  how  narrow 
this  distribution  is  because  the  uncertainties  in  the  TOF  measurements  become 
very  significant  when  dealing  with  these  fast-moving  particles.  More  accurate 
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measurement  of  these  velocities  is  required  before  anything  conclusive  can  be 
stated  about  the  distributions. 

For  neutrals,  the  TOF  measurements  are  extremely  accurate  because  the 
particles  are  traveling  at  about  1/1000  the  speed  of  the  ions.  For  a  neutral 
Ar  beam  vith  no  voltage  on  the  capillarltron,  the  translational  temperature  is 
2.8K.  For  CF^,  Tg  Increases  to  17K  because  the  Internal  degrees  of  freedom  of 
the  polyatomic  do  not  cool  quite  as  well.  When  a  discharge  is  put  on  the 
capillarltron,  the  velocity  distribution  of  the  neutral  component  of  the  beam 
does  not  change  significantly.  As  can  be  seen  in  Figs.  3  and  4,  the  average 
velocity  of  the  beam  speeds  up  slightly  and  for  Ar,  the  translational  tempera¬ 
ture  increases  by  2K.  For  CF^  the  temperature  does  not  appear  to  change  at 
all.  We  saw  no  warming  upon  turning  on  the  discharge.  At  first  this  seems 
surprising  since  the  temperature  of  the  electrons  in  the  capillarltron  plasma 
is  several  thousand  °K  and  the  ions  are  hot  as  observed  previously.  The 
electrons  are  light,  however,  and  cannot  transfer  much  momentum  to  the  other 
species;  the  ions  are  apparently  a  small  enough  portion  of  the  beam,  1-2Z, 
that  they  do  not  dominate  it. 

The  importance  of  this  discovery  is  twofold.  First,  since  radicals  are 
neutrals,  if  their  velocity  distribution  were  significantly  different  from  the 
other  neutrals,  we  would  have  observed  them  as  an  additional  peak  or  as  a 
broadening  of  the  main  peak.  Secondly,  this  means  the  bath  in  which  the 
radicals  are  expanding  remains  cold  even  if  the  radicals  themselves  are  formed 
hot.  The  OH  radical'  '  requires  500  to  1000  collisions  for  the  rotational 
temperature  to  equilibrate  to  the  temperature  of  the  buffer  gas.  In  our  beam, 
we  have  more  collisions  than  that  between  the  plasma  and  the  collisionless 
regime  and  the  buffer  gas  is  at  5K.  This  means  we  should  be  able  to  cool 
radicals,  at  least  rotationally,  after  they  have  been  formed  in  the  plasma. 
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TRANSLATIONAL  TEMPERATURE  (K) 


Future  experiments  on  this  project  will  involve  the  determination  of 
rotational  and  vibrational  temperatures  of  the  radicals  in  the  beam  with  laser 
induced  fluorescence  (LIF),  measuring  the  TOF  of  the  radicals  directly  and 
determining  the  energy  of  the  ions  with  an  electrostatic  energy  analyzer. 
Scanning  a  dye  laser  brings  various  transitions  into  resonance  which  will  take 
molecules  from  various  levels  in  the  ground  state  to  some  particular  excited 
level.  By  then  monitoring  the  total  fluorescence  from  the  excited  state  and 
knowing  the  Franck-Condon  factors  for  the  laser  transition,  one  can  determine 
the  population  in  the  various  rotational~vlbrational  levels  of  the  ground 
state.  This  population  distribution  can  then  be  analyzed,  assuming  it  is 
Boltzmann;  to  get  a  rotational  and  vibrational  temperature.  The  CN  radical, 
produced  by  a  discharge  of  acetonitrile  (CH^CN)  in  Ar,  has  a  simple  LIF 
spectrum  which  can  be  probed  by  the  fundamental  of  the  BBQ  dye  at  387nm.  This 
experiment  will  be  done  in  the  near  future. 

Becuase  the  ionization  potential  of  CF^  is  about  8  eV  lower  than  that  of 
CF^,  it  is  possible  to  ionize  only  CF^  by  turning  down  the  electron  energy  on 
the  ionizer  of  the  mass  spectrometer.  This  would  eliminate  fragmentation  of 
the  parent  and  allow  direct  observation  of  CF^  created  by  the  capillar itron. 
Another  strategy  would  be  to  modulate  the  current  on  the  caplllaritron  and 
look  for  a  modulated  signal  in  the  time  of  flight.  This  method  seems  more 
difficult  but  either  method  should  allow  direct  determination  of  the  velocity 
distribution  of  the  radicals  from  the  caplllaritron.  Determination  of  the 
velocity  distribution  of  the  ions  is  simple  and  involves  merely  observing 
their  deflection  in  an  electric  field. 

In  conclusion,  the  creation  of  cold  radicals  and  possibly  cold  ions  with 
the  caplllaritron  appears  to  be  possible.  Further  measurements  are  needed  to 
confirm  this  conclusion  and  these  measurements  should  be  done  in  the  near 
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future.  Then  if  the  capillarltron  does  prove  to  be  a  source  of  cold  radicals 
and  Iona,  It  should  prove  to  be  an  important  tool  in  deciphering  spectra  and 
kinetic  schemes. 

*This  research  vas  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78ER04940. 
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IV.  ENERGY  TRANSFER  PROCESSES  IN  QUANTUM  ELECTRONICS 

I  * 

A.  PHOTOCHEMICAL  AND  PHOTOPHYSICAL  DYNAMICS  OF  I2  IN  A  RARE  GAS  MATRIX* 

(P.  B.  Beeken,  E.  A.  Hanson,  G.  V.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  U.  Flynn  (212)  280-4162) 

I 

The  iodine  atom  laser  is  rapidly  gaining  importance  for  use  in  projects 
such  as  Inertial  confinement  fusion.  One  potential  drawback  to  obtaining  high 

*  power  output  from  this  laser  is  the  potential  for  Interference  from  the  low 
lying  electronic  states  of  the  byproduct,  I2«  Previous  work  on  Br2  Isolated 
in  rare  gas  matrices  has  helped  characterize  the  relaxation  times  and  energy 
transfer  pathways  among  the  low  lying  electronic  states  of  diatomic  halogens. 

3 

Specifically,  the  emission  from  the  lowest  state.  A'  tt^,  was  viewed 

3 

directly.  This  state  and  the  one  slightly  higher  in  energy,  A  ir2u,  in 
I2  could  prove  to  be  efficient  quenchers  of  excited  iodine  atoms.  A  complete 
analysis  of  the  relaxation  pathways  for  I2  could  be  of  great  Importance  for 
optimising  the  I  laser. 

f  The  following  abstract  is  from  a  paper  to  appear  in  the  Journal  of 

Chemical  Physics  in  May  of  1983  describing  our  work  on  energy  transfer  in 
matrix  Isolated  iodine  molecules. 

*  The  photochemical  and  photophysical  properties  of  I2  in  Ar,  Kr,  and  Xe 

matrices  have  been  Investigated  and  are  reported  here.  Following  laser 

3  +  1 

excitation  high  into  the  B  it(Ou)  and  dissociative  state,  prompt  emission 

3+3  3 

I  is  observed  from  the  B  w(0u),  A  ir^,  and  A'  *2u  electronic  states  depending 

on  the  matrix  host.  In  all  matrices  the  emission  attributed  to  the  B  elec¬ 
tronic  state  lies  200  nm  to  the  red  of  the  excitation  wavelength  and  is 

*  generally  featureless.  The  time  resolved  emission  for  this  state  has  a 

3 

detector  limited  rise  and  fall.  The  emission  attributed  to  the  A  r,  elec- 

lu 


> 
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tronic  state  was  resolved  with  a  monochromator  and  the  source  determined  to  be 


the  v'  •  0  level  of  the  A  state.  The  lifetime  for  this  emission  is  260  usee 
and  was  found  to  be  roughly  the  same  in  all  matrices.  Emission  from  the  A' 
*tate  was  also  detected  but  only  in  Xe  matrices.  The  measured  lifetime 
for  this  state  was  6.3  msec.  Data  from  a  number  of  previous  experiments  on 
matrix  isolated  halogens  in  Ar,  Kr,  and  Xe  are  reviewed  and  a  consistent  model 
for  the  relaxation  pathways  for  lntra-molecular  electronic  energy  transfer  is 
presented.  The  implications  of  this  model  for  studies  of  the  solvent  cage 
effect  on  1^  photodissociation  and  the  development  of  the  I*  chemical  laser 
are  also  discussed. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE-80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78E*04«  1. 
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B.  STUDIES  OP  IODINE  ISOLATED  IN  LOW  TEMPERATURE  MATRICES  OF  CARBON 
TETRACHLORIDE  * 

# 

(K.  Kasturi,  P.  Beeken,  G.  W.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

9  In  the  past  few  months  we  have  been  Involved  in  studies  aimed  at 

exploring  the  effect  of  the  matrix  cage  on  the  photodissociation  of  iodine  in 
carbon  tetrachloride  at  low  temperatures.  Iodine  is  excited  by  a  laser  into  a 

*  superposition  of  the  bound  B  state  and  the  repulsive  continuum  (Fig.  1). 

As  a  result  of  confinement  by  the  frozen  cage,  most  of  the  iodine  atoms  formed 
upon  dissociation  are  prevented  from  escaping  and  recombination  is  facilitated. 

The  effect  of  the  cage  on  the  photodlssoclatlon  of  iodine  in  liquid 
CCl^*’)  and  Xenon^^  has  been  studied  by  picosecond  techniques.  The  rate  of 
recombination  was  studied  by  probing  the  recovery  of  the  ground  state 

*  population  of  iodine  molecules  and  two  rates  were  reported: 

1)  a  fast  rate  (timescale:  100  picoseconds)  due  to  the  geminate  recombination 

of  iodine  atoms  derived  from  the  same  iodine  molecule.  2)  a  slow  rate 
“8 

t  (timescale  >10  aec)  attributed  to  the  recombination  of  iodine  atoms  that  have 

escaped  from  their  cages.  Theoretical  calculations  indicate  that  only  low 
lying  vibrational  levels  are  probed  in  these  experiments,  and  further  suggest 

*  that  alow  relaxation  down  the  vibrational  manifold  of  the  ground  state  makea 

(3) 

an  Important  contribution  to  the  recombination  dynamics. 

Studies  of  iodine  in  low  temperature  matrices  performed  in  our  labora- 
t  (4) 

T  tory  have  shorn  that  upon  photodissociation  and  subsequent  recombination 

the  excited  state  population  undergoes  relaxation  and  is  partitioned  between 
the  different  electronic  states  B,  A  and  A' .  (Relaxation  to  the  ground  state 

*  could  not  be  monitored  in  these  experiments.)  The  states  accessed  depend  upon 
the  nature  of  the  matrix,  the  A'  state  being  observed  only  in  Xe  matrices. 
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(  A\ 

It  has  been  suggested'  that  the  low  lying  A*  state  can  be  populated  by 
dissociation  along  the  state  which  correlates  at  long  lnternuclear 
distances,  attainable  only  In  the  large  Xe  cage,  to  the  A'  state.  The  higher 
energy  A  state  Is  accessible  even  In  small  cages.  CCl^  Is  expected  to  fora 
large  cages,  therefore  one  aay  speculate  that  excited  Iodine  can  cross  over  In 
a  CCl^  aatrix  to  the  A'  state,  which  has  a  long  llfetlae.  Consequently,  a 
possible  explanation  for  the  alow  recovery  of  the  ground  state  population  of 
Ig  In  CCl^  could  be  slow  relaxation  from  the  A'  state. 

Initially,  our  goal  was  to  study  the  channeling  of  energy  Into  the 
different  electronic  states,  and  to  Investigate  the  variation  of  the 
relaxation  dynamics  with  temperature.  To  achieve  this  end,  the  emission  of 
I 2  from  CCl^  matrices  was  subjected  to  spectral  and  temporal  analysis. 

Spectral  resolution  of  the  emission  has  yielded  intriguing  results  which  are 
described  below. 


Summary  of  results  obtained  with  I^  1°  2®  matrices. 

For  purposes  of  comparison  it  is  useful  to  summarize  the  Xe  matrix 
14) 

results.  Analysis  of  the  IR  fluorescence  arising  from  Ij  in  Xe  matrices 
reveals  two  signals: 

l)  a  fast  decay  with  a  lifetime  of  240  usee  +  20  usee,  which  originates  from 
the  1-2  u  spectral  region.  When  resolved  the  spectrum  displayed  a  progression 
of  16  peaks  due  to  transitions  from  the  v-0  level  of  the  A  state  to  various 
vibrational  levels  of  the  X  state;  2)  a  slow  decay  in  the  2-3  u  spectral 
region  with  lifetime  of  6  msec  +  1  msec.  The  signal  was  too  weak  to  be 
resolved  with  a  monochromator  but  was  presumed  to  originate  from  the  A'  state. 
A  near  red  signal  was  observed  with  a  very  fast  decay  (15-40  nsec)  and  a 
featureless  spectral  profile  peaking  at  740  nm. 
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Experimental 

Matrices  were  formed  by  depositing  mixtures  of  Iodine  and  carbon 
tetrachloride  on  a  cesium  iodine  window  held  at  12°R  by  an  Air  Products  202A 
helium  refrigerator.  The  composition  of  the  samples  of  iodine  in  carbon 
tetrachloride  was  varied  between  1:1000  and  1:2000.  The  matrices  were 
irradiated  with  the  output  of  an  excimer- pumped  dye  laser  (Lambda  Physik  EMC 
101,  FL2002)  in  the  wavelength  range  500-540  nm.  The  fluorescence  obtained 
from  the  samples  was  resolved  spectrally  as  well  as  temporally.  The  IR 
signals  from  1.0  to  6.0  p  were  detected  using  a  cooled  InSb  (77°K)  detector 
(risetime:  1  psec).  Signals  in  the  near-red  up  to  900  nm  were  monitored  with 
an  RCA  C31034A  GaAs  photomultiplier  tube  (rlsetlme:  15  ns).  The  signals  were 
resolved  with  either  narrow  band  filters  or  an  Instruments  SA  320 
monochromator . 

The  time  resolved  signals  were  amplified  with  a  Tektronix  AM502  amplifier 
(bandwidth:  1  Hz-lMhz)  and  signal  averaged  with  a  Biomation  8100-Nlcolet  1170 
combination.  Spectral  scans  were  taken  with  a  PAR  CW-1  boxcar  Integrator  and 
recorded  on  a  Houston  Instruments  strip  chart  recorder. 

Liquid  phase  absorption  spectra  were  taken  with  a  Cary  Model  17  Spectro¬ 
photometer  in  the  wavelength  range  900-1800  nm,  using  quartz  cuvettes  of 
pathlength  1  cm. 

Results  and  Discussion 

Upon  irradiating  matrices  of  pure  CCl^,  no  emission  was  observed. 
Irradiation  of  matrices  of  iodine  in  CC1&  yielded  emission  in  the  near  red  and 
in  the  Infrared.  Emission  in  the  near  red  was  monitored  with  the  photo¬ 
multiplier  tube  in  the  region  700  to  900  na,  but  no  spectrally  resolved  lines 
could  be  obtained  upon  observation  through  a  monochromator.  The  signal 
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decayed  very  fast  but  the  lifetiae  could  not  be  measured.  This  fluorescence 
probably  arises  from  the  B  state  since  It  corresponds  in  wavelength  range  and 
fast  decay  to  that  obtained  froa  iodine  in  rare  gas  aa trices. 

The  total  signal  in  the  1. 0-6.0  u  region  was  characterized  by  three 
lifetiaes:  1)  a  very  fast  decay  with  a  lifetiae  of  25  Msec  which  was  probably 
due  to  formation  of  an  ealssion  froa  an  (^2  <^mer*  38  suggested  by  studies 
of  the  temperature  dependence  of  the  emission.  2)  a  fast  decay  with  a 
lifetiae  of  148  Msec.  3)  a  slow  decay  with  a  lifetime  of  227  usee.  The 
rlsetiae  of  the  emission  was  detector  limited  in  all  cases. 

The  emission  was  resolved  with  various  narrow  band  filters,  and  the 
tiae-reaolved  signals  analyzed.  From  the  fitted  curves,  the  relative 
intensities  of  the  different  decays  were  determined  and  plotted  as  a  function 
of  the  bandpass  of  the  filter  (see  Fig.  2).  The  nature  of  the  signals  did  not 
change  with  the  concentration  of  iodine  in  the  matrix,  and  the  emission 
intensity  scaled  roughly  linearly  with  laser  power. 

The  emission  in  the  1-2  m  region  was  resolved  with  a  monochromator.  The 
ealssion  spectrua  so  obtained  (Fig.  3)  exhibited  two  peaks  centered  at 
approximately  1.30  and  1.34  m  against  a  featureless  background  which  extended 
beyond  2  m> 

The  ealssion  at  various  wavelengths,  resolved  with  the  aonochroaator, 
corresponding  to  different  regions  of  the  spectrua,  was  subjected  to  teaporal 
analysis.  In  the  vicinity  of  the  1.30  and  1.34  u  peaks  the  ratio  of  the 
intensity  of  the  slow  to  the  fast  lifetiaes  decreased  dramatically  (see 
Fig.  4),  as  coapared  to  other  regions  of  the  spectrum,  thereby  implicating  the 
fast  decay  as  the  source  of  the  spectral  peaks.  This  was  confirmed  when  the 
saas  spectrua  was  obtained  with  the  bandwidth  of  the  aapllfer  decreased  (to 


lKhz-lOOkhz),  thereby  filtering  out  the  slow  lifetiae.  A  scan  was  performed 
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2  Ratios  of  the  intensities  of  the  slow/fast  and  very  fast/fast 
decays  plotted  against  the  bandpass  of  the  filter. 


Fig.  3  Emission  spectrum  of  ^  in  a  CCl^  matrix  held  at  12  K. 
The  composition  of  the  matrix  was  1:1000  I„:CCl,. 


The  ratios  of  the  intensities  of  the 
a)  Slow/fast  and  b)  very  fast/fast  decays  measured  at 
various  wavelengths.  The  fluorescence  was  resolved  with 
a  monochromator,  signal  averaged  and  analyzed.  The  ratios 
of  the  fitted  intensities  were  then  determined. 


with  the  delay  on  the  boxcar  integrator  increased  to  500  paec,  thereby 
integrating  the  intensity  of  the  eaisslon  in  the  tail  of  the  fluorescence 
signal,  and  discriminating  against  the  fast  decay.  The  profile  of  the 
spectrum  changed,  with  the  intensity  of  the  peaks  decreasing  relative  to  the 
background.  This  further  suggests  that  the  peaks  in  the  spectrum  are  due  to 
the  fast  emission  signal. 

Comparing  these  results  with  those  obtained  for  iodine  in  rare  gas 

matrices,  the  fast  decay  is  seen  to  resemble  the  emission  from  the  A  state  of 

iodine  in  both  wavelength  range  and  lifetime.  (Assuming  purely  radiative 

decays,  the  extrapolated  gas  phase  lifetime  of  the  slow  decay  is  315  Msec 

while  that  of  the  A  state  emission  in  the  rare  gas  matrices  varies  between 

426  and  514  Msec;  the  calculated  radiative  lifetime  in  the  gas  phases  is 

220  usec.^^  In  the  case  of  rare  gas  matrices  the  differences  between  the 

extrapolated  and  calculated  lifetimes  can  be  attributed  to  changes  in  the 

(4) 

Franck-Condon  factors  induced  by  the  matrix.  However,  no  conclusions  can 
yet  be  drawn  as  to  the  identity  of  the  state  giving  rise  to  the  fast  decay, 
since  the  spectrum  obtained  does  not  correspond  to  emission  from  the  A  state 
of  iodine,  as  obtained  in  rare  gas  matrices. 

An  absorption  spectrum  of  a  concentrated  solution  of  iodine  in  liquid 
CCl^  was  taken,  to  see  whether  the  solution  spectrum  displayed  features 
similar  to  that  of  the  emission  spectrum  of  the  matrix  sample.  Scanning  in 
the  .9-1.8  m  range  no  spectral  features  could  be  discerned,  except  for  a  very 
weak  absorption  at  1.14  u.  Further  scans  are  planned  with  cuvettes  of  longer 
pathlength  (up  to  10  cm.).  If  any  absorption  is  discovered  in  this  region 
a  study  of  the  concentration  dependence  of  the  absorbance  would  help  to 
elucidate  the  nature  of  the  species  responsible  for  the  emission.  If  the  same 
transition  is  responsible  for  the  absorption  as  well  as  the  emission  spectra. 


solution  phase  studies  could  help  to  cast  light  upon  the  nature  of  the 
emitting  species  in  the  matrix  sample.  However,  a  survey  of  the  literature 
has  not,  as  yet,  revealed  the  existence  of  any  such  absorption  in  solutions  of 
iodine  in  CCl^. 

Temperature  dependent  studies 

A  temporal  analysis  of  the  total  emission  intensity  in  the  1-6  p  region 
was  performed  as  a  function  of  temperature.  The  temperature  of  the  matrix  was 
raised  in  increments  of  10°K  to  160°K;  after  annealing  the  temperature  was 
reduced  to  12°K.  The  matrix  was  once  more  heated  up  to  160°K  in  the  case  of 
1: 1000: I^tCCl^  matrices. 

At  each  temperature,  the  ratios  of  the  initial  Intensities  of  the 
slow/f ast  (Fig.  5)  and  the  very  fast/fast  emissions  (Fig.  6)  were  determined. 
It  was  found  that  the  ratio  of  the  slow/fast  lifetimes  increased  to  a  maximum 
at  60°K,  stayed  roughly  constant  till  140°K,  and  then  began  to  drop  sharply. 
Similar  behavior  was  observed  for  the  ratio  of  the  very  fast/fast  lifetimes. 
Thus,  with  an  Increase  in  temperature,  the  intensities  of  the  slow  and  fast 
lifetimes  decrease  relative  to  that  of  the  very  fast  lifetime.  In  the  case  of 
1:1000  matrices  it  was  seen  that  after  annealing  and  cooling  down  the  matrix, 
the  relative  intensities  of  the  decays  at  12°K  were  nearly  the  same  as  those 
seen  at  160°K.  For  matrices  of  other  composition,  the  relative  intensities  at 
12°K  before  and  after  annealing  are  not  the  same,  but  the  very  fast  lifetime 
intensity  is  still  large  compared  to  each  of  the  others.  Upon  a  further 
Increase  in  temperature  the  relative  intensity  of  the  very  fast  lifetime 
continues  to  Increase  at  the  expense  of  the  others.  Thus  annealing  causes  an 
Irreversible  change  in  the  emitting  species,  suggesting  that  the  iodine  dimer 
might  be  responsible  for  the  very  fast  decay.  As  the  temperature  increases. 


132 


Fig.  6  Ratio  of  the  intensities  of  the  very  fast/fast  decays  plotted 
against  temperature  a)  before  annealing  b)  after  annealing 
a  matrix  of  I2  in  CCl^  (Composition  of  matrix:  1:1000  1„:CC1.) 


the  matrix  softens ,  permitting  the  diffusion  of  Iodine  molecules  which  can 
clump  together  to  form  dimers*  Upon  cooling  the  matrix  after  annealing,  this 
state  of  the  matrix  is  frozen  out.  As  the  temperature  increases,  fluctuations 
in  the  cage  size  may  be  expected  to  become  larger  due  to  increased  molecular 
motion.  As  a  result,  the  A*  state  is  expected  to  become  more  accessible  with 
an  increase  in  temperature.  Since  the  relative  intensity  of  the  slow  decay 
increases  with  temperature  this  emission  might  arise  from  the  A'  state  as  in 
the  case  of  iodine  in  rare  gas  matrices. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78ER04940. 
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C.  ENERGY  TRANSFER  AMONG  THE  SYMMETRIC  STRETCH  AND  THE  BENDING  MODES  OF  CO.: 

AN  INFRARED  DOUBLE  RESONANCE  STUDY*  * 

(C. F.  {food,  J.O.  Chu,  and  G.W.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  V.  Flynn  (212)  280-4162) 

Vibrational  energy  transfer  in  CO^  has  been  extensively  studied  and 
accurate  rates  for  both  intermolecular  and  intramolecular  processes  have  been 
measured  as  well  as  the  temperature  dependence  of  these  rates. ^  jn 

spite  of  the  vast  amount  of  information  obtained  on  energy  transfer  rates  in 
COj,  there  is  an  Important  question  that  remains.  The  fundamental  of  the 
bending  mode  is  known  to  be  rapidly  coupled  to  its  overtone  as  well  as  to  the 
10%  level.  (See  Fig.  1)  Furthermore,  the  overtones,  and  the  stretch 
are  less  than  100  cm  *  apart  in  energy,  and  the  precise  coupling  among  these 
levels  has  long  been  unknown.  Huddleston  and  Heitz  measured  the  "ladder- 
climbing"  rate 

Ky 

2CO.(0110)  C0,(10%/0220/02%) 

2  >S)  2 

and  found  it  to  be  387  +  75  msec  *  torr  in  good  agreement  with  the  rate 

measured  by  Javan,  et  al.  ^  However,  due  to  various  problems  inherent  in  the 
fluorescence  emission  technique  that  was  used,  they  were  unable  to  perform  a 
complete  kinetic  analysis.  Poor  signal  to  noise  ratio  of  the  emission  out  of 
these  low  frequency  states,  limitations  in  detector  response  time,  and  the  low 
resolution  of  infrared  filters  make  this  type  of  experiment  particularly 
difficult.  Furthermore,  the  symmetric  stretch  is  infrared  inactive  and 
cannot  be  directly  observed  in  emission  studies.  Whereas  exact  information 
could  not  be  obtained,  Huddleston  and  Weitz  were  able  to  place  a  lower  limit 
of  about  1100  msec  *  torr  *  on  the  rate  of  energy  transfer  among  the 
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(10°0/0220/02°0)  states. 

Diode  laser  Infrared  double  resonance  Is  a  perfect  technique  to  obtain 
the  desired  Information  on  this  system.  The  high  resolution  (about  .001  cm  *) 
makes  it  possible  to  monitor  population  in  individual  rotational/vlbratlonal 
states  of  CO^*  By  selectively  monitoring  transitions  of  the  type 
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(011!) 

it  is  possible  to  observe  the  time  evolution  of  the  population  in  these  states 
after  perturbing  the  population  with  a  Q-switch  COj  laser.  Not  only  is  there 
no  question  as  to  which  state  is  being  observed,  but  by  "riding”  the  tran¬ 
sition,  the  strong  oscillator  strength  of  can  be  used  to  advantage  to  gain 
Information  on  the  modes  with  weak  oscillator  strength  and  to  directly  monitor 
the  Infrared  inactive  symmetric  stretch. 

Experimental : 

A  diagram  of  the  apparatus  is  shown  in  Fig.  2.  A  home  built  Q-swltch 
laser  was  used  to  pump  the  COj  sample,  using  either  the  10.6  micron  or  the  9.6 
micron  lines.  The  laser  pulse  shape  was  carefully  monitored  using  a  AuGe 
Infrared  detector;  the  pulse  width  was  never  more  than  1.5  microsec  FV  at  the 
10Z  amplitude  points.  The  COg  laser  beam  was  directed  through  the  sample  cell 
and  could  be  easily  aligned  by  means  of  an  NRC  optical  periscope;  the  final 
mirror  used  to  send  the  beam  through  the  cell  consisted  of  a  2  inch  Ge  beam 
splitter,  with  the  front  side  coated  for  80Z  reflectivity  at  10.6  microns  and 
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che  back  side  coated  for  antiref lection  at  10.6  microns*  Care  had  to  be  taken 
that  the  samll  portion  (about  20Z)  of  the  10.6  micron  beam  which  passed 
through  the  flat  rather  than  being  reflected  was  not  picked  up  by  the  probe 
beam  optics  and  sent  back  to  the  diode  laser.  A  4  to  5  micron  bandpass  filter 
was  used  to  block  the  10.6  or  9.6  micron  laser  beam  at  the  output  of  the 
sample  cell. 

The  probe  beam  at  4.3  microns  was  generated  by  an  infrared  diode  laser 
purchased  from  Laser  Analytics.  A  closed  cycle  helium  compressor  maintained 
the  diode  at  approximately  12K.  Current  and  temperature  were  controlled  by 
units  also  supplied  by  Laser  Analytics.  The  diode  laser  beam  was  collimated 
by  means  of  an  FI  CaF^  lens,  and  made  collinear  with  the  Q-switch  CO^  laser 
beam  by  passage  through  the  flat  (although  the  back  of  the  flat  was  AR  coated 
for  10  microns  this  coating  had  little  or  no  effect  on  the  4.3  micron  beam). 
After  passing  through  the  cell,  the  probe  beam  was  focussed  onto  the  slits  of 
an  F4.5  Bausch  and  Lomb  monochromator  by  means  of  an  F4.5  ZnSe  lens.  The  low 
power  (about  1  mwatt)  probe  signal  was  detected  by  means  of  an  InSb  photo¬ 
voltaic  detector  with  a  3  am  square  element.  The  response  time  of  the 
detector/ampllf ier  assembly  was  approximately  1.5  microsec.  The  signal  was 
preamplified  with  a  Perry  Amp,  model  730,  further  amplified  by  a  Tektronix 
Am502  1  Hz  to  1  MHz  wideband  amplifier,  digitized  with  a  Biomatlon  8100  tran¬ 
sient  digitizer,  and  signal  averaged  with  a  Nicolet  1170  signal  averager. 

After  averaging  the  signal  was  sent  to  a  PDP11/34  or  a  Vax  11/780  computer  for 
further  analysis.  Spectral  scans  were  taken  by  chopping  the  signal  at  about 
400  Hz  and  using  a  Keithley  Lock-in  Amplifier  whose  output  went  to  a  Houston 
2000  chart  recorder.  Relative  frequency  assignments  were  aided  by  use  of  a 
1  inch  Ge  etalon  whose  fringe  spaclngs  were  on  the  order  of  .047  cm  *  per 
fringe.  Absolute  assignment  of  one  point  was  necessary  in  order  to  convert 


140 


> 

relative  frequency  to  absolute  frequency*  This  was  done  by  comparison  of  the 

t 

relative  peak  positions  and  Intensities  with  those  reported  by  the  Air  Force 

Cambridge  Research  Laboratories  Atmospheric  Absorption  Line  Parameters 

Compilations. Additional  confirmation  was  possible  by  comparing  a  ^C0_ 

1  12 

spectrum  with  the  C(>2  spectrum  taken  over  the  same  current  range. 

Results: 

Transient  absorption  signals  were  observed  for  the  transitions 

(10°0) 

(0220) 

(02°0) 

(01l0) 

* 

using  both  the  9.6  micron  and  the  10.6  micron  CO^  laser  lines  as  pump  sources. 
In  the  case  of  the  (10^0)**  (10^1)  P(20)  transition  with  the  10.6  micron  pump, 

♦ 

as  well  as  the  (02  0) (02  1)  transition  with  the  9.6  micron  pump,  care  was 
taken  to  ensure  that  the  pump  laser  was  not  directly  pumping  the  same 
rotational/vibrational  state  that  was  being  probed.  In  the  few  instances  in 

» 

which  the  pump  laser  was  directly  pumping  the  same  rotational  state,  rapid 

rotational  relaxation  occurred  in  the  fall  of  the  signal,  and  the  decay  curve 

was  more  difficult  to  analyze.  Fig.  (3A)  shows  a  typical  signal  observed  when 

pumping  with  the  10.6  micron  line  and  observing  the  (10^0)-'-  (10^1)  P(20) 

transition.  The  rise  time  is  detector  limited  at  all  pressures,  and  occurs 

^  within  the  laser  pulse,  corresponding  to  rapid  depletion  of  the  (10^0)  state 

as  population  is  shifted  to  the  (00^1)  level  with  the  Q-switch  laser.  The 

2  0 

fall  is  rather  fast,  and  corresponds  to  the  adjacent  levels,  (02  0),  (02  0), 


(10ul)  P(20) 
(022l)  P(17) 
(02°1)  P(20) 


(ori)  p(3i) 
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Time 


(10°0)  -*  (10°1)  Probe 
.5  torr  CO2 

Full  scale  .52  msec  (.02  msec/. 05  msec  dual  time) 

(0220)  ■+  (0221) 

.8  torr  CO2 

Full  scale  .824  msec  (.024  msec/. 08  msec  dual  time) 

(0220)  ->■  (02C1) 

.5  torr  CO2 

Full  scale  1.05  msec  (.05  msec/1  msec  dual  time) 
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as  well  as  (01*0),  sharing  their  population  with  the  now  deficient  (10^0) 

2 

level.  Observation  of  the  (02  0)  P(20)  level  with  the  9.6  micron  pump  results 
in  a  qualitatively  similar  signal. 

Observation  of  the  (0220)  -*■  (0221),  (02°0)  -*•  (02°1),  and  (01°0) -*(01ll) 
transitions  with  the  9.6  micron  pump  can  also  give  information  on  the  coupling 
among  these  four  states.  In  this  case  however,  the  rise  time  should  not  be 
occuring  within  the  laser  pulse,  but  rather  will  reflect  the  time  dependent 
coupling  of  the  states.  Depending  on  the  actual  coupling  among  the  levels  and 
the  level  probed,  the  observed  fall  may  or  may  not  have  information  on  the 
coupling. 

The  curves  were  computer  fit  using  a  nonlinear  least  square  fitting 
routine.  The  following  is  a  summary  of  the  results: 

10.6  micron  Pump: 

(10°0)  -(10°1)  Probe: 

The  detector  limited  rise  was  followed  by  a  fall  which  appeared  to  have 
three  components:  1)  an  extremely  rapid  fall,  also  detector  limited,  at  all 
but  the  lowest  pressures  (less  than  .1  torr);  2)  a  second  fall  was  measured  to 
have  a  rate  of  about  700  msec  *  torr  *;  3)  a  third  slow  fall  was  thought  to  be 
either  the  V-T,  R  rate  or  the  rate  of  diffusion  out  of  the  diode  laser  beam. 
All  the  signals  displayed  this  slow  fall.  While  this  rate  has  no  real 
significance  in  the  population  equilibration  of  the  four  levels  of  interest, 
it  had  to  be  subtracted  off  before  the  other  rates  of  interest  could  be 
analysed.  This  rate  will  be  ignored  in  the  description  of  the  other  signals. 
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(0220)  -*(0221)  Probe: 

The  signal  is  shown  in  Fig.  (3B).  Studies  at  low  pressures  indicate  that 
there  is  a  single  fast  rise  having  a  rate  of  2500  +  500  msec  1  torr  The 
fall  has  only  one  relevant  component,  a  small  amplitude  decay  (less  than  25Z 
of  the  total  signal),  having  a  rate  of  about  400  +  40  msec-*  torr-1. 

(02°0)  -*-(02°l)  Probe: 

The  signal  is  shown  in  Fig.  (3C).  Again,  there  is  only  one  discernable 
rise  and  no  fall  other  than  the  slow  V  -+-T, R  rate.  The  rise  was  the  slowest 
observed  having  a  rate  between  700  and  1000  msec  1  torr  1. 

9.6  micron  Pump: 

(10°0)  -(10°1)  Probe: 

The  signal  is  shown  in  Fig.  (4A).  Pressure  studies  from  .5  torr  down  to 
•3  torr  place  this  rate  at  1000  +  300  msec  1  torr  1.  Again  no  fall  was 
clearly  observed;  however  at  higher  pressures  an  extremely  small  amplitude 
"wave"  was  discernable  in  the  fall.  The  small  amplitude  of  this  portion  of 
the  signal,  along  with  the  noise  made  it  Impossible  to  fit.  This  fall  is 
thought  to  be  due  to  the  ladder  climbing  process. 

(0220)  -*-(022l)  Probe: 

The  signal  is  shown  in  Fig.  (4B).  The  single  fast  rise  has  a  rate  of 
2300  +  500  msec  1  torr  1.  Again  no  fall  was  observed  other  than  the  slow 
V-*-T,R  rate. 
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Time 


Fig.  4 


A) 

(10°0)  -*•  (10°1)  Probe 
1  torr  CO2 

Full  scale  .824  msec 

(.024  msec/. 8  msec  dual  time) 

B) 

(02°0)  (0221) 

1  torr  CO2 

Full  scale  2.05  msec 

(.05  msec/2  msec  dual  time) 

C) 

(02°0)  (02°1) 

.5  torr  CO2 

Full  scale  1.05  msec  (.05  msec/1  msec  dual  time) 
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(02°0)-*-  (02°1)  probe: 

Here  we  are  observing  the  state  that  Is  being  directly  pumped  by  the 
laser.  Thus,  all  klnetically  Interesting  processes  will  occur  in  the  fall, 
the  rise  being  detector  United  and  within  the  laser  pulse  width.  The  fall 
clearly  has  two  conponents:  1)  a  fast  conponent  which  is  detector  United  at 
all  but  the  lowest  pressure  studied;  2)  a  second  conponent  with  a  decay  rate 
of  390  +  40  nsec  ^  torr  \  and  a  moderate  amplitude. 

All  of  these  signals  should  have  three  conponents  due  to  the 
equilibration  of  four  levels  (see  kinetics  section),  as  well  as  the  V^T,R 
rate.  In  no  case  were  all  three  rates  clearly  discernable.  The  reason  for 
this  seealng  reduction  of  observed  number  of  rates  is  simply  an  amplitude 
problem:  one  (or  more)  of  the  amplitudes  nay  be  relatively  small  in  comparison 
to  the  other(s).  The  relatively  small  amplitude  together  with  the  fact  that 
all  three  rates  are  rather  close  in  magnitude  may  make  the  small  amplitude 
rates  difficult  to  detect.  This  problem  was  especially  true  at  the  extremely 
low  pressures  that  were  required  to  ensure  that  the  faster  rates  were  not 
detector  limited.  At  these  low  pressures,  the  noise  was  on  the  same  order  of 
magnitude  as  the  small  amplitude  of  some  signals. 

The  technique  of  adding  Ar  to  the  sample  cell  in  order  to  distinguish 

between  resonant  versus  non-resonant  rates  was  invaluable  in  helping  us  to 

(14} 

distinguish  the  ladder  climbing  rate.  The  only  two  states  in  which 

extensive  Ar  pressure  studies  were  done  were  the  (10%)  level  with  the  10.6 
micron  pump,  and  the  (02%)  level  with  the  9.6  micron  pump.  Qualitative 
results  were  obtained  for  some  of  the  other  levels,  and  the  results  were  quite 
interesting. 
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Argon  Studies: 


10.6  micron  pump,  (10^0)  Probe: 

The  pure  C02  signal  exhibited  only  two  rates,  an  extremely  fast  fall  (too 
fast  to  measure  under  the  conditions  used),  and  a  slower  fall  which  had  a  rate 
of  at  least  700  msec  *  torr  * .  At  .125  torr  C02>  addition  of  up  to  10  torr 
Ar,  this  intermediate  rate  of  700  steadily  increased.  However,  addition  of 
12  to  30  torr  Ar  caused  the  observed  rate  to  slow  down.  (See  Pig.  (5)).  What 
appears  to  be  occuring  is  that  the  Ar  causes  the  rate  to  speed  up  and  become 
indistinguishable  from  the  detector  limited  fast  rate,  while  the  small 
amplitude  ladder  climbing  rate  begins  to  appear.  However,  at  .125  torr  C02, 
this  slower  ladder  climbing  process  has  an  amplitude  too  small  with  respect  to 
the  noise  to  fit  well.  At  .2  torr  C02  and  30  torr  Ar,  the  ladder  climbing 
rate  appeared  as  380  msec  1  torr  1  of  C02. 

9.6  micron  pump,  (02^0)  Probe: 

The  pure  C02  signal  again  exhibited  only  two  rates:  an  extremely  fast  one 
as  well  as  a  “slow"  one  with  a  rate  of  380  +  40  msec  *  torr  Upon  addition 
of  Ar  the  slower  rate  remained  approximately  the  same.  This  fact  clearly 
demonstrated  that  the  rate  under  study  was  the  resonant  ladder  climbing 
process  < 

9.6  micron  pump,  (10^0)  Probe: 

The  pure  C02  signal,  which  exhibited  a  rise  of  about  700  to  1000  msec  * 
torr~*  and  no  perceptible  fall  other  than  the  V  -*-T,R  rate,  upon  addition  of 
between  10  and  40  torr  of  Ar  exhibited  a  detector  limited  rise,  followed  by  a 
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fast  fall  whose  rate  Increased  with  Increasing  Ar  pressure!  The  source  of 
this  Ar  pressure  dependent  rate  is  a  puzzle  since  the  pure  COj  signal  shows  no 
such  process  at  all. 

Kinetic  Scheae: 

Since  the  four  levels  in  question  are  fairly  well  Isolated  from  other 
levels  (see  Fig.  1),  the  system  can  be  treated  as  a  four  level  system  at  short 
times.  The  validity  of  this  approximation  is  further  assured  by  the  fact  that 
the  state  to  which  ?<te  population  is  being  pumped,  the  (00%)  fundamental, 
has  an  exceedingly  long  lifetime  (on  the  order  of  30,000  collisions) 

Since  the  slowest  rate  in  the  four  level  system,  the  ladder  climbing  rate,  is 
at  least  1000  times  faster  than  this,  the  four  level  system  will  have  equili¬ 
brated  with  Itself  before  it  begins  to  receive  population  from  the  (00%)  or 
other  higher  energy  levels.  The  ultimate  goal,  of  course,  is  not  only  to 
determine  the  individual  coupling  rates  of  the  levels,  i.e.,  the  kinetic  rate 
constants,  but  concomitantly  to  unambiguously  determine  the  coupling  scheme. 
Since  the  experiment  directly  measures  the  time  dependent  populations,  the 
rate  Information  that  is  measured  corresponds  to  convolutions  of  the  kinetic 
rate  constants.  In  order  to  deconvolute  the  measured  rates  into  kinetic  rate 
constants,  it  is  necessary  to  solve  a  set  of  coupled  differential  equations. 

To  set  up  the  equations  a  kinetic  coupling  scheme  must  be  assumed.  Two 
kinetic  schemes  were  tried,  and  through  comparison  of  both  rates  and  relative 
amplitudes,  one  scheme  was  choosen  as  the  best  fit  to  the  data.  The  two 
kinetic  schemes  are  shown  in  Fig.  6.  The  only  difference  between  the  two 
schemes  was  as  follows.  In  scheme  I,  the  fundamental  of  the  bend  was  assumed 
to  couple  directly  to  the  (02%)  state  and  only  indirectly  to  the  other  states 
through  the  (02%)  state.  In  scheme  II,  the  (01%)  level  was  assumed  to 
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Kinetic  Scheme  II 


Figure  6 


2 

couple  to  the  (02  0)  level  and  only  indirectly  with  the  other  two  levels*  Of 
course,  if  the  ladder  climbing  rate  is  slow  in  comparison  to  the  rate  of 
energy  transfer  among  the  upper  three  levels,  the  two  schemes  will  be  vir¬ 
tually  indistinguishable. 

Since  the  approach  to  solving  either  scheme  is  exactly  the  same,  only 
scheme  II  will  be  presented  in  detail  as  it  is  the  one  which  best  fits  the 
data.  The  following  equations  describe  the  coupling  of  kinetic  scheme  II: 


CO.CO^O)  +  C0,(01l0)  C0„(0220>  +  C0_(00°0) 

4  /  y  i.  L 

k31 

£ 

C02(02°0)  +  C02(00°0)  ^  C02(10°0)  +  C02(00°0) 

K42  ‘  C 

K 

C02(02°0)  +  C02(00°0)  C02(0220)  +  C02(00°0) 

K32 

K 

C02(0220)  +  C02(00°0)  C02(I0°0)  +  C02<00°0) 

K43 


Solution  of  the  coupled  differential  equations  for  this  kinetic  scheme  yields 
3  eigenvalues  which  correspond  to  the  three  observed  rates  of  the  system.  The 
initial  conditions  are  then  used  to  evaluate  the  eigenvectors  of  the  system 
which  correspond  to  the  observed  amplitudes.  Since  the  elgnvalues  do  not 
depend  on  the  initial  conditions,  the  three  observed  rates  will  be  the  same 
regardless  of  which  laser  line  is  used  to  pump  the  system;  however,  the 
amplitudes  do  depend  on  the  initial  conditions,  and  therefore  will  be 
different  for  the  10.6  pump  versus  the  9.6  pump. 
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Discussion/ Conclusion: 


The  three  measured  rates  were: 


390  + 

80 

-1 

-1 

msec 

torr 

-1 

-1 

1000  + 

300 

msec 

torr 

700 

-1 

-1 

2000  + 

msec 

torr 

The  errors  on  the  two  faster  rates  were  set  somewhat  arbitrarily  at  approxi¬ 
mately  30Z  of  the  measured  rates.  This  value  was  chosen  because  there  are 
always  large  systematic  errors  in  measuring  rates  from  fluorescence  curves 
where  extensive  deconvolution  is  required.  This  was  particularly  apparent  in 
the  case  of  the  two  faster  rates  since  the  pressure  range  over  which  they 
could  be  measured  was  limited  on  the  high  side  by  detector  response  time,  and 
on  the  low  side  by  poor  signal  to  noise  ratios. 

The  first  question  to  be  addressed  is  whether  the  two  coupling  schemes 
can  be  distinguished  at  all.  The  two  schemes  were  examined  and  their 
associated  rate  constants  were  evaluated  by  assuming  values  for  kinetic  rate 
constants  and  then  solving  the  coupled  equations.  The  resulting  rates  and 
amplitude  factors  were  compared  to  the  experimental  results  and  the  kinetic 
rate  constants  were  varied  until  there  was  good  agreement  with  the  experi¬ 
mental  data.  Although  the  data  does  not  absolutely  eliminate  Kinetic  Scheme 

I,  the  results  do  seem  to  favor  Kinetic  Scheme  II.  In  scheme  II  the 

2 

fundamental  couples  directly  to  the  (02  0)  levels  with  which  it  is  nearly 
resonant,  rather  than  the  Fermi  mixed  (10^0)  and  (02^0)  levels.  The  kinetic 
rate  constants  which  appear  to  best  fit  the  observed  data  are: 
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K-.  1  400  asec  *  torr  * 

-  -1  -1 

Kjj  “  600  asec  torr 

K-.  •  600  msec  1  torr  * 

34  -1  -1 

=  100  msec  torr 

In  most  fluorescence  studies  direct  comparisons  of  amplitude  ratios  Is 

difficult  If  not  Impossible  because  of  the  large  errors  associated  with  their 

measurement*  Fortunately,  the  10*6  and  the  9.6  micron  pump  yield  different 

amplitude  factors  which  makes  qualitative  comparisons  possible.  The  fact  that 

2 

the  ladder  climbing  rate  appeared  quite  clearly  In  the  (02  0)  state,  but  not 
0  2 

in  the  (10  0)  state  or  the  (02  0)  state  when  pumped  at  10.6  microns  is  a  clue 

2 

that  the  fundamental  appears  to  preferentially  couple  to  the  (02  0)  state 
rather  than  either  the  (10%)  or  the  (02%). 


*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78ER04940. 
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D.  COMPETITION  BETWEEN  INELASTIC  AND  REACTIVE  SCATTERING:  A  STUDY  OF  N.O  AND 

<V  2 

(J.  Subblah,  G.  W.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Studies  of  energy  transfer  dynamics  using  hot  species  produced  by 
photofragmentation  in  the  gas  phase  are  valuable  because  they  allow  simulation 
of  the  chemistry  of  the  upper  atmosphere  to  determine  the  role  of  critical 
species  in  the  oxygen  cycle;  they  contribute  to  a  better  understanding  of  the 
dynamics  of  thermallzation  of  hot  species,  including  possible  competition 
among  elastic,  inelastic  and  reactive  channels  on  potential  energy  surfaces; 
they  can  be  used  to  measure  the  dependence  of  reaction  cross  section  on  the 
relative  translational  energy;  and  they  provide  valuable  tests  of  theories  of 
inelastic  scattering  and  translation-to-vibratlon  (T-V)  energy  transfer. ^ 

We  have  studied  collisions  of  hot  atoms  in  the  gas  phase  produced  by  the 

photolysis  of  suitable  precursors  at  193  nm.  The  first  system  studied 

involves  the  photolysis  of  N^O,  which  is  interesting  because  of  its  similarity 

to  the  linear  synsnetrlc  polyatomic  CO^,  whose  excitation  by  hot  atoms  has  been 
(2) 

studied.  Although  CO^  is  interesting  because  of  its  simplicity  and  huge 
oscillator  strength,  its  symmetric  stretching  mode  is  IR  inactive  and  this 
makes  it  difficult  to  compare  the  relative  degree  of  excitation  of  the  two 
stretching  modes.  N^O,  on  the  other  hand,  being  a  linear  but  unsymmetric 
trlatomlc,  does  not  have  this  problem,  and  is  a  good  test  case  to  determine 
the  posslblity  of  mode-selective  vibrational  excitation  by  a  translatlonally 
hot  projectile. 

The  first  experiments  were  performed  on  the  HBr-N^O  system.  The  relevant 


processes  are: 


hv  (X“193nm)  +  HBr  - >  H*  +  Br  Hot  Atom  Generation 

H*  +  N20  - >  H  +  N20  (v1)  T-V  Excitation 

H  +  HBr  - >  H2  +  Br 

H*/H  +  N20 - >  OH  +  N2  Reaction 

(AH  ■  -62.4  Real /mole) 
k  »  4x10  **  ccmol  *s  *  at 
20,000°K,  the  temperature  of  the 
hot-hydrogen  atoms 

H/H*  +  N20  - >  NH  +  NO  (AH  -  +30.8  Kcal/mole) 

Irradiation  of  HBr-N20  mixtures  produced  strong  IR  emission  from  the  N-N 

stretch  at  2154  cm  *  (see  Fig.  1),  while  the  emission  from  the  N-0  stretch  at 

-1  (3) 

1265  cm  is  very  weak.  Similar  work  by  Quick  et.al.  has  shown  that  the 

high  frequency  mode  is  slgnlf icantly  more  populated  than  the  v^  mode.  This 

result  is  totally  unexplainable  by  the  refined  impulse  approximation  approach 

to  T-V  transfer  which  predicts  a  more  favorable  excitation  of  the  lower 

frequency  mode. 

Irradiation  of  N20  alone  also  gives  rise  to  IR  fluorescence  from 
vlbratlonally  excited  N20,  though  the  signal  is  much  weaker  than  in  the  case 
of  HBr-N20  mixtures  (see  Fig.  2).  The  mechanism  is  postulated  to  be: 


1) 

N20  +  hv  (X 

-  193  nm) 

- >  N2  +  0  (lD) 

Photolysis 

2) 

0(lD)  +  N20 

- > 

0  +  N20  (v1) 

T-V  transfer 

3) 

0(lD)  +  N20 

- > 

2N0  (AH  »  -80.95  Kcal/mole) 

Reaction 

kj  ■  10***  ccmol  *s  * 

- > 

N2  +  02  (AH  -  -124.05  Kcal/mole) 

kj  ■  10***  ccmol  *s  *) 
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Fig.  1  IR  fluorescence  emission  from  N2O 
excited  by  irradiation  of  1:1 
mixtures  of  RBr  -f  N£0  at  193  nm 

total  pressure  :  0.5  torr 

time  base  :  100  ys  full-scale 

detector  :  InSb 

filter  :  4  to  5  micron 
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Fig.  2  IR  fluorescence  emission  from  N£0 

excited  by  irradiation  of  neat  N2O  at  193  nm 

total  pressure  :  0.5  torr 

time  base  :  10  ms  full-scale 

detector  :  InSb 

filter  :  4  to  5  microns 
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» 

As  In  the  case  of  HBr-N^O  mixtures  it  was  again  found  that  the  high 
I  frequency  N-N  stretch  was  much  more  strongly  excited  than  the  lower  frequency 

N-0  stretch. 

Fig.  2  shows  the  IR  emission  of  the  mode  of  Nj,0  when  the  neat  gas  is 
t  irradiated  at  193  nm.  One  of  the  most  obvious  differences  from  the  same 

emission  in  HBr-NjO  mixtures  is  that  the  risetime  is  not  detector-limited,  but 
is  rather  slow.  A  pressure-dependent  study  shows  that  the  rate  constant  on 
f  the  rise  is  about  7 (ms  torr)  *  corresponding  to  about  3900  collisions.  It  may 

be  that  this  excitation  corresponds  to  energy  transfer  from  (vibratlon-to- 
vlbration  energy  transfer)  or  from  0(*D)  (electronlc-to-vlbrational  energy 
transfer).  Further  work  is  required  to  verify  this  mechanism. 

As  Eq.  3  shows,  one  possible  excitation  channel  Involves  reactive 
encounters.  This  was  monitored  by  looking  for  vibratlonally  excited  NO,  which 

*  is  observed  (see  Fig.  3).  At  the  photolysis  wavelength  used,  it  is 
energetically  impossible  for  N20  to  dissociate  to  N  +  NO,  thus  the  NO  is 
probably  formed  by  the  reactive  channel  (Eq.  3).  This  reaction  is  known  to 

^  occur  at  gas  kinetic  rates,  accounting  for  the  fast  rise,  as  well  as  by  energy 

transfer  from  other  species  (slow  rise).  This  is  one  example  of  a  case  where 
reactive  scattering  actually  predominates  over  inelastic  scattering. 

^  Another  system  studied  was  C12~CH^,  which  can  provide  information 

relevant  to  the  possible  effects  of  the  presence  of  chlorine-containing 
species  of  natural  or  anthropogenic  origin  in  the  stratosphere,  and  their  role 

*  in  the  ozone  balance.  The  basic  features  are  similar  to  those  for  HBr-N20  and 
consist  of  the  following  steps: 

*  Cl2  +  hv  (X  -  193  nm)  - >  2C1*  (15000  cm  *)  Photodissociation 

Cl*  +  CH^  - >  Cl  +  CH^  (v£)  Inelastic  scattering 


« 
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3  IR  fluorescence  of  NO  following 

irradiation  of  neat  N£0  at  193  nm 

total  pressure  :  0.5  torr 

time  base  :  20  ms  full-scale 

detector  :  InSb 

filter  :  5.362  micron  narrow  band 


Reactive  channel 


> 


I 


t 


t 


f 


V 


I 


Cl*  +  CH4  - >  HC1+  +  CH3 

(AH  ■  -1.186  Real /mole) 

Cl*  +  CH^  - >  CHjCl^  +  H  Reactive  channel 

(AH  »  +21.4  Kcal/mole) 

Both  these  channels  are  accessible  because  of  the  excess  energy  of  the  hot 
chlorine  atom  produced  in  the  photolysis.  1R  fluorescence  was  monitored  from 
both  1R  active  fundamentals  of  CH^  a  degenerate  stretch  at  3019  cm  *  and  a 
degenerate  deformation  at  1306  cm  *  to  determine  the  distribution  of  energy 
among  different  modes  of  CH^.  Some  of  the  more  important  observations  are: 

(1)  Under  comparable  conditions,  the  IR  emission  is  much  stronger  in  the 
C^-CH^  system  than  in  the  HBr-CH^  system;  (2)  Strong  IR  emission  was  observed 
in  the  3  to  4  micron  spectral  region  when  5:1  mixtures  of  Cl^-CH^  were 
irradiated  at  193  nm  (see  Fig.  4).  This  emission  could  be  either  from 
vibrationally  excited  HC1  or  from  vibratlonally  excited  CH^.  Further  work 
will  be  done  to  assign  the  emission  by  wavelength  resolved  detection  using 
either  a  circular  variable  filter  or  a  cold  gas  filter. 

Future  experimental  work  will  Include  the  use  of  tunable  excitation 
sources,  better  wavelength-resolved  detection,  and  precursor  variation  to 
obtain  a  better  understanding  of  the  momentum  and  velocity  dependence  of  T-V 
energy  transfer. 


*Thls  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78ER04940. 
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LR  emission  following  irradiation  of 
1:5  mixtures  of  CH^rC^  at  193  nm 

total  pressure  :  0.8  torr 

time  base  :  2  ms  full-scale 

detector  :  Au:Ge 

filter  :  3  to  4  microns 


II 


» 


» 
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» 


» 
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E.  VIBRATIONAL  ENERGY  RANDOMIZATION  IN  S02  FOLLOWING  ELECTRONIC  EXCITATION* 

(B.  Lawrence,  P.  Beeken,  G.  W.  Flynn) 

(JSEP  work  unit  7,  1982-1985) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Despite  much  work  in  recent  years,  the  chemical  and  physical  processes 

of  electronically  excited  S02  are  not  well  understood.  The  physics  is  of 

particular  Interest  because  this  molecule  is  an  example  of  an  intermediate 

case  in  radiationless  transition  theory.  The  physics  and  the  chemistry  are, 

of  course.  Interrelated.  Spontaneous  radiationless  transitions  occur  in  large 

molecules  when  excited  to  an  electronic  state  which  is  coupled  in  some  way  to 

another  electronic  state  with  a  high  density  of  states  at  the  same  energy. ^ 

In  small  molecules,  the  density  of  states  is  not  large  enough  to  allow  for 

this  crossover  to  occur.  Collision-induced  radiationless  transitions  do  occur 

in  cases  such  as  S02  where,  perhaps,  colllsional  broadening  of  the  vlbronlc 

levels  increases  the  number  of  accepting  states  available  so  that  the 

(2) 

probability  of  curve  crossing  also  Increases.  If  an  intermediate-case 
molecule  is  electronically  excited  at  a  pressure  high  enough  for  the  molecule 
to  undergo  collisions  before  it  radistes,  it  may  cross  over  to  another  excited 
electronic  state  or  to  high  vibrational  levels  of  the  ground  state.  In  any 
case,  the  molecule  will  eventually  relax  to  the  ground  state  -  probably  with 
vibrational  excitation. 

One  might  ask  whether  the  excess  vibrationsl  energy  remaining  in  the 
molecule  when  it  returns  to  the  ground  state  will  be  randomly  distributed 
among  the  normal  modes  or  whether  one  mode  will  be  preferentially  excited  over 
the  others.  The  answer  to  this  question  should  give  some  information  about 
the  electronic  and  vibrational  relaxation  processes  occurring  after  the 
initial  excitation. 


163 


We  have  been  exciting  SO^  electronically  and  observing  tlae-resolved 

eaisslon  froa  the  aode  of  the  ground  state.  Vibrational  relaxation  of  the 

(3) 

ground  state  of  SOj  has  previously  been  studied  in  this  laboratory  by 
exciting  the  aode  and  observing  emission  from  the  aode.  Because  the 
rates  of  interaode  vibrational  energy  transfer  are  now  known,  observed 
emission  from  one  aode  will  yield  Information  about  the  initial  conditions  of 
vibrational  excitation. 

S(>2  has  three  major  absorption  regions  between  200  nm  and  400  nm  (see 
Pig.  1).  A  weak  absorption  between  390  nm  and  340  nm  has  been  assigned  to  a 

3 

spin-forbidden  transition  froa  the  ground  state  to  the  state.  This  region 
is  thought  to  be  perturbed  by  the  presence  of  two  other  triplet  states.  A 
second,  somewhat  stronger  absorption  region  froa  340  nm  to  260  nm,  called  the 
first-allowed  band,  has  been  assigned  to  two  overlapping  transitions  from  the 
ground  state:  an  electric-dlpole-allowed  transition  to  the  ^Bj  state  and  an 
electric-dlpole-forbidden  but  vibronically  allowed  (with  excitation) 
transition  to  the  ^2  state.  These  two  states  are  coupled  to  each  other  and 
to  the  lower  triplet  states.  Therefore,  this  region  is  very  complex.  A 
third,  very  strong  absorption  region,  called  the  second  allowed  band,  starts 
at  23S  nm  and  is  assigned  to  the  allowed  transition  from  the  ground  state  to 
the  *^2  state. 

In  these  initial  studies,  we  have  been  electronically  exciting  SO2  into 
the  first  allowed  band  using  the  frequency-doubled  output  of  a  Lambda-Physik 
exclaer- pumped  dye  laser  using  Rhodamlne  6G  dye.  The  SO2  has  been  at 
pressures  of  500  mtorr  to  8  torr.  The  time-resolved  7.4iia  emission  froa  the 
vibrational  aode  is  detected  through  a  7-8 ia*  bandpass  filter  by  a  AuGe 
(77K)  detector.  The  signal  is  collected  using  a  transient  digltlzer/dlgltal 
signal  averager  coabination  and  transferred  to  a  PDP11  coaputer  for  later 
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Energy  level  diagram  for  low-lying  electronic  states  of  SO2 
(EKC  Lee  and  G.  Loper  in  "Radiationless  Transitions"  (SH  Lin,  ed.) 
p.22.  Academic  Press,  New  York) 


Low  resolution  absorption  spectrum  of  SC^ 

(e  in  units  of  £mol-1cm~l  to  base  10) 

From  Strickler  and  Howell,  J,  Chem.  Phys.  49,  1949  (1969) 


analysis 


Preliminary  results  show  a  fast  exponential  rise  and  a  slower  exponential 

fall  in  the  emission,  (see  Fig.  2)  Data  taken  at  several  different 

pressures  was  fit  to  a  double-exponential  decay  curve,  after  the  subtraction 

of  a  thermal  tail.  Using  decay  curves  obtained  at  an  excitation  wavelength  of 

305  nm,  Stern-Vollmer  plots  of  rate  vs.  pressure  have  rates  of  80.7  (msec  * 

torr  *)  for  the  rise  and  3.5  (msec  ^  torr  *)  for  the  fall.  Comparisons  with 

(3) 

ground  state  energy  transfer  studies  from  this  laboratory  show  that  the 

rates  are  roughly  comparable  to  those  obtained  with  direct  laser  excitation  of 

the  ground  state  mode.  This  suggests  that  the  fast  rise  corresponds  to 

the  population  of  the  mode  by  the  rapid  equilibration  of  the  and 
(3) 

modes.  The  slower  fall  appears  to  be  the  result  of  vibrational  relaxation 

(3) 

through  the  bending  mode  to  the  translational  bath. 

(4) 

In  previous  experiments,  S02  was  excited  at  248  nm.  (This  is  50  nm 
higher  in  energy  than  the  present  experiments  in  a  region  between  the  second 
and  third  absorption  bands.)  The  resolved  emission  showed  a  fast  exponential 
fall  followed  by  a  slower  exponential  fall.  The  rates  were  roughly  comparable 
to  the  rise  and  fall  observed  at  the  longer  wavelength.  These  results  seem  to 
Indicate  that  at  305  nm  excitation  the  mode  is  initially,  preferentially 
excited,  while  at  248  nm  excitation  the  mode  is  initially  excited.  It 
should  be  emphasised  that  these  results  and  their  interpretation  are 
preliminary. 

In  the  immediate  future  we  plan  to  perform  the  same  experiments  using 
different  excitation  wavelengths  while  looking  at  emission  from  the  mode 
(with  a  AuGe  detector  and  an  8-9u  filter)  as  well  as  from  the  mode.  Using 
an  excimer  laser  at  308  nm  we  would  be  able  to  use  much  lower  pressures  of 
S02  in  order  to  learn  more  about  the  initial  rise  in  population  of  the 


Fig.  2  SO^Cv^)  emission 

pso,  ■  8-°  torr 

Nshots  ■  400 

Time  base  >0.2  msec  (full  scale) 
Excitation  \  -  305.0  nm 
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vibrational  nodes  and  to  rule  out  other  processes  such  as  radiational 
relaxation  or  Raman  scattering.  A  careful  study  of  the  variation  of  the  decay 
rates  with  excitation  wavelength  for  the  first  allowed  band  should  be  done. 

We  would  also  like  to  try  exciting  into  the  second-allowed  and  the  triplet 
region.  The  highest-energy  region  would  be  especially  Interesting  because  it 
is  much  less  complicated  and  the  results  should  therefore  be  easier  to 
interpret. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  contract  DE-AC-02- 
78ER04940. 

(1)  J.  Jortner,  S.  Rice,  and  R.  M.  Hochstrasser,  Adv.  Photochem.  T_,  149 
(1969). 

(2)  K.  F.  Freed,  Chem.  Phys.  Lett.  37,  47  (1976). 

(3)  J.  L.  Ahl,  Ph.D.  Thesis,  Columbia  University,  1982. 

(4)  Columbia  Radiation  Laboratory  Report  31,  1980-1981,  pp.  104-106. 
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F.  LASER  VS.  THERMAL  EXCITATION  OF  SF,  IN  A  MOLECULAR  BEAM  AND  THE  USE 
OF  LASER  EXCITED  SFfi  IN  A  CHEMI- IONIZATION  REACTION  WITH  X2* 

(JSEP  work  unit  7,  1982-1985) 

(G-  Spector,  G.  W.  Flynn) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

During  the  last  year,  work  was  undertaken  to  complete  two  projects,  the 

chemi-ionlzation  reaction  of  X2  and  SFg  studied  in  a  crossed  molecular  beam 

machine  and  the  study  of  laser-  vs.  thermal-excitation  of  SF,  molecular  beams. 

o 

A  manuscript  has  been  prepared^ ^  for  publication  which  describes  the  laser 
heating  experiments. 


(1)  Chemi-ionlzation 

Earlier  work  in  this  laboratory  had  demonstrated  that  internal  energy  can 

be  deposited  in  a  supersonic  molecular  beam  of  SF^  with  a  lower  power,  cw  CO,, 

(21 

laser  [P(16)10.6pm,  6.5W].  Later  studies  showed  that  laser  irradiation 

inside  a  transparent  nozzle  enhanced  that  internal  energy  of  SF,  molecules 

o 

(3) 

significantly.  In  early  1982,  an  application  study  commenced,  using  the 

laser-excited  SF^  beam  with  effusive  beams  of  K  and  Na.  E.W.  Rothe,  of  Wayne 

(31 

State  University,  had  *;4ied  several  M  +  SF^  systems  '  determining  the 
relative  chemi-ionlzation  cross  sections  of  electronically  excited  M*  vs. 
ground  state  M.  We  have  studied  a  similar  reaction,  but  with  a  vibrationally 
excited  SFg  beam. 

An  electron  multiplier  was  employed  to  detect  cations  at  the  scattering 
center.  Unfortunately,  the  multiplier  detects  all  positive  ions,  so  a 
definitive  prediction  of  particular  product  formed  is  Impossible.  One  can, 
however,  ascertain  from  Fig.  1  that  only  the  formation  of  K^F*  *•» 

energetically  favored,  the  positive  ion  is  formed  with  dimer  (K^)*  There  is 
no  stable  ionic  product  formed  with  monomer  (X).  The  apparatus  used  in  these 
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experiments  is  described  elsewhere'.^ 

Although  the  availability  of  mass  selective  detection  would  provide  more 
information,  it  is  the  relative  (laser  excited  SP^  vs.  ground  state  SFg) 
reaction  cross  section,  rather  than  the  identity  of  the  products,  in  which  we 
are  most  Interested.  More  than  thirty  scattering  experiments  were  performed 
with  several  modifications  and  adjustments  in  procedure,  including  changes  In 
positioning  of  the  electron  multiplier,  changes  in  collimation  (of  the  K  beam) 
and  the  insertion  of  a  repeller  plate  and  a  scattering  chamber. 

Fig.  2  shows  the  ratio  Il/Iq,  the  laser  on/laser  off  electron  multiplier 
signal  of  reaction  product  as  a  function  of  the  irradiation  position  inside 
the  noxsle  capillary.  From  these  values,  using  the  equation: ^ 


a 


rel 


tL  0  L  /  E°  X1^2 
I  Ap  <v  >  I  tr  1 

1°  APL  <v°>  \e"  / 


(1) 


the  relative  cross  section,  a  can  be  determined  and  is  plotted  in  Fig.  3 

as  a  function  of  average  Internal  energy,  As  Elnt  increases  (laser 

on),  the  cheml-ionlzation  cross  section  decreases.  This  is  a  most  Interesting 

result;  while  there  are  examples  of  reactions  with  energy  barriers  (endoerglc 

or  activation)  whose  cross  sections  increase  with  increasing  internal 
(8)  (9) 

energy'  ' ’ '  and  exothermic  processes  in  which  o  decreases  with  kinetic 
energy, few  if  any  exoergic  reactions  have  been  reported  with 
decreasing  cross  section  resulting  from  internal  excitation. 


(2)  SFj  thermal  heating  experiments 

In  earlier  experiments  we  determined  the  relative  effect  of  1)  laser 
excitation  and  2)  simple  thermal  heating  on  the  internal  energy  of  an  SF^ 
molecular  beam.  Energy  fluxes  (W  n<v>E{ q  •  FEtQt)  were  measured  using  a 


f 


171 


IRRADIATION  POSITION  (mm) 


Figure 


liquid  helium-cooled  bolometer;  number  densities  (n)  end  aversge  bean 
velocities  (<v>)  were  found  using  the  quadrupole  mass  filter.  An  ion  gauge 
flux  meter  provided  an  alternate  method  for  determining  the  beam  flux,  F. 


The  internal  energy  of  SF^  at  room  temperature,  Ejnt,  was  determined 
relative  to  Ar,  a  monatomic  gas  with  no  internal  energy,  using  bolometer,  flux 
meter  and  TOP  results.  Energies  at  elevated  temperatures  were  then  found 


using  the  following  equation: 


.(D,(12) 


wtA  .  i  -<»j>  +  (Et„t  +  Ea)<y 

J»'vo>  +  (4t  +  V<V 


(2) 


Internal  energies  of  laser-excited  SFg  beams  were  found  in  the  same  manner. 

These  experiments  showed  that  CO^  laser  excitation  of  SF^  beams  enhances 

the  internal  energy  far  more  than  simple  heating  of  the  nozzle.  There  was, 

however,  confusion  in  that  the  Eint  values  obtained  In  simple  thermal  heating 

experiments  were  considerably  less  than  those  obtained  using  an  enthalpy 
(13) 

balance  equation:' 

T 

/*nz  12  3 

4  kT  +  /  C  (vib)dT  *  ^  m<vS  +  ~  kT  +  E.  (3) 

nz  Jq  v  2  2  s  lnt 


Experiments  were  performed  to  examine  this  discrepancy.  First,  TOF 
measurements  of  heated  Ar  beams  were  studied;  they  concurred  with  previous 
results.  Then,  SFg  thermal  TOF  studies  were  performed,  indicating  different 
velocities  at  elevated  temperatures  from  those  measured  earlier.  It  was 
discovered  thst,  while  other  measurements  (W,n)  were  taken  while  the  nozzle 
was  ’’cooling’*  in  heated  experiments,  the  SF^  TOF  parameters  (vg,  o^)  had  been 
measured  during  "heating"  runs.  The  latter  suffer  from  large  uncertainties  in 
temperature  because  of  a  rapid  change  in  temperature  upon  heating.  During 
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"cooling"  runs  the  teaperature  changes  aore  slowly.  The  new  SF^  velocities 
(v0)  are  plotted  in  Fig.  4  as  a  function  of  nossle  teaperature.  Similarly, 
the  width  paraater  (o^)  is  shown.  These  velocities,  considerably  slower  than 
those  aeasured  earlier,  enhanced  soaewhat  the  values  of  Internal  energy  (cf ., 
Eq.  2). 

In  a  related  experiaent,  the  entire  oven  assembly  was  heated  as  well  as 
the  nozzle  to  extend  the  residence  tlae  of  the  SF,  molecules  in  a  heated 

D 

chamber.  It  had  been  hypothesized  that  the  molecules  might  not  have 

equilibrated  to  the  nozzle  teaperature.  There  was ,  however,  no  difference  in 

the  TOF  paraaeters  with  and  without  the  heated  oven,  indicating  that  SF^ 

molecules  had  Indeed  reached  T  • 

nz 

Once  the  SF^  velocities  had  been  measured,  energies  were  calculated. 

Fig.  5  shows  the  internal  energy  F^nt»  the  kinetic  energy  and  the  total 

energy  Ftot»  graphed  as  a  function  of  nozzle  temperature.  Elnt  remains 
essentially  constant  in  the  range  300-400K,  then  Increases  by  ~45Z  in  the 
range  400-544K.  These  Internal  energies  can  be  compared  to  those  determined 
using  the  enthalpy  conservation  equation,  Eq.  3.  The  two  curves  are  shown  in 
Fig.  6  as  a  function  of  nozzle  teaperature.  There  is  still  a  wide  discrepancy 
between  thea  but  there  is  no  doubt  that  laser  excitation  provides  far  aore 
Internal  enhancement  than  does  thermal  heating.  At  the  optimum  laser 
irradiation  position  inside  the  nozzle  capillary,  ~0.38  eV  molec  *  are  found 
in  Internal  degrees  of  freedom,  while  at  544K  (thermal  heated),  only  ~0.03eV 
aolec  1  of  internal  energy  are  measured  with  similar  translational  energy 


distributions 


IRATURE  (K) 


*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  CHE  80-23747  and  the  Department  of  Energy  under  Contract  DE-AC-02- 
78ER04940. 
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(12)  <v  >  measured  from  TOF.  E  is  heat  of  adsorption  of  molecules  on 
bolometer  surface  (taken  to  be  -AH  ). 

I  vap 

(13)  Adapted  from  references  1  and  2.  T  is  stream  temperature,  measured  from 

TOF  parameter  a  .  9 
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V.  PICOSECOND  ENERGY  TRANSFER  AND  PHOTOFRAGMENTATION  SPECTROSCOPY 

A.  PHOTOCHEMISTRY  AND  PHOTOPHYSICS  OF  DIPHENYLCARBENE* 

(Y.  Hang,  F.  Novak,  E.  Sitzmann,  J.  Langan ,  K.  B.  Eisenthal) 

( JSEP  work  unit  8,  1982-1985) 

(Principal  Investigator:  K.  B.  Eisenthal  (212)  280-3175) 

In  recent  years,  much  theoretical  and  experimental  work  has  centered  on 
attempting  to  understand  the  chemical  and  physical  properties  of  divalent 
carbon  compounds  (carbenes),  an  Important  class  of  short-lived  and  reactive 
intermediates.^^  Of  particular  Interest  is  the  presence  of  an  unshared  pair 
of  electrons  which  results  in  two  low-lying  electronic  spin  configurations, 
singlet  and  triplet.  EPR  spectroscopy  has  been  extremely  useful  in  establish¬ 
ing  the  ground  state  multiplicity  of  most  aryl  and  alkyl  carbenes  to  be  the 
(2) 

triplet.  Valuable  structural  Information  for  these  carbenes  is  also 

(2) 

available  from  EPR  and  ENDOR  experiments,  but  little  is  known  concerning 

the  energetics  and  dynamics  between  the  two  spin  states.  An  Important  and 

controversial  aspect  of  this  is  the  magnitude  of  the  singlet-triplet  energy 

gaps.  Numerous  calculations  of  singlet-triplet  energy  splittings  have 
(3) 

appeared;  however,  little  quantitative  experimental  information  is 

available.  Recent  reports  of  the  experimental  determination  of  the  singlet- 

triplet  energy  gap  in  the  parent  methylene  have  placed  the  value  at  <V8 
(A) 

kcal/mole;  '  however,  the  corresponding  value  for  other  carbenes  has  not  been 
determined. 

In  addition  to  the  dearth  of  experimental  information  concerning  singlet- 
triplet  energy  gaps,  little  is  known  about  the  dynamics  of  the  singlet-triplet 
Interconversion.  Generation  of  a  ground  state  carbene,  either  photochemically 
or  thermally  from  an  appropriate  precursor  should  lead  initially  to  the 
singlet  state  because  of  spin  conservation  restrictions.  The  singlet  may  then 
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undergo  reaction  or  intersystem  crosa  to  the  ground  state  which  subsequently 
t  reacts.  As  a  result,  interpretation  of  carbene  chemistry  is  often  complicated 

by  products  originating  from  both  spin  states.  It  is  Important,  therefore,  to 
determine  the  rate(s)  of  Interconversion  of  the  spin  states  and  in  addition 
t  the  rate  constants  for  reactions  of  each  with  various  substrates. ^ 

He  have  undertaken  an  in-depth  study  of  diphenylmethylene  (DPM)  and  wish 
to  report  the  determination  of  the  absolute  rate  constants  for  singlet-triplet 
4  interconversion  and  the  energy  spacing  between  these  two  states  for  diphenyl- 

methylene  by  extension  of  the  pioneering  work  of  Closs  and  Rabinow.^ 
Diphenylmethylene  is  a  particularly  attractive  candidate  for  study  becaus'  * 
f  should  be  prototypic  of  arylcarbenes  and  furthermore,  it  has  been  shown  t  «s 

amenable  to  flash  spectroscopic  techniques. 

We  have  combined  competitive  experiments,  picosecond  and  nanosecond 
4  flash  spectroscopy,  and  triplet  sensitization  experiments  to  evaluate  the 

rates  of  singlet-triplet  Interconversion  of  DPM  in  acetonitrile  at  ambient 
temperature.  From  these  rate  constants,  K£  and  AGgT  are  determined  for 
*  DPM.*7) 


The  sequence  of  energy  decay  steps  following  excitation  of  dlphenyldiazo- 
methane  with  an  ultraviolet  light  pulse,  ultimately  leading  to  diphenylcarbene 
in  its  triplet  ground  state,  we  find  can  be  written  as: 
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1  3 

The  spin  conversion  from  DPC  -*•  DPC  occurs 
(In  collaboration  with  Professor  N.J.  Turro 
equilibrium  constant 


in  llOps.  From  measurement 
and  his  group)  we  obtain  the 


of  \s 


K  -  7~  -  (2.2  +  0.6)  x  103. 

TS 

The  free  energy  difference  at  2S°C  is  calculated  to  be  4.6  +  0.2  kcal/mole. 
Assuming  that  the  entropy  difference  between  the  singlet  and  triplet  states  is 
due  chiefly  to  the  multiplicity  difference  we  obtain 

AH  -  AG  +  TAS  -  3.9  +  0.2  kcal/mole. 

We  therefore,  have  succeeded  not  only  in  obtaining  the  rate  of  singlet-triplet 
equilibration,  but  also  an  estimate  of  the  energy  gap  between  these  states. 


*Thi8  work  is  also  supported  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AFOSR-81-0009-A  and  the  National  Science  Foundation  under  Grant 
NSF-CHE81-21945. 
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(5)  («)  Early  reporta  of  such  studies  for  fluorenylldene  have  appeared: 
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Hetherlngton,  G.  M.  Korenowski,  and  M.  J.  McAuliffe,  J.  Am.  Chem.  Soc. 
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B.  CHEMICAL  INTERMEDIATES  IN  THE  PHOTOCHEMISTRY  OF  AROMATIC  ENDOPEROXIDES 
SINGLET  OXYGEN* 

(C.  Dupuy,  S.  Y.  Hou,  E.  Sitzmann,  M.  McAullffe,  K.  B.  Elsenthal) 

(JSEP  work  unit  8,  1982-1985) 

(Principal  Investigator:  K.  B.  Elsenthal  (212)  280-3175) 


In  our  recent  studies  ^  of  the  photodisaodation  of  oxygen  containing 
(2) 

molecules,  and  in  particular  anthracene  endoperoxide ,  we  discovered  a  new 
dissociative  pathway  which  produced  not  only  oxygen  in  an  excited  state  but 
also  an  excited  anthracene  fragment.  A* 


Ph 


Ph  CHj 


To  determine  the  importance  of  this  pathway  for  bond  breaking,  we  measured  the 
formation  time  of  the  ground  state  anthracene  fragment,  A.  If  the  principal 
route  is  via  the  excited  state,  i.e.. 


AO ^  :A*  +  ^2 

+ 


1 


A 


then  *A  would  appear  as  the  excited  state  ^A*  decayed.  We  found  that  this  was 
not  the  case.  The  ground  state  appearance  time  was  found  to  be  45ps,  a  much 
longer  time  than  the  appearance  time  of  *A*  which  was  less  than  5ps.  This 
result  suggested  that  a  principal  path  for  bond  rupture  involves  not  excited 
anthracene,  but  a  two-step  process  involving  a  chemical  intermediate.  A 
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sequence  showing  two  possible  Intermediates  Is: 


t 


t 


We  ere  presently  seeking  to  trap  end  Identify  the  chemical 


Intermediate  in 


this  adiabatic  photodlssoclatlon. 


*Thls  work  Is  also  supported  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AFOSR-81-0009-A  and  the  National  Science  Foundation  under  Grant 
NSF-CHE81-21945. 
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C.  PICOSECOND  DYNAMICS  OF  TWISTED  INTERNAL  CHARGE  TRANSFER  PHENOMENA* 

(Y.  Wang,  M.  Crawford,  K.  B.  Elsenthal) 

(JSEP  work  unit  8,  1982-1985) 

(Principal  Investigator:  K.  B.  Elsenthal  (212)  280-3175) 

The  dual  fluorescences  phenomena  of  DMABN  and  related  molecules  la  a 

subject  that  has  received  intensive  studies  in  the  past/*^  ^  In  nonpolar 

solvents,  excitation  of  the  DMABN  molecule  results  in  a  single  fluorescence 

band  centered  around  340  nm.  However,  in  polar  solvents  a  new,  Stokes-ahifted 

fluorescence  appears,  which  is  associated  with  a  large  dipole  moment  of  -16 
(4) 

D.  To  explain  the  origin  of  this  new  emission  band,  a  number  of  mechanisms 
have  been  proposed.  From  the  results  accumulated  over  the  past  several  years, 
we  conclude  that  both  the  bond  twisting  (between  the  amino  group  and  the 
benzene  ring)  and  solvation  play  important  roles  in  causing  the  charge 
redistribution^  ^ 


In  spite  of  the  numerous  studies  on  the  subject,  the  relationship  between 
the  two  fluorescence  bands  was,  surprisingly,  not  fully  understood  until 
recently.  By  using  picosecond  fluorescence  spectroscopy/^  we  have  directly 
demonstrated  the  parent-daughter  relationship  between  the  short  wavelength 
emission  band  and  the  long  wavelength  emission  band.  We  have  shown  that  the 
two  emitting  species  reach  equilibrium  very  rapidly  (~20  pa)  with  an 
equilibrium  constant  30  in  propanol  solution /^ ^ 

To  study  the  role  of  solvent  molecules  in  the  intramolecular  CT  process 
of  excited  DMABN  molecules  and  hopefully  to  gain  a  better  understanding  of  the 
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coupling  between  the  twisting  notion  end  solvent  relaxation,  the  dynamics  of 
the  CT  transfer  process  in  nixed  solvents  were  investigated.  The  basic  idea 
is  to  control  the  local  nunber  density  of  the  polar  solvent  nolecules  without 
greatly  altering  the  viscosity  of  the  solution.  The  effect  of  the  change  in 
local  dipole  density  on  the  dynamics  and  yield  of  the  CT  process  could  then 
yield  information  on  the  role  of  the  polar  solvent  nolecules.  A  convenient 
way  to  adjust  the  local  dipole  density  is  by  diluting  the  polar  solvent 
nolecules  with  nonpolar  nolecules.  We  have  chosen  to  study  the  nixed  solvent 
system  butanol /hexadecane  because  of  their  similar  viscosities  3.0  and  3.3  cp 
at  20°C,  respectively. 

The  evolution  of  the  photoexclted  DMABN  molecule  from  the  initial  planar 
structure  to  the  final  perpendicular  structure  was  followed  by  monitoring  both 
the  singlet  fluorescence  decay  at  350  nm  and  the  CT  fluorescence  rise  at  470 
nm  as  the  butanol  solution  is  successively  diluted  by  the  nonpolar  hexadecane 
molecules . 

The  key  observations  made  were:  (1)  the  appearance  of  a  fast  and  a  slow 
kinetic  process  upon  the  dilution  of  butanol  solutions,  (2)  the  dynamics  of 
the  fast  process  is  Independent  of  the  dilution  while  the  rate  of  the  slow 
process  shows  a  linear  dependence  on  the  butanol  concentration,  and  (3)  the 
relative  Importance  of  the  two  processes  is  affected  by  the  dilution,  the  fast 
process  dominates  at  high  butanol  concentration. 

The  origin  of  the  fast  process  (i.e.,  the  30  ps  rise  observed  for  all 
butanol  concentrations),  we  believe,  is  due  to  the  excitation  of  the  complexes 
already  formed  in  the  ground  state  between  DMABN  and  butanol  molecules.  This 
30ps  rise  is  faster  and  larger  than  can  be  accounted  for  by  a  diffusion 
controlled  process  for  the  lower  butanol  concentrations.  The  formation  of  the 
ground  state  complexes  is  really  not  surprising  in  light  of  the  large  dipole 


According  to  our 


moment  7D  of  the  ground  state  DMABN  molecule. ^ 
interpretation,  the  measured  30  ps  risetime  for  the  fast  process  is  then 
mainly  due  to  the  Intramolecular  bond  twisting  motion  along  the  C-N  bond  of 
the  excited  DMABN-butanol  complex.  This  Intramolecular  twisting  motion  can  be 
impeded  by  the  frictional  force  of  the  solvent  molecules,  thus  accounting  for 
the  observed  viscosity  dependence. 

As  che  rapid  rise  component  (30  ps)  reflects  the  population  of  ground 
state  complexes,  the  slow  component  of  the  rise  reflects  those  excited  DMABN 
molecules  which  either  do  not  have  neighboring  butanol  molecules  or  do  not 
have  the  proper  arrangement  of  butanol  molecules  for  reaction  to  occur.  Even 
at  the  higher  alcohol  concentrations,  for  which  there  are  several  alcohol 
molecules,  on  the  average,  surrounding  a  ground  state  DMABN  molecule,  a  slow 
component  is  observed.  It  is  not  until  the  appropriate  configuration  is 
achieved,  as  given  by  the  slower  time  component,  that  the  rapid  30  ps  reaction 
can  occur.  The  key  feature  in  the  stabilization  of  the  twisted  intramolecular 
charge  transfer  structure  appears  to  be  the  prior  formation  of  a  complex 
between  excited  DMABN  and  some  butanol  molecule  or  molecules. 

Our  proposed  mechanism  can  be  briefly  summarized  in  the  figure  shown 
below,  noting  that  the  single  butanol  entity  shown  in  the  figure  can  be 
monomeric  or  part  of  some  higher  aggregate. 
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VI.  GENERATION  AND  CONTROL  OF  RADIATION 
A.  OPTICAL  COHERENT  TRANSIENT  SPECTROSCOPY* 


(R.  Beach,  B.  Brody,  T.  J.  Chen,  D.  DeBeer,  M.  Click,  R.  Klnchlnski, 
E.  Xu,  S.  R.  Hartmann) 

(JSEP  work  unit  9,  1982-1985) 

(Principal  Inveatigator:  S.  R.  Hartmann  (212)  280-3272) 


The  channeling  of  our  two  YAG  laser /amplifier  outputs  into  four  separate 
laboratories  has  been  completed  and  we  can  now  pump  dye  lasers  simultaneously 
at  all  our  experimental  stations.  A  centralized  computer  facility  has  been 
established  for  control  and  analysis  of  all  experiments.  Each  station  has  its 
own  computer  terminal  to  allow  for  local  interaction.  Analysis  is  performed 
while  data  is  acquired  and  the  results  are  displayed  in  graphical  form  at  the 
local  terminals. 

Photon  echo  signals  have  been  obtained  in  Li  vapor  over  an  unprecedented 
dynamic  range  exceeding  1010.^  Echoes  were  obtained  from  atoms  which  were 
excited  more  than  15  natural  lifetimes  earlier.  These  results  are  shown  in 
Fig.  1.  It  should  be  noted  that  our  largest  echoes  were  only  a  factor  of  25 
weaker  than  the  ir/2  excitation  pulse  which  first  excited  them.  The  purpose  of 
this  experiment  was  to  study  echo  degredation  due  to  nonparallel  excitation 
beams.  The  two  sets  of  data  in  Fig.  1  correspond  to  excitation  beam  angling 
of  0.70  and  1.40  mrad.  The  effect  of  beam  angling  above  is  obtained  by  taking 
the  ratio  of  these  data  and  is  presented  in  Fig.  2.  By  taking  the  ratio  we 
get  rid  of  radiative  lifetime  and  hyperflne  structure  effects.  For  beam 
angling  of  9,  the  echo  intensity  when  the  excitation  pulses  are  separated  by  a 
time  x  is  just 


echo 


'hkqQ0  2 

-2(— — )t  -  =— 

*  I0  e  e  T1  M(t) 


where  T^  is  the  natural  lifetime,  k  is  the  wavevector  associated  with  the 
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TIME  BETWEEN  1st  EXCITATION  (ir/2)  PULSE  AND 
ECHO  IN  NATURAL  LIFETIMES  (27  nsec)  OF  THE 
22P 3/2  STATE  OF  Li. 


Fig.  1  Photon  echo  intensity  vs  time  between  first 

excitation  pulse  and  echo  for  two  values  of  beam  angling. 
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ECHO  AMPLITUDE  {0:l.4Omrod)/ECHO  AMPLITUDE  (ds.70mrod) 


ANO  ECHO  IN  NATURAL  L IFE  TIMES  (27ns«c> 
OF  THE  22P3/zSTATE  OF  Li 


Fig.  2  Photon  echo  intensity  ratio  corresponding  to  9  *  14 

and  6  ■  .7  vs  time  between  first  excitation  pulse  and  echo. 
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optical  transition,  qg  ^ —  la  the  Inverse  thermal  de  Broglie  wavelength 

h 

and  M(t)  describes  the  aodulation  caused  by  the  hyperfine  structure.  The  ratio 


R(t) 


Iecho(T>  lA  arad) 

WT’  °*7  mrad) 


R- 


is  Independent  of  Tj  and  M(t).  For  our  experimental  conditions  -  55 

nsec.  The  solid  line  in  Fig.  2  uses  this  value  for 

The  Initial  motivation  of  this  experiment  was  provided  by  our  extension 

( y ) 

of  the  Billiard  Ball  Echo  Theory  to  Include  the  effects  of  using  sub- 

f  3\ 

Doppler-linewldth  lasers  in  photon  echo  experiments .  Sub-Doppler-llnewldth 
lasers  generate  nonspherical  excited  state  wave  packets,  and  photon  echo 
experiments  such  as  described  above  measure  the  transverse  wave  packet  shape. 
As  explained  in  the  last  progress  report,  the  Billiard  Ball  echo  model  begins 
by  localizing  gaseous  atoms  with  atomic  wave  packets.  When  these  wave  packets 
are  Irradiated  by  a  series  of  light  pulses  they  split  into  separate  parts 
corresponding  to  each  of  the  atomic  states  generated.  Each  resulting  wave 
packet  recoils  according  to  whether  it  absorbs  or  emits  a  photon  while  it  is 
being  formed.  Whenever  overlap  occurs  between  two  daughter  wave  packets  of 
the  same  parent,  coherent  radiation  ensues  and  photon  echoes  are  formed. 

The  Billiard  Ball  echo  model  was  developed  to  obtain  a  simple  analytical 
tool  which  would  work  equally  well  for  the  two  pulse  photon  echo  and  multiple 
pulse  -  multiple  level  echoes.  It  was  specialized  to  short  pulse  excitations 
which  uniformly  excited  all  atoms  independent  of  their  velocity  thus  giving 
rise  to  spherical  wave  packets  which  were  called  Billiard  Balls.  Long  pulse 
excitation  leads  to  a  more  complicated  shape;  however,  just  as  in  the  case 
with  the  spherical  Billiard  Balls,  wave  packet  spreading  can  be  ignored.  We 
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have  demonstrated  that  long  pulse  excitations  do  not  change  the  essential 

(3) 

features  of  the  Billiard  Ball  echo  theory.  Even  though  the  shape  of  the 
excited  state  wave  packets  Is  no  longer  spherical,  the  recoil  or  billiard  ball 
analysis  remains  the  same. 

In  the  most  general  case,  the  excited  vave  packet  shape  can  be  quite 
complicated.  However,  in  all  cases  it  is  characterized  by  a  length  to  width 
ratio  L/U  of 


L/W-Ti/./T,, 


for 


tD  <  T  1/e 


where  Tjye  characterizes  the  width  of  the  excitation  pulse  electric  field 
amplitude 


E-E  exp  {-l/2(t/xw  )2} 
o  1/  e 

whereas  is  a  measure  of  the  width  of  the  thermal  atomic  velocity 
distribution. 


•fik 

m  tD 


It  is  the  time  it  would  take  the  excited  state  wave  packets,  recoiling  with 
velocity  iik/m,  to  separate  from  the  ground  state  wave  packet  if  their  spatial 
extent  were  given  by  the  thermal  de  Broglie  wavelength. 

For  short  pulse  excitation  where  >  x ^  »  L/W  -  1  and  all  wave  packets 
are  spherical.  For  but  null  area  excitation  pulses,  the 

excited  wave  packet  shape  simplifies  and  becomes  ellipsoidal. 

The  longitudinal  extent  of  the  excited  wave  packet  can  be  easily  measured 
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as  It  is  equal  to  the  duration  of  the  photon  echo.  Our  echoes  have  in  fact 
always  been  as  long  as  our  excitation  pulses,  in  agreenent  with  our  calcula¬ 
tion.  The  transverse  dimension  of  the  wave  packet  can  be  measured  by  angling 
the  excitation  beams  as  then  the  wave  packet  collision  is  glancing  rather  than 
head  on.  This  is  clear  from  the  recoil  diagram  of  Fig.  3.  Our  experiments  in 

Li  vapor  Figs.  1  and  2  show  good  agreement  with  theory.  Our  experimental  and 

(4} 

theoretical  work  has  recently  been  summarised 

It  should  be  noted  that  analysis  of  photon  echoes  generated  by  angled 
excitation  beams  allows  one  to  obtain  a  measure  of  transient  atomic  velocity 
distributions.  This  may  prove  valuable  in  studies  involving  nonresonant 
photon  assisted  collisions  which,  depending  on  which  side  of  a  resonance  a 
laser  is  tuned,  excited  state  thermal  distributions  are  produced  which  are 
either  hotter  or  colder  than  that  associated  with  the  ground  state  atoms. 

In  setting  up  the  experiment  to  generate  photon  echoes  using  angled 
excitation  beams,  it  was  decided  to  exploit  the  superradiant  character  of  the 
photon  echo  to  obtain  signals  over  a  large  dynamic  range.  One  advantage  of 
using  angled  excitation  beams  is  that  the  echo  is  radiated  in  the  direction 
2k2~itj  where  and  k^  are  the  wave  vectors  of  the  first  and  second  excitation 
pulses  respectively.  Thus  the  echo  is  spatially  separated  from  the  excitation 
beams*  At  short  pulse  separation  times  the  echo  was  several  orders  of  magni¬ 
tude  greater  than  any  leakage  from  the  excitation  pulses  and  optical  shutters 
were  unnecessary.  At  long  separation  times  an  AOM  was  used  to  eliminate  stray 
leakage  so  that  our  photo-tube  response  was  quiet  enough  to  detect  .02  photons 
per  echo.  The  experimental  set  up  is  shown  in  Fig.  4. 

The  advantage  of  being  able  to  work  with  a  large  dynamic  range  is  that  it 
allows  precise  studies  of  relaxation.  Spectroscopic  measurements  are  also 
facilitated  as  the  limitations  imposed  by  the  uncertainty  principle,  although 
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Recoil  Diagram  for  nonparallel  pulse  excitation. 


Fig.  4  Schematic  diagram  of  photon  echo  beam 
angling  apparatus. 


still  present,  are  determined  by  the  total  time  an  atom  exists  in  its  excited 
state. 

The  limiting  value  of  the  dynamic  range  over  which  echoes  can  be  observed 
is  determined  by  the  competition  between  the  echo  signal  from  atoms  excited  by 
the  first  excitation  pulse  and  the  incoherent  fluorescence  signal  from  atoms 
excited  by  the  second  excitation  pulse  applied  a  time  t  after  the  first.  The 
former  decays  as  exp(-2T/T^)  the  latter  as  exp(-r/T^),  where  T^  is  the 
excited  state  lifetime.  Once  these  signals  become  equal,  the  signal  to  noise 
of  the  echo  varies  only  with  the  square  root  of  the  number  of  measurements. 

One  can  of  course  continue  to  average  but  the  payoff  is  greatly  diminished. 

We  estimate  this  condition  as  follows. 

The  number  of  photons  in  an  echo  ne  is  given  by^^ 

232  -2t/T 

n  =  (8/4;  n*  L)  (nAL)  Wt  e 

e  air  p 


where  n  is  the  number  density,  L  is  the  sample  length,  A  is  the  sample  area, 

.  *  * 

W  ■  1/T. ,  t  is  the  pulse  width,  and  @  ■  (T«/t  )  if  t  <  T„  otherwise  B  *  1. 

1  p  Z  p  p  2 

A"  J 

is  the  inhomogeneous  linewidth- 

2 

The  number  of  photons  incoherently  radiated  in  a  solid  angle  X  /A  during 

a  time  t  at  the  time  of  the  echo  is  given  by 
P 


n 


inc 


SnAL  t  Utp 


Equating  with  n^nc  and  relating  the  number  density  to  the  absorption 
coefficient  a  by  the  formula 


a 


* 
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yields  the  condition 


IT  Ti  A  T/Tl 

T  6  71  aL  *  e 

T2  X 

At  this  value  of  T  the  number  of  lifetimes  an  echo  atom  has  lived  is  given  by 


N  -  2  £n  g  4  aL} 
6  T*  X2 


If  we  work  with  Na  vapor  using  a  laser  such  as  ours  which  puts  out  5  nsec 
pulses  then 

t  ■  5  x  10  ^  sec  A  ~  3  cm2  0 

P 

A  *w  -in  **  -5  o 

T2  -  10  sec  1  ~  6  X  10  3  X 

Tj  ■  16  x  10  ^  sec  aL  ~  1 


and  we  obtain  N  -  42. 

Picosecond  excitation  of  long  lifetime  lnhomogeneously  broadened  lines 
can  increase  this  estimate  considerably.  He  are  now  setting  up  to  make  a 
serious  measurement  in  Na  vapor.  The  key  to  success  is  good  photon  counting 
technique  and  good  echo  collection  efficiency. 

Our  application  of  the  tri-level  echo  technique  to  relaxation  studies  in 
atomic  sodium  has  in  the  past  enabled  us  to  study  Noble— gas-induced  broadening 
of  transitions  to  Rydberg  S  and  D  states.^  He  are  now  in  the  process  of 
performing  tri-level  echo  experiments  in  Li  vapor  in  order  to  measure  the 
relaxation  of  S-S  and  S-D  superposition  states  due  to  collisions  with  various 
buffer  gases.  This  experiment  utilizes  three  dye  lasers.  One  laser  excites 
the  2S-2P  transition  at  t-0  while  a  second  laser,  also  at  t-0,  excites  a  2P-nO 
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or  2P-nS  transition.  At  some  delay  time  later  a  third  laser  excites  the  2P-nD 
or  2P-nS  transition.  These  experiments  will  be  able  to  measure  collision 
cross  sections  well  into  the  Rydberg  regime. 

A  second  kind  of  two-photon  tri  level  echo  is  also  being  investigated, 
this  time  in  Na  vapor.  We  first  resonantly  excite  the  3S-3P  transition  and 
then  at  a  later  time  we  two-photon  excite  the  3S-4D  transition.  We  observe  an 
echo  at  the  4D-3P  transition.  We  also  observe  a  coherent  burst  of  radiation 
at  the  4D-3P  transition  shortly  after  the  3S-4D  transition  is  excited.  This 
we  interpret  as  delayed  superradiance.  So  far  we  have  only  been  able  to 
generate  echoes  at  high  atomic  densities  where  the  absorption  length  of  the 
3S-3P  radiation  is  much  smaller  than  the  length  of  the  sample.  This  means 
that  the  4D-3P  echo  is  radiated  into  a  high  gain  medium.  As  a  consequence  it 
may  sometimes  be  radiated  early,  before  the  3D-3P  rephasing  is  complete.  We 
have  evidence  of  this  behavior  which  we  will  study  further.  We  are  now  in  the 
process  of  replacing  our  old  Hinch  type  lasers  with  our  newer  Thompkln's 
design  1GHz  dye  lasers.  This  will  allow  us  to  excite  4D-3P  superposition 
states  with  much  less  laser  power.  We  have  some  evidence  of  ionization 
generation  which  would  complicate  our  analysis.  This  will  be  much  less 
Important  with  our  new  lasers . 

Our  4D-3P  echoes  have  been  rather  weak.  Nevertheless  they  are  large 
enough  to  enable  us  to  make  colllslonal  cross  section  measurements.  We  find 
the  cross  section  for  4D-3P  superposition  state  dephasing  is  comparable  to 
that  for  3S-3P  dephasing.  By  using  our  improved  lasers  we  will  be  able  to 
obtain  a  precise  comparison. 

In  a  series  of  three  papers  written  in  1973  Aihara  and  Inaba  discussed 
the  general  problem  of  Multiple  pulse  excitation  of  three-level  systems. 

They  found  that  under  suitable  conditions  photon  echoes  were  generated  whose 


200 


character  differed  dramatically  from  those  previously  observed.  It  la 
generally  accepted  that  echoes  formed  In  a  two-level  system  appear  at  a  time 
after  an  excitation  pulse  equal  to  Its  separation  from  a  previous  excitation 
pulse  or  the  separation  of  two  previous  excitation  pulses.  This  now 
unremarkable  effect  Is  a  consequence  of  the  fact  that  two-level  atoms  In  a 
static  environment  have  resonance  frequencies  whose  fractional  deviation 
•  («i»o-Uj)/<uo  from  the  central  frequency  «o  of  the  transition  being 

(9) 

excited  after  each  excitation  pulse  is  a  constant.  '  In  other  words  two 

pulse  photon  echoes  generated  by  pulses  it  t  ■  0  and  t  •  t  form  at  t  ■  2t 

because  they  rephase  after  the  second  excitation  pulse  at  exactly  the  same 

rate  at  which  they  dephased  after  the  first  excitation  pulse.  Echoes  do  not 

form  at  all  If  Aw./w  changes  sign  after  an  excitation  pulse,  as  happens  In  F 
jo 

spin  echo  experiments  In  CaF^  when  the  second  pulse  has  area  ir.  On  the  other 
hand  echoes  generated  In  an  electric  quadrupole-broadened  magnetic  resonance 
experiment  can  appear  at  anomolous  delays  due  to  the  unequal  but  correlated 
line  broadening  of  the  energy  eigenlevels/10^ 

The  simplest  situation  In  which  one  obtains  high  correlation  In  the 
Inhomogeneous  broadening  of  electronic  energy  levels  Is  found  In  a  gas.  The 
source  of  Inhomogeneous  broadening  Is  the  Doppler  shift  which  scales  directly 
with  the  central  frequency  of  each  electronic  transition.  This  essentially 
perfect  correlation  has  made  possible  a  series  of  experiments,  done  in  this 
laboratory,  to  investigate  the  colllslonal  relaxation  of  states  of  similar 
parity. The  technique  is  also  applicable  to  studying  relaxation  of  states 
of  opposite  parity  as  we  are  now  doing  In  Na  vapor  by  observing  an  echo  on  the 
4D-3P  transition  generated  by  two-pulse  excitation  of  the  3S-3P  and  3S-4D 
transitions . 

Photon  echoes  In  solids  of  the  same  character  as  those  in  gases  should 
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also  be  possible.  Phonon  Induced  relaxation  however  puts  a  severe  limitation 

on  the  number  of  transitions  which  may  be  studied.  By  working  at  liquid  He 

temperatures  some  transitions  may  be  frozen  in  but  most  remain  unavailable  to 

an  echo  type  experiment  because  they  decay  too  fast  by  spontaneous  emission  of 
3+ 

a  phonon.  In  Pr  doped  LaP^  several  transitions  are  available  and  ordinary 

3  3 

two-pulse  photon  echoes  have  been  observed  on  both  the  -  Pq  and 
3  1 

transitions.  As  we  have  studied  echo  behavior  on  both  these 

transitions  it  would  seem  natural  for  us  to  extend  our  work  to  perform  a 

multiple-level  echo  experiment  in  this  medium.  This  could  be  accomplished  by 

3  1  3  3 

either  sequentially  exciting  the  -  D,,  and  -  Pq  transitions  and  then 

1  3 

looking  for  the  echo  formed  on  the  D2  -  Pq  transition  or  by  performing  an 

3  1 

ordinary  two-pulse  echo  experiment  on  the  -  B,,  transition  while 

3  3 

performing  a  single-pulse  excitation  of  the  -  Pq  transition.  This 

3  3 

single-pulse  excitation  of  the  -  Pq  transition  will  generate  an  echo 

3  1 

since  the  two-pulse  excitation  of  the  -  D2  transition  regenerates  a 

3 

state  component  whose  state  vectors  are  modulated  according  to  the 

Inhomogeneous  broadening  associated  with  the  transition.  It  is  best  to  think 

of  the  echo  formation  process  here  in  terms  of  a  stimulated  echo  where 

modulated  ground  state  components  are  recombined  after  a  single  pulse  to 

produce  an  echo.  Ordinarily  this  pulse  would  be  applied  at  the  same  frequency 

as  the  pulse  pair  which  produced  the  modulated  ground  state  components  but 

this  is  not  necessary.  The  only  requirement  is  that  the  excitation  pulse 

produce  a  superposition  state  which  rephases  at  some  rate  proportional  to  the 

dephaeing  rate  of  the  primary  transition  and  this  will  happen  according  to  the 

degree  to  which  the  Inhomogeneous  broadening  is  correlated. 

Of  the  two  experiments  discussed  above  the  first  is  more  difficult  for 
1  3 

three  reasons.  (1)  The  02  -  Pq  transition  is  in  the  infrared  where 
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sensitive  detectors  are  lacking.  (2)  The  transition  does  not  appear  to  be 
strong  enough  even  to  deserve  mention  in  the  literature.  (3)  Phase  matching 
may  be  a  problem  as  the  k  vector  induced  in  the  -  ^Pq  transition  may  not 
correspond  to  2 w/ A  where  A  is  the  wavelength  associated  with  this  transition. 
The  single  positive  feature  is  that  the  echo  forms  at  a  frequency  different 
from  either  excitation  pulse  so  that  excitation  pulse  rejection  by  the 
detection  system  would  be  facilitated.  This  latter  feature  together  with  the 
fact  that  only  two  lasers  are  required  prompted  us  try  this  experiment.  We 
failed  for  some  combination  of  the  first  three  stated  objectives. 

We  plan  therefore  to  work  with  the  second  scheme.  Phase  matching  is 
automatic  as  we  observe  the  same  wavelength  at  which  we  excite.  Excitation 
pulse  rejection  is  accomplished,  rather,  by  angling  the  excitation  beams. 

Calling  kj,  k^ ,  and  k^  the  wave  vectors  associated  with  the  first  and  second 

3  1  3  3 

pulses  on  the  -  D2  transition  and  the  pulse  on  the  -  PQ  transition, 

respectively,  we  make  f  (k^,  k2>  *  Xk^,  k2>  *  J  (k2>  k^).  The  ordinary 

two-pulse  echo  at  3H^  -  1D2  will  form  along  2k2  -  kj  while  our  new  tri-level 
3  3 

echo  at  -  PQ  will  form  along  ^3  ^2  ”  ^1*  t',e  echo*8  *re  separated 

spatially  from  each  other  by  the  same  amount  as  they  are  separated  spatially 

from  the  immediately  preceding  excitation  pulse.  Our  experience  in  observing 

3  1  3  3 

ordinary  echoes  on  both  the  -  D2  and  -  PQ  transitions  assures  us 
that  no  problem  will  be  encountered  in  isolating  the  echo.  This  is  important 
for  this  experiment  as  we  cannot  be  certain  that  our  Pockels  cell  shutters 
will  operate  quickly  enough  to  separate  the  echo.  Actually,  the  question  of 
being  able  to  see  an  echo  is  moot  as  in  the  worst  case  we  still  would  be  able 
to  follow  the  degradation  of  the  four-wave  mixing  signal  at  k^  +  kj  -  k^ 
formed  by  the  instantaneous  application  of  the  three  excitation  pulses.  Echo 
generation,  after  all,  is  no  other  than  a  form  of  delayed  four-wave  mixing. 
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The  question  of  whether  rephaaing  occurs  does  not  arise  when  the  pulses  are 
siaultaneous . 

The  degradation  of  the  four-wave-nixed  signal  with  pulse  separation 
should  nanifest  itself  in  a  change  of  echo  anplitude  and  shape.  Echo  position 
is  related  to  echo  shape  and  yields,  as  in  tri-level  echo  experiments  in 
gases,  the  scaling  factor  associated  with  any  correlated  inhomogeneous 
broadening.  Ve  also  expect  that  signal  modulation  will  be  present  which  would 
allow  an  alternative  means  of  studying  nuclear  hyperfine  structures. 

Observation  of  the  position  and  shape  of  echoes  generated  as  in  the  above 
scheme  will  provide  information  regarding  atoms  whose  transitions  are 
broadened  in  a  correlated  way. 

*Thls  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-DMR80-06966  and  the  U.  S.  Office  of  Naval  Research  under  Contract 
No.  N00014-78-C-0517 . 
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SIGNIFICANT  ACCOMPLISHMENTS  AND  TECHNOLOGY  TRANSITION  REPORT 


Significant  Accomplishments 

I.  Quant an  Detection  end  Sensing  of  Radiation 

Ve  have  provided  a  new  underatandlng  of  the  way  In  which  Individual 
atonic  photon  emissions  combine  together  to  generate  chaotic  light  at  high 
average  levels.  The  previous  theory  was  phenonenologlcal . 

We  have  solved  an  Important  branching  process  In  statistics,  a  cascade  of 
Poissons,  useful  In  applications  from  avalanche  photodiodes  to  cosnlc  ray 
nolae. 

We  have  defined  the  role  of  signal-dependent  noise  In  VLSI  system 
reliability. 


II.  Physical  and  Photochemical  Properties  of  Electronic  Materials 

For  the  first  time  high  quality  germanium  nitride  Insulating  films  have 
been  grown  successfully.  It  now  appears  that  a  germanium  insulated  gate 
field-effect  transistor  can  be  fabricated  using  the  technology  developed  In 
our  laboratory. 

A  new  tunneling  rectifier  has  been  fabricated  using  a  metal-amorphous 
silicon-silicon  sandwich  structure. 

For  the  first  time,  grain  boundary  parameters  have  been  measured 
unambiguously  using  a  laser-spot  photoconductivity  technique  developed  In  our 
laboratory.  These  parameters  Include  barrier  height,  interface  state  density, 
energy  level  and  capture  cross  section. 

Models  Including  both  Coulomblc  scattering  and  phonon-assisted  processes 
at  a  semiconductor-semiconductor  Interface  have  been  completed. 

The  ultraviolet  spectra  of  physlsorbed  metallorganlcs  has  been  measured 
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for  che  first  time.  Metellorgenics  are  a  prevalent  metal-carrier  species  in 
many  phases  of  microelectronics  fabrication,  including  laser  direct-writing. 
Further,  physlsorbed  layers  of  these  molecules  are  believed  to  play  an 
Important  role  in  many  phases  of  microelectronics  processing. 

The  first  observation  of  laser-enhanced  surface  chemistry  has  been  made. 
In  these  experiments,  we  have  found  that  submicrometer  metal  spheres  can  be 
used  to  enhance  the  surface  photochemical  decomposition  rate  of  adsorbed 
molecules.  By  arranging  these  spheres  in  an  organized  pattern  on  the  surface, 
it  may  be  possible  to  further  enhance  the  surface  reaction  rate.  Thus 
photolithographic  fabrication  techniques  which  have  been  developed  for 
fabrication  of  semiconductor  integrated  circuits  may  prove  useful  for 
producing  catalytic  surfaces  for  photochemical  reactions. 

The  direct  fabrication  of  submicrometer  diffraction  gratings  has  been 
achieved  using  laser  photochemical  Interferometry.  These  results  have  shown 
that  high  quality  optical  gratings  with  periodicities  as  small  as  130  nm  can 
be  made  readily  with  this  laser  technique.  The  "conventional*’  approach  uses 
electron  beam  lithography  and  is  thus  slower  and  much  more  expensive.  Five 
period  gratings  on  compound  semiconductors  are  Important  as  components  for 
electroptical  devices  and  integrated  optical  circuits. 

III.  Characteristics  and  Properties  of  Energetic  Ion  and  Radical  Beams 

The  energy  distribution  of  both  ions  and  neutrals  emanating  from  a 
caplllarltron  ion  source  have  been  determined  using  time  of  flight  mass 
spectrometry.  The  neutrals  have  been  found  to  be  extraordinarily  cold  while 
the  ions  are  quite  hot. 
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IV.  Energy  Transfer  Processes  in  Quantum  Electronics 

Collisional  redistribution  of  electronic  energy  in  sulphur  dioxide  has 

» 

been  found  to  preferentially  excite  the  symmetric  stretching  mode  of  the 
ground  electronic  state.  The  antisymmetric  stretching  mode  appears  to  be 

filled  only  indirectly  by  energy  transfer  from  the  symmetric  stretch  and  not 

t 

directly  from  the  electronically  excited  state. 

Photochemically  generated  hot  hydrogen  atoms  have  been  found  to  produce 
bend-stretch  excitation  preferentially  over  pure  stretch  excitation  in  CO. 
during  H/CO^  collisions.  In  H/CO  collisions  large  amounts  of  rotational 
excitation  are  produced  in  CO  along  with  vibrational  excitation  of  the 
molecule. 

3  * 

The  lifetime  of  the  I^  n2u^  ^  state  has  been  measured  for  the  first 
time.  A  matrix  Isolation  technique  (10°K),  which  takes  advantage  of  the  cage 

effect  to  produce  these  triplets  following  laser  excitation,  has  been  used  in 

f 

conjunction  with  IR  fluorescence  detection  to  determine  the  lifetime  of  this 

level.  This  state  may  be  Involved  in  limiting  the  maximum  power  output  of  the 

chemical  I  atom  laser, 
t 

V.  Picosecond  Energy  Transfer  and  Photofragmentation  Spectroscopy 

Ve  have  established  that  the  stabilization  of  Intramolecular  charge 
transfer,  which  is  the  origin  of  the  anomalous  fluorescence  in  dimethyl- 
amlnobenzonitrile,  is  due  to  formation  of  an  electrostatic  ground  state 
complex  with  a  polar  solvent  molecule  or  molecules. 

We  have  carried  out  the  first  measurement  of  the  equilibration  time  for 
spin  conversion  in  a  carbene.  The  value  for  dl phenyl carbene  is  llOps.  From 
the  kinetic  measurements  we  also  obtained  an  energy  splitting  of  3.9  kcal/mole 
between  the  ground  triplet  state  and  the  neighboring  singlet  state. 
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VI.  Generation  and  Control  of  Radiation 


Photon  echoes  have  been  generated  In  atoms  excited  more  than  IS  lifetimes 
earlier.  These  experiments  were  performed  in  Li  vapor  where  it  was  possible 
to  observe  echo  signals  over  an  unprecedented  dynamic  range  exceeding  lO1**. 


Technology  Transition  Report 

1.  Bell  Laboratories  has  initiated  several  programs  in  laser  chemical 
processing  such  as  those  being  studied  at  Columbia.  As  a  result  of  this 
interest  they  have  recently  awarded  a  Bell  Laboratory  Graduate  Fellowship  to 
Esaul  Sanchez.  His  work  will  stress  the  use  of  surface  spectroscopy  In 
diagnosing  surface  reactions. 

2.  Allied  Chemical  has  recently  become  active  in  developing  new  electronic 
chemicals.  They  have  had  a  significant  interest  in  finding  new  applications 
based  on  laser  processing.  With  this  in  mind,  they  have  assisted  members  of 
Prof.  Osgood's  group  in  developing  a  new  technique  for  writing  micrometer 
scale  oxide  patterns  and  silicon  wafers.  The  technique,  which  uses  a  spin-on 
silicate  material,  is  expected  to  be  useful  for  directly  writing  custom  rib 
waveguides . 

3.  Although  too  numerous  to  mention,  there  are  many  examples  of  laser 
applications  to  defense,  industrial,  and  research  uses.  One  of  the  most 
potentially  interesting  is  the  use  of  lasers  in  photochemical  direct  writing 
of  microcircuits.  The  original  research  leading  to  the  development  of  the 
laser  was  begun  under  the  JSEP  program  at  Columbia  in  1951  when  work  on  the 
maser  started  under  the  direction  of  Professor  Townes. 
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U.  S.  Army  Electronics  RAD  Command 
Night  Vision  and  Electro-Optics  Labs 
ATTN:  Dr.  Randy  Longshore,  DELNV-IT 
Fort  Bel voir,  VA  22060 

Commander 

U.  S.  Army  Research  Office 
ATTN:  DRXRO-EL  (Dr.  James  Mink) 

P.  0.  Box  12211 

Research  Trlangls  Park,  NC  27709 
Commander 

Harry  Diamond  Laboratories 
ATTN:  DELKD-RT-A  (Nr.  J.  S.l.rno) 

2800  Povdar  Mill  Road 
Adelphl.  MD  20753 

Director 

U.  S.  Army  Electronic,  RAD  Commend 
Night  Viiion  end  Electro-Optic.  Lab. 
ATTN-  DELNV-IRTD  (Dr.  John  Pollard) 
Port  Relvolr,  VA  22060 

Commander 

L'.  S.  Army  Reaeerch  Office 

ATTN:  DRXRO-EL  (Dr.  William  A.  Sander 

P.  O.  Box  12211 

Reaearch  Triangle  Park,  NC  27709 


Director 

U.  S.  Army  Electronlce  Technology 
and  Devices  Laboratory 
ATTN:  DELET-E3  (Dr.  A.  Tauber) 

Fort  Monmouth,  V  07703 

Director 

Division  of  Neuropsychiatry 
Walter  Read  Army  Inst,  of  Res. arch 
Washington.  DC  20012 

Cownndar 
USA  AMADCOH 

ATTN:  DEDAR-SCF-IO  (Dr.  J.  lavnda) 
Dover.  NJ  07*01 

Director 

U.  S.  Army  Signals  Warfare  Lab 
ATTN:  DELSW-D-OS 
Vint  Hill  Parma  Station 
Warranton,  VA  22166 

Director 

U.  S.  Army  Electronics  Technology 
and  Devices  Laboratory 
ATTN:  DEIET-ED  (Dr.  E.  H.  Poindexter) 
Fort  Moomouth,  NJ  07703 

Commander 

U.  S.  Army  Research  6  S t ends rd lee t ion 
Group  (Europe) 

ATTN:  (Dr.  f.  Rothwarf) 

Box  65 

FPO  NY  09510 

U.  S.  Army  Reaeerch  Office 
ATTN :  Library 
P.  0.  Box  12211 

Reaearch  Triangle  Perk,  NC  27709 
Cowendcr 

U.  S.  Army  Communications  Command 
ATTN:  DRSEL-COM-RF  (Dr.  T.  Rleln) 

Fort  Monmouth,  NJ  07703 

Mr.  Jerry  Brookshire 
Guidance  and  Control  Directorate 
U.  S.  Army  Missile  Cowand 
ATTN:  DRSMX-RGC ,  Bldg.  4381 
Radatonc  Arsenal,  AL  35898 

Dr.  Michael  Fahey 
Advanced  tensors  Directorate 
U.  S.  Army  Missile  Cowand 
ATTN:  DRDMI-RER 
Redstone  Arsenal,  AL  »  35898 

Dr.  Charles  Bowden 
U.  S.  Army  Missile  Command 
Research  Directorate 
ATTN:  DRSM1-RFD 
Redstone  Arsenal,  AL  35898 

Dr.  Arthur  R.  Slndoris 
Harry  Diamond  Laboratories 
ATTN:  DELKD-PO-P 
2800  Powder  Mill  Road 
Adelphi ,  MD  20783 

t 

Dr.  Horst  R.  Wittmann 
U.  S.  Armv  Research  Office 
P.  0.  Box  12211 

Reaeerch  Trijiu!** Park,  NC  27709 


Dr.  Jlalt  R.  $uccl« 

U.  S.  Anry  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709 

Hr.  Charles  Graff 
U.  S.  Army  Com.  -  Electronic# 
CoMand 

ATTN ;  DRSEL-COM-RF-Z 
Fort  Monmouth,  NJ  07703 

Mr.  Edward  Herr 

U.  S.  Arsy  Com.  -  Electronics 
Coaesnd 

ATTN;  DRSEL-COM-RX-4 
Fort  Monsouth.  NJ  07703 

Mr.  Roland  Wright 

Night  Vision  &  Electro-Optics 

Labs 

Fort  Belvoir,  VA  22060 

Dr .  Robert  Rohde 

Night  Vision  4  Electro-Optics 

Labs 

Fort  Belvoir,  VA  22060 

Dr.  Donn  V.  Campbell 

V.  b.  Armv  Com.  -  Electronics 
CotMwnd 

ATTN:  DRSEl-COM-RN-4 
Fort  .Monmouth,  NJ  07703 

Dr.  Nick  Karaylanls 
Harry  Diamond  Laboratories 
ATTN:  DEUCD-RT-CA 
2800  Powder  Hill  Road 
Adelphi.  MD  20783 

Dr.  T.  N.  Chin 
U.  S.  Arsy  ARRADC0M 
ATTN:  DRDAR-SCF-10 
Dover,  NJ  07801 

Dr.  John  Malatsus 

Night  Vision  6  Electro-Optics 

Labs 

Fort  Belvoir,  VA  22060 

Dr.  Rudolf  G.  Burner 
Night  Vision  6  Electro-Optics 
Labs 

ATTN :  DELNL-L 

Fort  Belvoir,  NJ  22060 

Dr.  W.  Ealy 

Night  Vision  6  Electro-Opc ics 

Labs 

ATTN:  DELNV-AC 

Fort  Monsouth,  NJ  22060 

Dr.  J.  Hall 

Night  Vision  6  Electro-Optics 
Labs 

ATTN:  DELNV-AC 
Fort  Belvoir,  NJ  22060 

Dr.  J.  Burgess 

Night  Vision  6  Electro-Optics 

Labs 

ATTN:  DELNV-RH-RA 
Fort  Belvoir,  NJ  22060 


DEPARKENT  OF  THE  AIR  FORCE 
Dr.  E.  Champagne 

afwal/aadd-i 

Wright-Pat tcreon  AFB,  OH  45433 

Mr.  V.  Edwards,  Chief 
AFWAL/AAD 

Wright-Pat terson  APB,  OH  45433 


Profesaor  R.  E.  Fontana 
Head,  Department  of  Electrical 
Engineering 
AF1T/ENG 

Wright-Pat terson  AFB,  OH  45433 

Dr.  Alan  Garscadden 
AFVAL/POOC-3 

Air  Force  Aeronautical  Labs 
Wright -Pat terson  APB,  OH  45433 

Mr.  Alan  R.  Barnum  (CO) 

Rose  Air  Development  Center 
Grlffiss  AFB,  NY  13441 

Chief,  Electronic  Research  Branch 
AFWAL/AADR 

Wright-Pat tertoo  AFB,  OH  45433 

Mr.  John  Hott-Smlth  (ESD/ECE) 

HQ  ESD  (AFSC),  Stop  36 
Hanscom  AFB ,  KA  01731 

Dr.  J,  Ryles 
Chief  Scientist 
AFWAL/AS 

Wrlght-Pattcrson  AFB,  OH  45433 
Dr.  Allan  Schell 

radc/ee 

Hanscom  AFB.  MA  01731 

Dr.  J.  Bran 
AFOSR/NN 

Bolling  AFB,  DC  20332 

Dr.  David  W.  Fox 
AFOSR/NN 

Bolling  AFB.  DC  20332 

Dr.  J.  Neff 
AF0SR/NE 

Bolling  AFB,  DC  20332 

Dr.  N.  H.  DeAngclis 
RADC/ESR 

Hanscom  AFB,  HA  01731 


Dr.  Gerald  L.  Witt 
Program  Manager 
Electronic  6  Material 
Sciences  Directorate 
Department  of  the  Air  Force 
AFOSR 

Bolling  AFB,  DC  20332 


Dr.  Tom  Welsh 
AF0SR/NE 

Bolling  AFB,  DC  20332 

Dr.  Edvard  Altshuler 
RADC/EE? 

Hanscom  AFB.  MA  01731 
department  or  the  navy 


Navel  Surface  Weapons  Center 
ATTN-  Technical  Library 
Code  0X-21 
Dahlgren,  VA  22448 

Dr.  Gernot  M.  R.  Winkler 
Director,  Tims  Service 
U.  S.  Navel  Observatory 
Massachusetts  Avenue  at 
34th  Street,  NW 
Washington,  DC  20390 

G.  C.  Dllvorth,  Jr. 

Technical  Director 

Navel  Coastal  Systems  Center 

Panama  City.  FL  32407 

Naval  Air  Development  Center 
ATTN:  Code  -  301  A.  Witt 
Technical  Library 
Warminster,  PA  18974 

R.  S.  Allgaler.  R-45 
Navel  Surface  Weapon#  Center 
Silver  Spring,  MD  20910 

Office  of  Naval  Research 
800  North  Quincy  Street 
ATTN:  Code  250 
Arlington,  VA  222 17 

Office  of  Naval  Research 
A00  North  Quincy  Street 
ATTN:  Code  414 
Arlington,  VA  22217 

Office  of  Navel  Research 
800  North  Quincy  Street 
ATTN:  Code  4 LIMA 

(Dr.  Stuart  L.  Brodsky) 
Arlington,  VA  22217 

CoMsndlng  Officer 
Naval  Research  Laboratory 
ATTN:  Dr.  S.  Teitler,  Code  1450 
Washington,  DC  20375 

CoManding  Officer 
Naval  Research  Laboratory 
ATTN:  Mrs.  D.  Folen,  Code  2627 
Washington,  DC  20375 

Commsndlnc  Officer 

Naval  Research  Laboratory 

ATTN:  Mr.  A.  Brodxinsky.  Code  5200 

Washington,  DC  20375 

Commanding  Officer 
Naval  Research  Laboratory 
ATTN:  Nr.  J.  E.  Davey,  Code  68 10 
Washington,  DC  20375 
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Coaaendlng  Officer 

Navel  Research  Laboratory 

ATTN:  Nr.  8.  0.  HeCoabe.  Code  6600 

Washington .  DC  20375 

Coaaand  lng  Officer 
Navel  Research  Laboratory 
ATTN:  Nr.  W.  1.  Feuec,  Code  650* 
Washington,  DC  20175 

Technical  Director 

Naval  Underwater  Systeas  Center 

New  London,  Cl  06220 

Navel  Research  Laboratory 
Underwater  Sound  Reference  Decachaenc 
Technical  Library 
p.  0.  Sox  8117 
Orlando,  FL  22856 

Naval  Ocean  Svstens  Center 
ATT*?  Dr.  F.  C.  Fletcher,  Code  92 
San  Diego.  CA  92152 

Naval  Ocean  Systaaa  Center 
ATTN:  Nr.  W.  J.  Delka,  Code  8302 
San  Diego,  CA  92152 

Naval  Ocean  Systeae  Center 
ATTN:  Dr.  Alfred  K.  Nedoluha, 

Code  922 
San  Diego,'  CA  92152 

Naval  Weapons  Center 

ATTN:  Dr.  C.  H.  Winkler.  Code  381 

China  Lake.  CA  93555 

Dr.  Donald  E.  Kirk  (62) 

Professor  end  Chelraen,  Electrical 
Engineering 
SF-104 

Naval  Postgraduate  School 
Monterey,  CA  939*0 

Dr.  D.  F.  Dence 

Naval  Underwater  Systeas  Center 
New  London  Laboratory 
ATTN:  Code  3* 

New  London,  CT  06320 

Director,  Technology  Assessaent 
Division  (OP-987) 

Office  of  the  Chief  of  Naval  Oper. 
Navy  Departaent 
Washington,  DC  20150 

Hr.  J.  W.  Willis 

Naval  Air  Systeas  Coses* nd 

AIR-310 

Washington,  DC  20361 

Naval  Electronics  Systeae  Coaaand 
MC  #  l 

ATTN:  Code  61R 

2511  Jefferson  Davis  Highway 

Arlington,  VA  20360 

Departaeac  of  the  Navy 

Naval  See  Systeas  Coaaand 

ATTN:  W.  W.  Elaine  (SEA-62R) 
Washington,  DC  20162 


David  Taylor  Naval  Ship  Research 
end  Developatmt  Canter 
ATTN:  Nr.  C.  H.  Glslssner,  Code  18 
Sathesda,  t®  2008* 

Hr.  Martin  Hendelberg 
Coast  Guard  RAD  Center 
Avery  Folnt 
Croton,  CT  06340 

Naval  Underwater  Systeae  Center 

New  London  Laboratory 

ATTN:  101R  (Dr.  Edward  S.  Eby) 

New  London,  CT  06320 

Nr.  Thoaee  J.  Manuccia,  Need 
Chealstry  end  Application  Section 
Coda  65*3 

Naval  lesesrch  Laboratory 
Washington,  DC  20375 

Dr.  Stephen  C.  Bishop,  Head 
Semiconductor  Branch 
Code  6870 

Navel  Reteerch  Laboratory 
Washington,  DC  20375 

Dr.  John  W.  Rockway 
Caaa.  Tachnology  Prog.  Off. 

Coda  8105 

Navel  Ocean  Systeas  Centar 
Sea  Diego,  CA  92132 

Dr.  Berry  F.  Shay 
Joint  Program  Office, 

ODUSD(P) 

The  Pentagon,  Ra  AD82S 
Washington,  DC  20301 

Dr.  Sydney  R.  Parker 
Professor,  Electrical  EnRlnaerlng 
Code  62FX 

Naval  Postgraduate  School 
Noneerey,  CA  939*0 

Dr.  CeorES  B.  Wright 
Office  of  Newel  Research 
Code  *27 

Arlington,  VA  22217 

OTHER  GOVERNMENT  AGENCIES 

Dr.  Ronald  E.  Reger tee 
Director 

Division  of  Materials  Research 
National  Sclance  Foundation 
1800  G  Straet 
Washington,  DC  20550 

Director 

Division  of  Electrical,  Coaputer 
end  Syeceas  Engineering 
National  Science  Foundation 
Weehington,  DC  20550 


Dr.  Dean  L.  Mitchell 
Section  Heed 

Condensed  Hatter  Sciences  Section 
Division  of  Materials  Research 
National  Science  Foundation 
1800  G  Street,  N.  W. 

Weehington,  DC  20550 

Judoon  C.  French,  Director 
Center  for  Electronics  and  Eleccrlcel 
Engineering 

8  358  Metrology  gulldlng 
National  Inreeu  of  standards 
Washington,  DC  20234 

NON -GOVERNMENT  AGENCIES 

Director 

Columbia  Radiation  Laboratory 
Columbia  University 
538  Waat  120th  Street 
New  York.  NY  10027 

Director 

Coordinated  Science  Laboratory 
University  of  Illinois 
Urbane,  1L  61801 

Associate  Director  pt  Materials 
end  Electronics  Research 
Division  of  Applied  Sciences 
McKay  Laboratory  107 
Harvard  Ualverslty 
Cambridge,  HA  02138 

Director 

Electronics  Research  Center 
University  of  Texes 
F.  0.  Box  7728 
Austin.  TR  78712 

Director 

Electronics  Research  Laboratory 
University  of  California 
Berkeley.  CA  9*720 

Director 

Electronics  Sclencee  Laboratory 
University  of  Southern  California 
Los  Angelas,  CA  90007 

Director 

Microwave  Research  Institute 
Polytechnic  Institute  of  New  York 
333  Jay  Street 
Brooklyn,  NY  11201 

Director 

Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Caabrldgo.  NA  02139 

Director 

Stanford  Electronics  Laboratory 
Stanford  University 
Stanford,  CA  9*305 

Director 

Edward  L.  Clnxton  Laboratory 
Stanford  University 
Stanford,  CA  9*305 
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Dr.  Letter  Eeetaen 

School  of  Electrical  Engineering 

Cornell  University 

116  Phillips  Hall 

Ithaca,  KT  14850 

Dr.  Carlton  Walter 
Electro  Science  Laboratory 
The  Ohio  Stete  University 
1320  Rlnnear  Road 
Coluabue,  OH  41212 

hr.  Richard  Sasha 
Dept,  of  Electrical  Engineering 
Texas  Tech  University 
Lubbock,  TX  79409 

Director 

School  of  Electrical  Engineering 
Ceorgla  Institute  of  Technology 
Atlanta,  CA  10132 

Dr.  John  F.  Walkup 
Dept,  of  Electrical  Engineering 
Texas  Tech  University 
Lubbock,  TX  79409 

Mrs.  Renats  D'Arcangelo 

Editorial  Office 

130  Pierce  Hall 

Division  of  Applied  Sciences 

11  Oxford  Street 

Caabrldge.  HA  02138 


